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Flooding and Drainage Effects on Slash 
Pine and Loblolly Pine Seedlings 


SLasH PINE (Pimus elliottii var. elliotti 
Engelm.) and loblolly pine (P. taeda L.) 
seedlings appear to have a high degree of 
tolerance to poor drainage and flooding. 
This tolerance is important in the silvicul- 
ture of the southeastern Coastal Plain flat- 
woods where imperfectly drained Bladen 
clay loam and related soils predominate. 
Approximately 75 percent of the area in 
these soils, in which both surface and in- 
ternal drainage are poor, is devoted to 
woodlands. Pines are the most valuable 
species, yielding high quality timber, pulp- 
wood, and naval stores. Generally, the 
pines are restricted to the better drained 
sites while less valuable hardwoods are the 
dominant species in low-lying, poorly 
drained zones. Throughout this region, ex- 
tensive areas occur in which, when drained, 
stands of pine develop from natural seed- 
fall. Both stocking and growth of pines 
may be substantially increased by providing 
adequate drainage. 

Trousdell and Hoover (1955) are 
mindful that “in years with heavy spring 
rainfall, it may be difficult to regenerate 
pine stands on some poorly drained sites 
by clearcutting systems because of stand- 
ing water on the prospective seedbed dur- 
ing the germination period.” Wet condi- 
tions also result in death of residual seed 
trees, 

In addition to silvicultural implications, 
high water causes costly and unpredictable 
logging operations. Heavy 
makes logging roads impassable, and the 
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puddling of soil by compaction reduces 
productivity of a site. 

No major research has been directed at 
determining effective and economical 
drainage coefficients for woodlands and, 
therefore, the full benefit of particular 
water removal rates to both tree growth 
and logging remain speculative. Agricul- 
tural and civil engineers requested to de- 
sign woodland drainage systems have little 
data on which such designs may be based. 
Nevertheless, there is much active interest 
in woodland drainage and engineers are 
receiving an increasing number of requests 
for assistance in the layout of water-re- 
moval systems. 

This paper presents the results of an 
effort to determine the effect of (1) con- 
tinuous flooding, (2) continuous drainage, 


L. C. Walker is Associate Professor, School 
of Forestry, University of Georgia; R. L. 
Green, now head of the Department of Agri- 
cultural Engineering, University of Mary- 
land, was Superintendent of the Southeastern 
Tidewater Experiment Station at Fleming, 
Georgia, during the course of this study. 
J. M. Daniels is an engineering aid at the 
station. Journal Series No. 119 of the College 
Experiment Station of the University of 
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Stations. The study, a cooperative project of 
the University and the Soil and Water Con- 
servation Research Division, Agricultural Re- 
search Service, U. S. Dept. Agric., was sup- 
ported in part by the Georgia Forest Research 
Council. Manuscript received Jan. 20, 1960. 
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and (3) three drainage coefficients upon 
survival and growth of slash and loblolly 
pine seedlings. A drainage coefficient is 
the amount of water, in inches, removed 
from a site during a 24-hour period. 


Literature Review 


The effect of deep channel drainage on 
slash and loblolly pines in the North Caro- 
lina coastal area, where ditching lowered 
water tables to a distance as great as 1000 
feet, has been reported by Pruitt (1947) 
and Schlaudt (1955). Height growth was 
increased from the ditch to a distance of 
500 feet, and soil physical properties were 
improved. The slightly detrimental growth 
effect of flooding upon shortleaf (Pinus 
echinata Mill. ), loblolly and pond pines (P. 
serotina Michx.) is reported by Hunt 
(1951). Although loblolly and shortleaf 
pines made more growth than pond pine 
in sand cultures flooded for various periods, 
in general, all seedlings proved unusually 
resistant to injury. After 12 weeks, slight- 
ly reduced growth was recorded for seed- 
lings continuously flooded with standing 
water. Flooding for three months did not 
permanently damage root systems; but 
after ten months, the roots appeared so 
badly injured that the plants would prob- 
ably have died if the soil dried below field 
capacity. Gaiser' also reported loblolly 
pine to be more resistant than pond pine 
to flooding injury. Both species made 
better growth in soil intermittently flooded 
for short periods, 2 days out of 9, than in 
soil maintained at field capacity. 

Prompt, removal of excess water from 
small stagnant bogs beneficially influenced 
leader growth of black spruce Picea mari- 
ana ( Mill.) B.S.P. and balsam fir (A dies 
balsamea L. Mill.) saplings. Acceleration 
of growth decreased with distance from 
the drainage ditch in this Michigan study 
(Satterlund and Graham 1957). 

The effect of flooding on seedling sur- 
vival is exhibited even for such hydrophytes 
as baldcypress, Taxodtum distichum (L.) 


'Gaiser, R. N. Unpublished thesis. Duke 
University, 1947. 


Rich. Ten to 12 days submergence may 
cause death. Demaree (1932) observed 
that this species must grow to sufficient 
height during the first year to stay above 
floods during the second year except for a 
very few days at a time. 

Other species detrimentally affected by 
high water include red pine (Pinus resinosa 
Ait.) (Stone et al. 1954), conifers of the 
Lake States’ region (Ahlgren and Han- 
sen 1957), and numerous hardwoods in 
the upper Mississippi River (Green 1947, 
Yeager 1949) and adjacent to northern 
Swiss lakes (Kuster 1948). Ahlgren and 
Hansen’s work showed no uniformity in 
the effect of duration of submergence on 
subsequent terminal growth. Trees flooded 
for less than 28 days in many cases showed 
as much growth reduction as trees flooded 
for longer periods. 

Redox potential measurements show that 
conditions favoring oxidation and/or re- 
duction in a soil in the absence of oxygen 
are not readily removed by additions of 
oxygen deficiency were considered to be 
waterlogging on plant growth may persist 
for some time after the soil is drained, and 
even temporary flooding may be injurious 
to plant growth (Scott and Evans 1955). 

The variation among species in ability 
to withstand flooding is, of course, ancient 
knowledge, but the mechanics of continued 
growth are speculative. Kramer (1949) 
and Leyton and Rousseau (1958) note 
that Salix roots, for instance, grow and 
absorb water in almost complete absence 
of oxygen. Heliotropium, on the other 
hand, ceased water absorption 15 hours 
after the soil atmosphere was replaced by 
nitrogen. Carbon dioxide toxicity and 
oxygen deficiency, were considered to be 
the causes. Reduced water absorption 
might be due to lessened metabolic activity 
of roots or physical changes in permeability. 
It is accompanied by retarded transpira- 
tion, photosynthesis, and mineral absorp- 
tion, and finally manifested as_ wilting. 
Kozlowski (1958) and Bergman (1920) 
reported that in flooded soils, absorption of 
water by roots often decreased and lagged 
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behind transpiration, causing leaves to dry 
out. 

Further evidence that absorption is not 
the only factor inducing injury was pre- 
sented by Kramer (1933) to the eftect 
that loblolly pine seedlings absorbed water 
and appeared uninjured for two weeks 
after roots were killed. As Kramer later 
stated (1951), lack of water cannot ex- 
plain the curvature of leaf tips, hypertrophy, 
nor development of adventitious roots. 
These changes take place more likely in 
turgid, growing tissues. He suggested that 
flooding injury may be related to disturb- 
ance of translocation of carbohydrates and 
possibly of hormones. Interference with 
downward translocation due to lack of 
oxygen “results in an accumulation of 
carbohydrates in the lower part of the 
stem, near the water line, causing hyper- 
trophy and the development of adventitious 
roots.” Kramer considers it unlikely that 
nutrient deficiencies are an important fac- 
tor in flooding injury but, rather, he sug- 
gests that poisoning by toxic substances 
moving up from dying roots may be the 
reason. Loblolly pine transpiration, ac- 
cording to Parker (1950), increased to 
125 percent of the expected rate for a 


week after flooding, then decreased con- 
stantly to 60 percent when drained 4 
weeks after flooding, then a week later to 


20 percent. 

Water-logged soils high in organic mat- 
ter, typical of many southeastern coastal 
forest sites, may also contain toxic quanti- 
ties of iron, sulphides, and manganese, 
built up as a result of the presence of COz 
released in biological activity (Kramer, 
1949). Chlorosis and wilt of leaves of 
flooded plants is known to decrease with 
removal of COz from soil solutions (Jack- 
son, 1956). 

Dean (1933) found aeration impor- 
tant in the production of root hairs and 
lateral roots, even for aquatic plants. In 
studying hydrophytes, Bergman (1920) 
reported that as soon as the plants grow- 
ing beneath the surface of swamp water 
showed signs of wilting, supplying air or 
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oxygen would result in’ restoration to 
normal.” 


Description of the Site 


The study was conducted at the South- 
eastern ‘Tidewater Experiment Station, 
south of Savannah, Georgia. The natural 
forests of the area, known as the flatwoods, 
are principally loblolly pine, intermingled 
with willow oak, longleaf pine, sweetgum, 
and reeds. 

The Bladen clay loam of the area is a 
low humic gley soil formed from thick 
beds of acid clay. It is noted for its slow 
internal drainage and plastic B horizon. 
The principal type of the series is fine 
sandy loam with slightly lighter color of 
the A horizon than the yellow-brown 
mottling typical of the underlying B. The 
A horizon is + to 10 inches thick and 
strongly acid. The texture of the B hori- 
zon varies from sandy clay to clay. When 
dry, it is extremely hard, shrinks and 
cracks, and breaks into angular blocks. 
The B horizon is about 6 feet deep in the 
study area. Bladen soils join areas of tidal 
marsh and occur on broad level flats with 
a few areas in ponded positions from 4 to 
30 feet above sea level. Drainage canals 
are considered essential for agriculture. 
The impermeability of Bladen soil is indi- 
cated by the very slight lowering of the 
water table at stand edges in contrast to 
the sharp decline one-half chain within a 
stand (Trousdell and Hoover 1955). 

Studies with replicated lines of draw- 
down wells, about 500 yards from the 
area of the experiment reported here, 
showed that a drainage ditch five feet deep 
had no effect on ground water profiles at 
distances greater than 12 feet from the 
ditch. Other studies indicate wide varia- 
tion in soils referred to as the “Bladen 
Series” along the South Atlantic Coast. 


7A new book “Physiology of Trees” by 
P. J. Kramer and T. T. Kozlowski (McGraw- 
Hill, 1960), published after submittal of this 
paper, gives an excellent literature review of 
this subject. 


in 
ins 
du 


e 


h 


O 


‘is 
of 


On Bladen and similar soils, differences 
in vegetative stocking are reported as affect- 
ing the drawdown rates of water tables 
during 2 rain-free summer period. The 
rates ranged from 0.14 to 0.09 feet per 
day, depending upon the proximity of the 
water table to root concentrations. Draw- 
down during the growing season, caused 
principally by evapotranspiration, appeared 
to be three to four times that of drainage 
alone.” 

Rainfall at the study area for 1957 was 
52 inches; and for 1958, 48 inches. 


Description of Main Experiment 


One-vear-old slash and loblolly pine seed- 
lings were planted at 2 x 2 foot spacing in 
mid-February 1957, in previously pre- 
pared 12 x 24 foot diked and ditched 
plots. The soil was undisturbed at time 
of planting. The whole plots were divided 
into 2 split plots, each randomly assigned 
25 trees of one species. Nursery-grown 
stock from seed originating in the coastal 
area Was used. 

A severe four-month drouth preceded 
planting. Therefore, to improve survival 
chances, all seedlings were watered within 
five days. “Twenty-four slash pines and 
13 loblolly pines, less than 3 percent of 
the 1350 trees planted, died and were 
replaced by early March. Between plant- 
ing and beginning of water treatments on 
May 31, average height of slash pine seed- 
lings had increased 3.5 inches (to 7.5 
inches) while loblolly had grown an aver- 
age of 5.9 inches (to 11.8 inches). All 
growth measurements were made to the 
nearest inch—from the ground to the 
terminal bud. The difference in_ initial 
size of the two species was of such magni- 
tude that subsequent results were expected 
to be strongly biased in favor of loblolly 
pine (Table 1). 

Initial watering following planting, re- 
planting, and favorable growing condi- 
tions for three and one-half months prior 


SGallup, L. E. Unpublished thesis. North 
Carolina State College, 1954. 





Figure 1. A drainage plot at the time treat- 


ments e@gan, 


to establishment of water treatments in- 
sured strong, vigorous plants for study. 
The design consisted of 3 blocks of nine 
randomized treatment plots each. 
Treatments (Fig. 1) were: 
(a) Water level maintained 8 inches 
above mean plot elevation (+8), 
(b) Water level maintained 4 inches 
above mean plot elevation (+4), 
(c) Water level maintained at mean 
plot elevation (0), 
(d) Water level maintained 4 inches 
below mean plot elevation (—4), 
(e) Water level maintained 8 inches 
below mean plot elevation (—8), 
(f) Water applied to a_ three-inch 
depth at three-week cycles. Out- 
lets were lowered daily 1% inch 
(drainage coefficient), plus evapo- 
ration (dr 4), 
(zg) Same as (f) with “%4-inch drain- 
age coefficient (dr 4 » 
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(h) Same as (f) with Y%-inch drain- 

age coefficient (dr %), 

(i) Uncontrolled check plot on un- 

disturbed soil (ck). 
Flooding treatments were begun in early 
June 1957 and maintained throughout 
the balance of the two-year experiment, 
except for (a), which was concluded at 
the end of the first year. Drainage treat- 
ments (f), (g), and (h) were suspended 
in} November 1957 and resumed on 
March 26, 1958, at the first indication of 
the breaking of dormancy. 

Water from an artesian well was sup- 
plied through plastic piping of one-inch 
mains and _ three-quarter-inch _ laterals. 
Water levels were controlled by inverted 
siphon outlets from dikes or ditches which 
encircled the plots. The water supply to 
each plot had independent manually- 
operated inlet valves. The constant water 
level treatments were maintained to offset 
evaporation, transpiration, and seepage. 
Valves were adjusted to provide slight 
water movement through overflow  out- 
lets which had sufficient capacity to re- 
move excess rainfall within 24 hours. 
Prior to lowering outlets in the drainage 
coefficient plots, evaporation was measured 
in a standard Weather Bureau pan for 
the preceding 24 hours and the amount 
added to the drainage coefficient. 

Check plots adjacent to the diked areas 
received no treatment other than water- 
ing at time of planting to insure survival. 
However, ponding occurred on these plots 
for periods of a week or more on several 
occasions after heavy rains, but the plots 
were never flooded to three-inch depths, 
as were the drainage coefficient plots at 
three week intervals. The check plots 
were then at an intermediate soil moisture 
level throughout most of the growing sea- 
son. The maximum variation in eleva- 
tions within plots was ().48 foot. Most 
plots did not vary more than 0.2 foot in 
elevation. 

Movement of water through Bladen 
soil is known to be very slow. Thus, it 
was not generally necessary to puddle core 
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trenches under levees, which were 18 
inches high with an 18-inch top width and 
a 1:1 side slope. Where fill material did 
not bond with top soil, causing seepage 
along lines of cleavage, cores were puddled 
through the top soil. Crayfish (Cambarus 
sp.) were a problem on some plots: they 
tunneled under levees, making holes over 
an inch in diameter through which water 
temporarily escaped. 

Survival counts were taken monthly and 
heights of seedlings measured at time of 
planting and periodically thereafter. 

Water temperatures on +4 and +8 
plots were recorded on five days each 
month and whenever atmospheric tem- 
peratures exceeded 90°F the first year, 
and about every three days during the 
second growing season. At the same time 
—about 3:30 pm—air temperatures were 
observed above the plots. 

Determination of dissolved oxygen in 
standing (semi-stagnant) water was made 
for samples withdrawn on three occasions 
from the +8 treatments. For compari- 
son, dissolved oxygen was determined for 
artesian water (source for treatments), 
water from a nearby running stream, stag- 
nant water ponded in the forest, and stag- 
nant water in an open field. The Solvay 
(1957) Method 22 (Winkler MnSO;) 
method was employed. 

‘The study was terminated after two 
years because it was evident that, on the 
plots with better growth, crowns of seed- 
lings would be interlocked before the end 
of the third growing season. 


Supplementary Study 


Because few seedlings survived in the 
three plots flooded to eight-inch depths 
from June 1957 to January 1958, further 
observations would serve no useful pur- 
pose. These plots were then drained on 
January 30, 1958, and a supplementary 
study installed in mid-February to deter- 
mine (a) whether slash and loblolly pine 
seedlings can survive if planted on water- 
logged Bladen soil immediately after sur- 
face drainage is provided, and (b) whether 
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398 Slash Ss ty 27 228 21.7 21.3 19 

Loblolly 8.7 et Be OS 24.3 22.3 20 
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431 Slash 7.7 18.0 7 21.8 21.7 21.3 19 
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459 Slash 7.3 17.7 21.7 21.0 21.3 21.3 19 

Loblolly S22 2h.7 Fas Fes 23.7 21.7 20 
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seedlings, subjected to 8-inch flooding for 
periods of 2, 4, and 8 weeks will survive 
after drainage. As only these three plots 
were available, range in duration of flood- 
ing was sought, rather than replication. 
One hundred seedlings of each species 
were planted in each plot at 1 x 1 foot 
spacing. During the two-week interim be- 
tween draining and planting, a_ crust 
formed on the soil which would tempo- 
rarily support a man; but the trampling 
during hand-planting left the soil mushy. 
Machine planting would have been im- 
possible under these conditions. As in 
1957, seedlings were permitted to grow 
until June before being subjected to flood- 
ing. Loblolly pine seedlings were from 
seed from a single tree, slash pine seed- 
lings were common stock. 


Results and Discussion 


Survival. Few seedlings died during the 
first 41 days, from May 31 to July 11, 
even on plots flooded to a depth of 8 inches 
above mean plot elevations. Nevertheless, 
highly significant differences were noted 
between species (Table 1). Appreciable 
mortality was first observed for slash pine 
on the +8 plots 68 days after treatment 
began. For loblolly pine, it was observed 
at 133 days on both +4 and +8 plots. 
Of the slash pine seedlings, 64 percent sur- 
vived at the end of the experiment. 
Seedling heights seemed important to 
survival. For the first year, loblolly pine 
seedlings survived better than slash pine 
under flooded conditions, possibly because 
the loblolly pines were taller when treat- 
ments began. By rearranging height and 
survival data, we found a relationship be- 
tween height of terminal bud above the 
water surface and the time required for 
flooding to kill pine seedlings during the 
first year (206 days), regardless of species 
(Fig. 2). Much of the difference in sur- 
vival occurred in the +8 treatments, 
where slash pine heights averaged 0.9 inch 
less than the water level. The loblolly 
seedlings for this treatment averaged 4.0 
inches above the water. Slash pine so 
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treated (+8) had survival rates signifi- 
cantly lower (1 percent level) than for 
all other treatments. Loblolly pine seed- 
lings averaging 6.6 inches above water 
level also had lower survival rates than 
those which averaged 11.0 inches above 
that point. Between other seedling heights, 
as related to water level, no significant 
differences either within or between species 
were apparent. The interaction between 
species and treatments, although statistically 
significant in August (68 days), September 
(101 days), and December (206 days) of 
the first year is, therefore, not of real 
importance, It is possible that the same 
reasoning would hold for the significant 
interaction found at the first measurement 
of the second year (245 days). 

Following the first 41 days of treat- 
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with respect to time, Treatments are shown 
in parentheses. The number of seedlings 
at the time flooding began was 75 in each 
case. LS.D.—5 percent= 8.6; 1 percent 
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ment, trees died at an increasing rate on 
plots flooded to depths of 4 and 8 inches 
and by 206 days practically all seedlings 
on the +8 plots were dead. Only 66 per- 
cent survived on the +4 plots. On the 
=() plots, survival was 87 percent, and 
all other treatments were above this ac- 
ceptable minimum. In the 27 days follow- 
ing the first +1, from July 11 to August 
7, the differences between treatments taken 
together became highly significant and re- 
mained so, At no time during the first 
crowing season Was there a significant 
difference in survival between the check 
plots, the three drainage coefficient treat- 
ments, or the three constant water level 
drainage treatments (+0, —4, —8). 

The amcunt of drainage appeared un- 
important to survival at the end of the ex- 
periment, as long as seedlings were not 
continuously flooded. In December 1958, 
550 days after establishing water control 
treatments, there were still no significant 
differences between treatments or species 
except for the +4 treatment. For slash 
pine, survival was definitely inferior on 
the +4 to the 4, 5, dr ! By and dr ! } 
plots. But no significance was shown be- 
tween these two drainage treatments and 
the dr ™% treatments. For loblolly pine, 
the +4 plots exhibited poorer survival than 
did all other treatments. (It will be recalled 
that the +8 plots were previously aban- 
doned.) Slash pine survival was_ better 
than loblolly pine (1 percent level) on the 
+4 treatments and inversely so on the 
+() plots, Yet, for all treatments com- 
bined, specits differences were not appar- 
ent. 

In the supplementary study average sur- 
vival 4 months after planting, when treat- 
ment began on June 6, was 69.3 percent 
for slash pines and 96.0 percent for lob- 
lolly pines. This is similar to the main 
study survival observation—that the slash 
pine stock employed here apparently re- 
quired a more favorable planting environ- 
ment than did the loblolly pine for equal 
chances of survival. Loblolly pine seed- 
lings temporarily withstood inundation 
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better than slash pine. On water-logged 
soil flooded to depths of 8 inches—approxi- 
mately the height of terminal buds—sur- 
vival was better for loblolly pine only as 
long as submergence continued (‘Table 3). 
This held especially for the 4- and 8-week 
inundation periods. The initial advantage 
of loblolly pine was then offset by its higher 
mortality after drainage began. In all 
three plots, survival of slash pine 6 weeks 
after the plots were drained slightly ex- 
ceeded that of loblolly pine, and this differ- 
ence persisted. Eighteen weeks after the 
2-week flooding period was over, slash 
pine seedling survival was 76 percent and 
loblolly pine 60 percent. This indicates 
that trees with an 8-inch inundation for 
2 weeks in June may have such a low 
survival as to require subsequent replanting 
(Fig. 3). It is possible, however, that 
flooding earlier in the year, when water 
temperatures are lower, would have less 


severe effects. 
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Growth, Height growth of pines in the 
flatwoods continues throughout the season 
in which pines are normally considered 
dormant. While it was first believed that 
early October approximates the inception 
of winter dormancy in the southeastern 
Tidewater area, trees in several plots 
measured on November 13, 1957 had 
grown twice as much in the 33 days since 
October 11 as in the previous 133 days. 
Because these data contradicted anticipa- 
tions, it was decided to make growth 
measurements monthly. As survival was 
low for plots flooded to a depth of 8 inches, 
growth data for these treatments are 
omitted to permit a more accurate analysis 
of variance between other treatments. 
On October 11, 133 days after estab- 
lishment of water treatments, the first 
height measurements were made and net 
growth of surviving trees determined. At 
that time slash pine height growth ranged 
from (0.3 inch on the +0 plots to 3.0 
inches on the —-8 plots. Loblolly pine 
averaged 0.7 inch on the +0 treatment 
and 4.0 inches on the check plots. Treat- 





Figure 4, Slash pine seedlings at conclusion 
of the experiment. From left to right: 
+0, +4, —8, and check treatments. 
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ment responses were significantly different 
(1 percent level), as was the interaction 
between species. Species x treatment inter- 
action, however, was significant at the 5 
percent level, while variation among repli- 
cates of the same treatment was not. 

The —8 treatment was best for slash 
and loblolly pine height growth combined. 
Highly significant differences were found 
between the +4 and —4, +4 and —§8, 
+4 and ck, +0 and —4, +0 and —8, 
+() and ck, —4 and dr “4, —8 and dr 
14, and ck and dr ™% treatments. For 
each succeeding date of measurement, and 
for each species, these statistical variances 
are readily determined from the figures on 
least significant difference (LSD) in 
Table 2. 

Treatment differences remained highly 
significant throughout the study. The sta- 
tistical difference between species dis- 
appeared after the first measurement; but 
early in the second year, when replication 
variance appeared, species difference be- 
came highly significant and remained so 
until the conclusion of the experiment. 
Near the end of the study replication vari- 
ance was found highly significant. No con- 
sistency was noted for the species x treat- 
ment interaction. 

At the end of the study, slash pine 
height growth ranged from 4.4 inches on 
the +0 plot to 19.7 inches where drain- 
age was to a depth of 8 inches (Fig. 4). 
Loblolly pines behaved similarly, being 
smallest (6.2 inches) on the £0 plot, and 
21.8 inches on the —8 plot. The trend 
in growth for slash pine from least to 
best was: +0, dr 4, +4, dr “4, dr kK, 
ck, —4, and —8. For loblolly, it was +0, 


/ 


dr 4, dr “%, dr YW, +4, ck, —4, and 
—§8. Extreme plot conditions at the end 
of the tests are shown in Figure 5. 
Throughout the first year loblolly pine 
seedling growth was best on the check 
plots. Not until over 15 months (459 
days) had elapsed was the —8 treatment 
found to be better. Perhaps the second 
growing season was more typical for judg- 
ing the effect of poor natural drainage 
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TABLE 2. 
constant water level treatments. 
ment began on May 31, 1957.) 








Average growth, in inches, of surviving seedlings with drainage and 
(Number of days refers to the time since treat- 


No Constant water levels Drainage coefficients heck LSD! Anal. of Varia 1 
lays Species +4 —0 4 8 le A 5% 1% Repl Trea Sp Sp. x 1 
133 Slash 0.6 0.3 2.0 3.0 1.4 06 0.6 1.6 
Loblolly 1.0 0.7 26 2.0 17 1.6 0.9 4.0 1.4 1.9 
§S+L 0.8 5 2.3 2.6 1.6 1.1 0.7 2.8 0.9 1.3 
ns ** * * 
166 Slash 2.0 1.7 3.7 4.1 2.7 1.6 2.0 2.4 
Loblolly Pe 1 4.1 3.1 ) 2.6 1.9 S.1 1.9 2.6 
S+ L 1.8 1.8 3.9 3.6 2.8 2.3 1.9 3.8 1.3 1.9 
s ** a 1s 
é Slash 2.9 2:3 4.2 5.0 3.7 2.5 27 3.3 
Loblolly 2.7 2.2 4.4 3.7 3.2 3.0 2.3 5.7 aes 2.1 
S+L 2.8 2.2 4.3 4.4 3.4 a9 2.5 4.5 0.9 1.2 
s ** “ 
245 ash 4.1 2.8 4.5 5.4 4.0 3.4 3.4 78 
Loblolly 3.0 2.7 4.7 2 3.4 3.3 2.6 6.1 1.4 1.9 
Ss I 3.6 2.7 4.6 4.8 oF 3.3 3.0 4.9 1.0 1.3 
s ** s * 
277 3.2 2.2 4.5 5.5 7 3.0 2.8 3.8 
1 5.5 2.6 4.7 4.2 3.4 3.3 2.6 6.1 ! o.3 
Ss I 4.3 2.7 4.6 4.8 3.6 3.1 2.7 5.0 1.2 1.7 
ns + s s 
304 Slash 3.2 2.8 4. 5.6 7 2.8 3.0 4.1 P 
Loblolly 5.3 2.7 4.8 3 3.4 3.3 2.9 6.2 1.5 2.1 
Ss I 4.2 a8 4.9 5.0 5 3.1 2.9 5.2 0.8 1.1 
s *+ s * 
335 Slash 3.5 3.1 8.2 9.7 4.7 3.2 i 6.9 - 
Lobloll 6.2 3.6 10.4 9 4.9 4.5 4.3 12.2 2.1 2.9 
S I 4.9 3.4 9.3 2.8 4.8 3.8 4.0 v5 1.7 2.3 
* * ° 
369 Slash 3.6 3.3 ) 11.5 §.2 3.4 4.1 7 P , 
Loblolly 6.3 4.2 13 13.3 5.6 3 §.2 14.3 «.) - 
S I 4.9 3.8 11.5 12.4 5.4 4.4 4.7 11.1 2.0 2.7 
* * “* * 
398 Slash 3.8 3.6 1.1 13.7 5.5 3.7 4.7 9.0 ie , 
Loblolly 6 5.1 5.8 16.1 ; §.7 6.2 16.5 - ’ 
S I $.2 4.3 3.5 14.9 3.9 4.7 5 12.8 2.6 3.6 
* * ° 
431 Slash 4.0 3.9 12.9 16.1 5.7 3.8 4.9 10.6 aa a 
Loblolly 7.1 5 18.2 18.6 6.8 6.0 6.7 18.8 .< 3.0 
s+ 1 5.6 4.7 5.5 17.3 2 4.9 5.8 14.7 2.8 3.8 
* * * . 
459 Slash 4.1 3.9 3.4 17.2 5.8 3.9 4.9 11.5 : = 
Loblolly 7.1 $.7 19.5 19.9 6.9 6.1 6.7 19.6 4.2 5.8 
Ss .. J 5.6 4.8 16.4 18.6 3 5.0 5.8 15.6 2.9 4.0 
* * . 
489 Slash 4.4 4.0 4.5 18.1 5.9 4.1 5.1 12.1 
Loblolly 8.9 5.9 0.2 21.0 7.0 6.5 5.8 19.9 3.6 4.9 
Ss L 6.7 5.0 7.4 19.5 6.5 §.3 5.9 16.0 2.3 31 
i . on” * on 
518 Slash 4.8 4.2 15.2 18.7 6.3 4.3 5.6 12.6 y : 
Lobloll 9.2 6.0 20.7 21.5 7.4 6.8 7.8 201 3.4 4.6 
s I 7 5.1 8.0 20.1 5.5 6.4 16.3 3.4 4.7 
* ** ne 
$50 Slash 4.8 4.4 15.6 19.7 6.4 4.6 6.3 13.0 s 
Lol | 1.6 6.2 24.2 21.8 7 6.9 7.3 20.4 3.2 4.4 
S I 7.2 5.3 18.4 20.8 6.9 $7 6.8 16.7 24 < 
ILSD Least Significant Difference at the respective levels c ing statistical iriance betwee iny two figures in the 
line, « ines, f which the LSD’s a gi O asterisk indicates 5 ent level a two asterisks indicate | percent level. 
since critical drouth conditions did not factory growth despite excess water the 


need to be overcome, as was the case the 
first year. If the first year had been nor- 
mal, it is doubtful that check plot growth 
would have been superior to the +() treat- 
ments. Firm seedling establishment dur- 
ing the first year, moreover, enabled satis- 


second year. 

The pretreatment first-year drouth left 
soil moisture depleted to the degree that 
abundant rainfall could be absorbed and 
retained by the soil without a serious flood 
condition occurring on the check plots. 
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Figure 5, An 8-inch plot, showing 


plot, loblolly pines at the end of the second growing season 
(left), and a +4-inch plot sith hydrophytic vegetation and stunted seedlings (right). 

Also, adequate rainfall well distributed drainage is needed other than the elimina- 
during the growing season favorably in- tion of permanent ponding. However, 
fluenced survival and growth. If the par- average growth on plots with water at 
ticular season sampled had been normal, —§ inches and —4 inches was almost as 
this might imply that for loblolly pine no much as for the check early in the second 
TABLE 3. Survival and growth of seedlings with eight-inch water level 


treatments in the supplementary study on water-logged soil. Growth calcula- 
tions (inches) are cumulative with time for surviving seedlings. Survival per- 
centages are based upon trees living at the begining of treatment. 


2-week 4-week 8-week 
treatment treatment treatment 
Slash Loblolly Slash Loblolly Slash Loblolly 
Sar. Gr. Ser. ‘Gr. Ser. Gr. Sur. Gr. Sur. Gr. Sur. Gr. 
Per- In- Per- In- Per- In- Per- In- Per- In- Per- In- 
cent ches cent ches cent ches cent ches cent ches cent che 
End of treatment 95 0.4 100 0.4 51 1.2 96 0.9 39 29 §3 £9 
End of treatment 
plus 6 weeks it 2. ae BS st. 22 2S Taz 29 3.4 24 2.3 
End of treatment 
plus 12 weeks 76 2.6 61 1.6 St 23 24 12 29 «(63S 22 23 
End of treatment 
plus 18 weeks 7 


6 27 GO 4.7 at. Gee. ele 7 ao 30 Zt. 23 
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growing season (335 days), and_ these 
three treatments were significantly better 
(1 percent level) than all others. 

Slash pine seedlings, in contrast, had 
their greatest growth during the entire 
experiment on —§ plots, although this was 
not significantly greater than for the —4 
treatment. It was, however, significantly 
more (5 percent level) than on the check 
plots. Slash pine growth on the check plot 
was at an intermediate level among all 
treatments. As stated previously, the vari- 
ation between species was probably due to 
the initial difference in the size and ap- 
parent vigor of the stock, the loblolly pine 
being much superior. This, of course, does 
not account for the disappearance of sig- 
nificance during the period from 166 to 
304 davs after treatments began. 

Growth was slightly better for both 
species on +4 plots than on the +0. 
While slash pine growth on —8 plots was 
slightly superior to loblolly pine the first 
vear, the opposite occurred the second year. 
The better growth of loblolly pine than 
slash pine on the drained plots may indi- 
cate that the former is benefited most by 
drainage in the seedling stage or, again, 
the larger stock may have been responsible. 

Another trend for both species under 
the coefficient treatments was towards 
higher average growth with the slowest 
drainage rate of “%-inch per day, though 
the differences were not significant at any 
time during the experiment. A coefficient 
of '-inch is lower than that now used by 
drainage engineers. 

Different water treatments appear to 
have had little effect on the duration of 
dormancy. A small amount of growth 
occurred in January, virtually none in 
February, and some in March. April was 
the period of spring flush, and through 
July of the second year, rapid growth 
continued on check plots and for —4 and 
—§ treatments. 

In the supplementary study, average 
growth of slash pine was greater than for 
loblolly pine at every period of measure- 
ment after plots were drained (Table 3). 


Both species averaged 7.2 inches in height 
at time of treatment on June 6. This was 
approximately 0.8 inch less than the depth 
of flooding. During each period of sub- 
mergence, growth of surviving slash pine 
seedlings was equal to, or greater than, that 
of loblolly pine ( Fig. $0 


Temperatures. High survival rate on 
one +8 plot the first year may have 
been due to its lower water temperatures. 
This plot also had abnormally favorable 
height growth, judging from the growth 
of seedlings in its replicates. Therefore, 
detailed temperature measurements were 
made about every three days during the 
second growing season to confirm the re- 
sults from limited data collected the first 
year. Should further testing substantiate 
that pines are so sensitive to water tem- 
peratures, a practical aspect would be the 
greater need for elimination of ponding 
in open fields than where partial shade 
occurs, 

Water temperatures were usually about 
to 5°F warmer than the air, but on 
occasion differences of 13°F were free- 


70 


corded. Prior to October, only for the 
single plot mentioned above was the water 
temperature cooler than the air, the differ- 
ence averaging 4°F. During October and 
November, ambient temperature generally 
exceeded water temperature. Tempera- 
ture in the air above flooded plots varied 
as much as 7-F while water temperatures 
differed as much as 11°F (Fig. 6). No 
reason is surmised for these differences 
during the first year of treatment other 
than the effect of some shading from a 
woodland on the west and a shrub hedge 
to the east of the plots. ‘The temperature 
variation between plots continued through 
1958, and again water temperatures tend- 
ed to be higher than air temperatures. 


Oxygen. Dissolved oxygen in the water 
of the +8 treatments, as well as in the 
comparison samples, was extremely in- 
consistent. In the flooded plots, it ranged 
from 2.8 to 8.3 mg per liter. Stagnant 
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Figure 6, Average air and water tempera- 
tures of plots flooded to depths of four and 
eight inches. 


ponds with timber contained about 5 mg, 
ponds in the open 7 mg, running streams 
7 mg, and the artesian water source ranged 
from 3.7 to 8.3 mg at the outlet. While 
poor aeration may have limited growth in 
the flooded plots, the amount of free oxy- 
gen present generally exceeded the re- 
quirements for growing plants hydro- 
ponically, 


Summary and Conclusions 


Slash and loblolly pine seedlings planted 
in a plastic clay loam soil were continuously 
flooded to + 0, 4, and 8-inch depths; con- 
tinuously drained to 4 and 8 inches below 
ground level; and flooded to a depth of 
three inches at 3-week intervals with sub- 
sequent drainage at 4g, 4%, and %-inch 
per day coefficients. 

A supplemental study tested survival 
ability of slash and loblolly pine seedlings 
planted on waterlogged soil immediately 
after drainage, but subjected to subsequent 
2-, 4-, and 8-week periods of inundation 
to 8 inches. 
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Results and suggested practical conclu- 
sions are: 

1. Mortality was inversely related to 
height of seedlings. Moreover, it appears 
that terminal buds must be above water to 
survive more than a few weeks of inunda- 
tion. 

2. Slash pine seedlings in the plots be- 
gan to die 68 days after initiation of flood- 
ing to 8-inch depths, and loblolly pine 
within 133 days. Flooding to 8-inch 
depths resulted in severe mortality, while 
flooding to 4-inch depths caused significant 
losses the first year. Slash pine inundated 
for periods longer than two months and 
loblolly pine for four months will probably 
require replanting. 

3. One-year-old seedlings of slash and 
loblolly pine averaging 7.2 inches tall en- 
dured a maximum of 2 weeks under 8 
inches of water without excessive mortality 
in the supplemental study. To permit seed- 
ling survival, previously water-logged sites, 
now drained, will require surveillance to 
guard against prolonged subsequent in- 
undation to depths of 8 inches and over. 

4. First year height growth of loblolly 
pine on plots continuously drained to + and 
8 inches below ground level considerably 
exceeded growth where water constants 
were maintained at ground level and at 
4 and 8 inches above that point. Best 
growth of slash pine was on plots with 
drainage to a depth of 8 inches, although 
it was not significantly different than for 
a 4-inch depth. By the end of the second 
year, greatest growth for both species was 
on the plots drained to 8 inches below 
ground level. Maintaining ground water 
at a depth of 4 inches, however, appears 
satisfactory for young seedlings. 

5. Among drainage coefficient plots, no 
real differences occurred between rates of 
water removal of 44, 44, and 4 inch per 
day. QOne-eighth inch coefficients appear 
satisfactory for seedling establishment. 

6. Growth of slash pine exceeded that 
of loblolly pine following drainage of all 
plots inundated from 2 to 8 weeks, indi- 
cating a greater ability of the former to 


M 


withstand unfavorable water conditions. 
7. Trees in several plots grew twice 
as much from October 11 to November 


13 as in the 133 days prior to October 11. 
Literature Cited 


AHLGREN, C. E., and H. L. Hansen. 
1957. Some effects of temporary flood- 
ing on coniferous trees. J. For. 55: 
647-650. 

BERGMAN, H. F. 1920. The relation of 
aeration to the growth and activity of 
roots and its influence on the ecesis of 
plants in swamps. Ann. Bot. 34: 13-33. 

Dean, B. Evizaperu. 1933. Effect of 
soil type and aeration upon root systems 
of certain aquatic plants. Plant Physiol. 
8: 203-222. 

DeMarEE, Devziz. 1932. Submerging 
experiments with Taxodium. Ecol. 13: 
258-262. 

GreEN, W. E. 1947. Effect of water 
impoundment on tree mortality and 
growth. J. For, 45: 118-120. 

Hunt, F. M. 1951. Effect of flooded 
soil on growth of pine seedlings. Plant 
Physiol. 26: 363-368. 

Jackson, W. T. 1956. The relative 
importance of factors causing injury to 
shoots of flooded tomato plants. Amer. 
& Bot. 43: 637-639. 

Koztowskl, T. T. 1958. Water rela- 
tions and growth of trees. J. For. 56: 
498-5()2. 

KRAMER, P. 5. 1933. The intake of 
water through dead root systems and 
its relation to the problem of absorption 
by transpiring plants. Amer. J. Bot. 
20: 481-492. 

1949. Plant and_ Soil 
Water Relationships. McGraw-Hill 
Book Co., N. Y. pp. 150, 220 et. fol. 
—————. 1951. Causes of injury to 
plants resulting from flooding of the 


soil. Plant Physiol. 26: 722-736. 

Kuster, A. 1948. Der Einfluss der 
Juraseen-Hochwasser auf die Strand- 
waldungen. Schweiz. Z. Forstw. 99: 
428-438, 

Leyton, L., and L. Z. Rousseau. 1958. 
Root growth of tree seedlings in rela- 
tion to aeration. In The Physiology of 
Forest Trees. Ronald Press, N. Y¥. 

PARKER, J. 1950. The effects of flood- 
ing on the transpiration and survival of 
some southeastern forest tree species. 
Plant Physiol. 25: 453-460. 

Pruirt, A, A. 1947. Study of effects of 
soils, water table, and drainage on the 
height growth of slash and loblolly pine 
plantations on the Hoffmann Forest. 
J. For. 45: 836. 

SATTERLUND, D. R., and §. A. GRAHAM. 
1957. Effect of Drainage on Tree 
Growth in Stagnant Bogs. Mich. For. 
No. 19. U. Mich., Ann Arbor. 

ScHLaAuptT, E. A. 1955. Drainage in 
forestry management in the South. In 
Water. Yrbk. of Agric., U. S. Dept. 
Agric. 

Scorr, A. D., and D. D. Evans. 1955. 
Dissolved oxygen in saturated soil. Soil 
Sci. Soc. Amer. Proc. 19: 7-12. 

Sotvay Co. 1957. Water Analysis. 
Technical and Engineering 
Bull. 11. 

Srong, E. L., R. R. Morrow, and D. S. 
We cu. 1954. A malady of red pine 
on poorly drained sites. J. For. 52: 
104-113. 

TrouspELi, K. B., and M. D. Hoover. 
1955. A change in groundwater level 
after clearcutting of loblolly pine in 
the Coastal Plain. i For. 53: 493-498. 

Yeacer, L. E. 1949. Effect of Perma- 
nent Flooding in a River Bottom Tim- 
ber Area. Ill. Nat. Hist. Surv. Bull. 
25, Art. 2. 


Se rvice 


volume 7, number 1, 1961 15 








Seasonal Variation in the Moisture 
Content of the Green Leaves 


of Mountain Laurel 


THE MOISTURE CONTENT of forest fuels 
is a major factor influencing the flammabil- 
ity of the forest fuel complex. Changes in 
moisture content exert a powerful and often 
controlling influence on fire behavior, Since 
most fires burn primarily in the dead fuels, 
it is understandable that fire control per- 
sonnel have been concerned almost solely 
with moisture content of the dead fuels. 
But live vegetation is also a factor in deter- 
mining how a fire behaves. As a fire be- 
comes more intense, more and more of the 
green understory vegetation becomes in- 
volved. Thus the moisture content of un- 
derstory vegetation may exert a controlling 
influence on fire intensity and may deter- 
mine whether or not a fire will pass from 
a surface fire to a crown fire. 

Little is known of how the moisture con- 
tent of green leaves varies with time 
through a fire season or a year. Unan- 
swered or partially answered questions in- 
clude these: How does the moisture content 
of new growth compare with the older 
leaves on evergreen species? How does 
dryness of site affect the moisture content 
of leaves? Are green leaves drier in dry 
years than in wet years? 

The present study was established to 
attempt to answer these questions for moun- 
tain laurel (Kalmia latifolia L.), an ever- 
green shrub common in the northeastern 
states. 

Mountain laurel is an understory shrub 
growing in dense thickets in the hardwood 
forests of the Northeast. As a forest fuel 
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type, it is classed as extreme in the rate of 
spread classification developed by Jemison 


and Keetch (1942). 


Previous Work 

Studies of seasonal variation in moisture 
content of evergreen shrubs have been 
carried out primarily in western forests and 
in Europe. New leaves of brush species of 
Arctostaphylos, Ceanothus, Adenostoma, 
and Quercus, growing in California and 
Idaho, contained moisture contents usually 
above 200 percent of their dry weight at 
the beginning of the growing season ( Fire- 
stop 1955, Richard 1940, Buck,’ ? and 
Fons*). But moisture contents decreased 


1Buck, he eas Progress report on Variations 
in the moisture content of green brush foliage 
on the Shasta Experimental Forest during 
1937. Calif. For. and Range Expt. Sta. unpub. 
mss. 1938. 

“Buck, C. C. Progress report on variations 
in the moisture content of green brush foliage 
on the Shasta Experimental Forest during 
1938. Calif. For. and Range Expt. Sta. unpub. 
mss. 1939. 

’Fons, W. L. Progress report on seasonal 
Variations in moisture content of chaparral on 
the San Dimas Experimental Forest during 
1942. Calif. For. and Range Expt. Sta. unpub. 
mss. 1943. 


The author is cn the faculty of the Yale 
Forestry School, Yale University, New Haven, 
Conn. He wishes to thank George M. Furnival 
for advice on statistical techniques. Manuscript 
received Jan. 11, 1960. 
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rapidly through the summer, reaching min- 
imum values of 65 to 120 percent in the 
fall, usually in September. Thereafter, leaf 
moisture increased slowly to a level of about 
100 percent in the winter and spring. 

These same studies showed that the old 
leaves of the same species reached a mini- 
mum moisture content in the late winter or 
spring, at values between 75 percent and 
100 percent. As the new flush of growth 
occurred, the old leaves also increased their 
moisture content, reaching a maximum of 
about 120 percent in the summer, then 
decreasing to 100 percent or less in the fall. 
There often followed a winter increase, 
with a second minimum in the spring just 
before the new growth appeared, 

Measurements on Rhododendron ferru- 
gineum, Calluna vulgaris, and Arctostaphy- 
los uva-ursi in Europe (Ulmer 1936) 
showed a similar cycle of leaf moisture 
content. 


Design and Procedures 


An even-aged oak forest of sprout origin 
with an abundance of mountain laurel in 
the understory was chosen as the experi- 
mental area. The oak stand, approximately 
57 years old at the time of the study, 
occupied the east-facing slope of West 
Rock Ridge, just above a small reservoir, 
Lake Wintergreen. The area, wholly 
within the Eli Whitney Forest, and just 
outside the northwest boundary of the city 
of New Haven, Connecticut, has been de- 
scribed in detail by Chapman.* , 

Three areas containing a heavy under- 
growth Of laurel and which appeared 
likely to have widely different soil moisture 
conditions were selected. One was on the 
flat at the bottom of the ridge, near the 
re servoir; this is referred to as the “‘bot- 
tom” site. Another, the “midslope”’ site, 
was about 400 feet distant on a gentle 
slope. The third, another 400 feet up the 

‘Chapman, G. L. The influence of moun- 
tain laurel undergrowth on environmental 
onditions and oak reproduction. Yale School 
of Forestry, Ph.D. dissertation. 157 pp. 1950. 


slope, on a rocky knoll with very shallow 
soil about 50 feet above the level of the 
reservoir, is referred to as the “‘knoll”’ site. 

Leaf samples were taken at irregular 
intervals from March 1957, through No- 
vember 1958, more frequently during 
spring and summer months than during 
fall and winter. Separate samples were 
taken from two clumps of bushes 25 to 100 
feet apart on each of the areas, and leaves 
were collected from various levels on the 
plants. Collections were made from dif- 
ferent plants each time, in order not to 
defoliate any one plant. New leaves were 
collected and analyzed separately from old 
leaves as long as it was possible to distin- 
guish them by their color and feel, and as 
long as analysis showed that their moisture 
contents differed measurably. 

Collections were made in early after- 
noon, as close to 2:00 pm as possible, in 
order to obtain leaves when they could be 
expected to have the lowest moisture con- 
tent of the day. Systematic bias was elimi- 
nated as much as possible by varying the 
sequence in which collections were made. 
Total time for collecton was about one- 
half hour. Collections were not made on 
the day of a rain, and care was taken to 
see that no leaves were taken that had 
obvious surface moisture. 

The samples, each of which consisted of 
approximately 40 grams of leaves, were 
analyzed for moisture content by the xylene 
distillation method described by Buck and 


Hughes (1939). 


In order to study possible effects of soil 
moisture on leaf moisture content, a run- 
ning soil moisture budget was maintained 
for the period of the study. Inasmuch as 
the main purpose of the soil moisture budget 
was to enable rough comparisons between 
the two years of the study, it was not con- 
sidered essential to have absolute soil mois- 
ture values. 

Accordingly, a climatological method 
was selected which maintained the budget 
in terms of potential evapotranspiration and 
rainfall. Although numerous such methods 
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have been developed, I used the one de- 
veloped by Hamon (Travelers Weather 
Research Center, 1957), which has given 
good results in Connecticut. 

The soil was considered to be at field 
capacity at the end of the winter, on 
April 1. For each day, estimated evapo- 
transpiration was subtracted and the rain- 
fall amount added to the budget. In Ha- 
mon’s method, potential evapotranspiration 
is considered to be a function of day length 
and the difference in water vapor density 
at daily maximum and minimum tempera- 
tures. Following Palmer and Havens 
(1958), I reduced computed potential 
evapotranspiration values by the same per- 
centage as the previous day’s soil moisture 
percent to gain an estimate of actual evapo- 
transpiration. It was considered that field 
capacity of the soil at all three locations was 
two inches of water. The assumption of 
constant field capacity had the effect of 
diminishing possible estimated differences 
between years on the dry sites as compared 
with the areas with deeper soil. This was 
not considered serious inasmuch as the two 
years showed marked differences in esti- 
mated soil moisture and the use of different 
field capacities for the three areas would 
have had the effect only of emphasizing 
these differences on the shallower soils, 

Rainfall amounts used in the soil mois- 
ture budget calculations were estimated 
from precipitation recorded at three near- 
by stations: Mt. Carmel (4% miles); New 
Haven Airport (734 miles); and Lake 
Dawson (34 mile). Despite the closeness 
of the Lake Dawson guage, it was not used 
alone because it was on the opposite side 
of the ridge from the experimental area. 
Weighted averages were calculated, using 
weighting factors inversely proportional to 
the distances of the several weather stations 
from the experimental area. 

Maximum and minimum temperatures 
used in computing the budget were those 
recorded at Mt. Carmel, the nearest tem- 
perature station in a similar topographic 
position. Rainfall and temperature data 
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were taken from Climatological data, New 
England (U. S. Weather Bureau, 1957 
and 1958). 


Results 


Individual moisture content values of the 
laurel leaves are contained in Table 1. The 
data are plotted in Figure 1, along with 
soil moisture values expressed as a percent 
of assumed field capacity. Each knoll site 
point is an average of two samples. These 
are labelled “dry site” in the figure. Other 
points are averages of four samples, as ex- 
plained in the next section, and are labelled 
“moist site” in the figure. 


Analysis 


In order to test the hypothesis that there 
were no differences between the three 
areas, an analysis of variance was performed 
on the old-growth leaves (‘Table 2). Dun- 
can’s (1955) test for the significance of 
differences in ordered means showed that 
the difference between means of the bot- 
tom and midslope sites was not significant, 
whereas the knoll mean was significantly 
different from both of the others. <Ac- 
cordingly, data from the two lower sites 
were grouped. In Figure 1, these data are 
labelled “‘moist site.” 

The analysis also disclosed that there 
was significant variation with time, a con- 
clusion that is evident from inspection of 
the graph. The non-significant interaction 
can be interpreted as indicating that the 
three sites behaved similarly throughout the 
year. 

A similar analysis performed on the new 
growth leaves (Table 3) yielded the same 
conclusions. In this case, the rapid decrease 
of moisture content as the leaves matured 
contributed to the large mean square for 
time. 

The year 1957 was one of extreme 
agricultural drought in Connecticut. Ha- 
mon’s analysis (cited previously) showed 
that, “. . . the drought of 1957 was more 
severe than any combination of the driest 
growing seasons during the past 30 years.” 
The following year, 1958, on the other 


TABLE 1. 


. : : 
Moisture content of laurel leaves in percent of dry weight. 




















Date Current year’s leaves Old-growth leave 
llected Knoll Midslope Bottom Knoll Midslope Bottom 
Percen 

3/29/57 104 132 126 117 29 114 
4/12/57 118 116 126 126 22 115 
4/19/57 19 114 138 119 36 118 
4/26/57 09 110 112 111 120 113 
§/ 3/57 106 97 112 97 04 112 
$/ 9/57 100 105 121 108 23 110 
$/24/57 299 287 351 347 288 301 23 115 133 76 24 121 
6/ 7/57 224 192 273 243 260 321 129 12 143 141 144 131 
6/18/57 160 166 203 194 218 202 118 132 146 152 $3150 
7/10/57 141 140 154 176 156 155 33-113 149 158 $4 144 
8/ 6/57 136 133 157 150 50 150 
8/13/57 08 114 143 142 52 149 
8/29/57 28 128 155 146 50 162 
9/24/57 41 32 157 155 61 162 
10/10/57 137 139 157 141 56 156 
11/14/57 130 134 140 134 30 138 
2/ 3/58 30 124 130 127 129 124 
3/ 7/58 26 128 134 133 130 138 
4/ 4/58 31 127 129 134 16 127 
4/18/58 121 12 127 126 25 125 
4/30/58 18 120 116 111 29 113 
5/12/58 170 211 239 227 229 222 100 107 96 103 100 110 
3/26/58 269 286 35 291 274 280 9S 97 98 99 106 100 
6/16/58 213 232 261 237 231 263 21 111 135 128 30 128 
7/ 3/58 203 351 204 218 203 230 56 154 170 163 40 174 
7/16/58 169 176 179 199 178 177 $3 160 160 144 17 51 
7/29/58 153 176 191 170 227 175 154 164 161 164 35 164 
8/11/58 143 135 171 154 67 166 
9/22/58 156 131 154 160 73 154 
10/31/58 134 129 140 128 123 144 


I1Two values shown for each category. 


hand, had a plentiful supply of precipitation : 
52 inches at New Haven, compared with 
a normal of 46 inches. This fortuitous cir- 
cumstance gave an excellent opportunity to 
assess the effect of two greatly different 
moisture regimes on leaf moisture content. 
It should be pointed out, however, that 
although the 1957 summer was very dry 
for southern New England, it cannot be 
compared with the typical long dry season 
experienced on the West Coast of the 
United States. 

In order to test the hypothesis that there 
was no correlation between soil moisture 
differences and leaf moisture differences 
between the two years, the following pro- 
cedure was used. Monthly averages of the 
laurel leaf moisture content were computed 
for the knoll site and for all sites grouped. 





Values for 1958 months were then sub- 
tracted from the corresponding 1957 
months to give a series of differences. A 
similar series of differences was calculated 
These differ- 
ences are compared in Table 4. 


from soil moisture values. 
From an 
inspection of the soil moisture values, it is 
evident that 1957 was indeed a drier year 
than 1958, at least during the months of 
March through July, confirming Hamon’s 
work cited earlier. 

The apparent lack of correlation be- 
tween soil and leaf moisture is confirmed 
by regression analyses of leaf moisture on 
soil moisture which yielded correlation co- 
efficients of 0.21 and 0.08, neither of which 
is significant at the 5 percent level. It may 
be argued that calendar months are not 
the appropriate time divisions inasmuch as 
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TABLE 2. 





Analysis of variance of old-growth laurel leaves. 


Source of Degrees of Sum of Mean Varian 
iriation freedom squares square 10 
Lime 29 52.156 1.798 23.57% 
Site 2 3,824 1,912 25.07** 
I S 58 5.746 99.1 1.30 
Error 90 6.864 76.3 
Total 179 68,590 

wes) le Be 
Site percent; midslope—134.7 ts boe 135.8 percent 
‘yy ) ? jcaenins , seein . ’ —_— arent Me — 2el2e 
TABLE 3. Analysis of variance of new-growth laurel leaves. 
Ss e of Degrees of Sum of Mean Varian 
ition freedom squares square ratio 
lime 9 134,395 14,933 23.49** 
Site 2 5,678 2,839 +.47* 
I S 18 18,807 1,045 1.64 
Error 30) 19,071 636 
Total 59 177,951 
*Significant at 5 percent level. 
**Significant at 1 percent level. 
Site means: Knoll—140.6 percent; midslope—155.4 percent; bottom—153.0 percent. 


plants are influenced in their growth by 
factors other than the calendar. However, 
growth stages in these two years did not 
differ by more than about two weeks. In 
contrast to this, the period of drought lasted 
five months. Furthermore, the 
drought was a severe one for the southern 


since 


New England region, any important leaf 
moisture differences would be expected to 
show up even when analyzed on a monthly 
basis. 

Surprisingly, and contrary to expectations, 
moisture content of laurel leaves was high 
in the driest months, May and June 1957, 
as compared with the same months in 
1958. This contrasted with soil moisture, 
which was low during May and June 1957 
as compared with the same months in 1958. 

A further attempt was made to corre- 
late leaf moisture content with potential 
evapotranspiration calculated for the day 
in which the sample was collected. These 
calculations, which are not shown, also 


failed to demonstrate a relationship between 
leaf moisture and drying power. 


TABLE 4. Comparison of leaf and 
soil moisture contents for a dry year 
and a wet year. 
1958 value). 


(1957 value minus 


Average Average laurel leaf 
soil moisture content 
moisture Dry site All sites 
Month (percent) (percent) ( percent) 
Mar. — 4 — 9 —12 
Apr. —12 —10 — 4 
May —24 8 8 
June —31 10 12 
July —16 —34 —16 
Aug. 5 —15 —Il4 
Sept. —11 — 8 — 4 
Oct. 8 6 15 
Correl. coefficients! 0.21 0.08 
Value for significance 
at 5 percent level 0.63 0.63 
1Coefficients of correlation between leaf moisture 


differences and soil moisture differences. 
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Although it is contrary to common sense 
and expectation, there appears to be no 
evidence in these data that there is a rela- 
tionship between laurel leaf moisture and 
the dry ness of the year. This is especially 
surprising in view of the fact that there 
apparently is a leaf-moisture difference be- 
tween plants growing on a rocky area and 
nearby plants growing on an area with 
deeper soils and presumably a greater mois- 
ture supply. It may be that plants on the 
rocky area have undergone some _ physio- 
logical conditioning or even a genetic 
change which is reflected in lower average 
leaf moisture contents as compared with 
nearby plants on better-watered areas, This 
internal conditioning may be _ powerful 
enough to mask year-to-year variations, at 
least with the moisture conditions prevail- 
ing in southern New England. 

It seems unlikely that there was no 
moisture stress on these plants inasmuch 
as agricultural plants throughout the area 
on much deeper soils did exhibit serious 
drought symptoms. However, and_ this 
may be important, the moisture stress on 
these understory laurel plants may be much 
less than the overstory vegetation because 
there is very much less net radiation avail- 
able for evapotranspiration under the shade 
of the overstory. There is no reason to be- 
lieve that the plants on the rocky area had 
any hidden moisture supply although it is 
possible that the plants on the lowest site 
were able to tap water from some depth. 
This seems unlikely, however, in view of 
the fact that all three sites responded 
similarly as indicated by the non-significant 
time-x-site interaction. It would certainly 
be unwise to make generalizations on the 
basis of this one example. However, the 
relationships (or lack of them) revealed in 
this study are sufficiently unorthodox to 
warrant further study. 

Conclusions 
The following specific conclusions may be 
drawn from the study: 

1. Green laurel leaves from plants 
growing on the rocky, shallow soil site were 
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consistently drier than those growing in 
deeper soils. The differential was about 
15 percentage points, greater in summer, 
less in winter and early spring. 

2. Old-growth laurel leaves varied in 
moisture content cyclically throughout the 
study period, reaching minimum values in 
May and maximum values in late summer. 
During the two years of measurements, 
moisture content of laurel leaves was about 
100 percent during the spring fire season 
and about 140 percent during the fall fire 
season, 

Ee No clear 
drought and leaf moisture content could 
be demonstrated. Visual inspection of the 


relationships between 


graphs of old-growth leaves seems to indi- 
cate slightly lower moisture contents in the 
dry year, but the differences do not stand 
up under statistical analysis. It should be 
remembered that 1957 was an extraordi- 
narily dry year for Connecticut (although 
much wetter than a typical California 
summer) whereas 1958 was well watered, 
with more than normal rainfall. 

4. New leaves of laurel started with a 
moisture content of about 300) percent, at 
least when they were large enough and 
numerous enough to collect. Their mois- 
ture content decreased steadily, approach- 
ing that of the previous year’s leaves 
asymptotically, reaching old-leaf values in 
about two month’s time. 
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Recently from East Germany has come a new of nitrogen, and resistance to brown scale (Eu- 
publication on fertilizing forests. This book, lecanium rni) of black locust; also the effect 
with an introduction by Professor W. Erteld of rabbit and h lamage on red oak as in- 
and a foreword by Professor W. Wittich, dis- fluenced by fertilization 
cusses the effect of potash and magnesium ad- B PP ee ee ' ‘ 5 
ee — I secause of the limited number of reports on 
litions on the development of Pinus sylvestris, aia’ i oF 
Sa J : potassium and magnesium fertilizer responses 
Robinia pseudoacacia, and Ouercu prealis. ea ‘ wi a 

f ~ in forestry, this book is a real contribution to 


Ihe three major sections of the book deal , ; : = 
. : ; ; ‘ > the literature. It substantiates previous work 
with (1) experiments with Scotch pine ir the : : Ae ie 
. oT 1; > . . that potassium and magnesium fertilization 
region of Templin, (2) experiments with i 





black locust, red oak and Scotch pine in the a : mo \ red in the amelioration of de- 
Schreen forest district, and (3) investigations I — er 

concerned with the phenomena and appear- Phe bok Comterne numerous tables, charts 
ances of nutrient deficiencies in the three tree and photographs. Of special note are the 33 
species. pages of ex ellent olor photogr iph plat S$ le- 
In addition to discussing growth of vegcta- picting nutrient deficiencies and growth re- 
tive responses of the investigated species erow- sponses oObtaine 1 as well as profile s of soils 
ing on degraded sandy soils, poor in nutrients supporting the nutrient-deficient trees. There 
and influenced by subsoil water, there is also are 138 literature citations. 

a discussion of generative development of black With its high standards of preparation, illus- 
locust. Mentioned are the effects of potash- tration, typography and paper, the book will 
magnesia fertilization on flowering, frost re- be of great value and interest to those con- 
sistance, defoliation, depth of root penetra- cerned with forest fertilization and tree nu- 
tion, abundance of root nodules and quantities trition, 
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Vegetative Indicators of Rate of 
Height Growth in Ponderosa Pine 


DurinG the course of a study aimed at 
delimiting the basic categories of forest 
vegetation of eastern Washington and 
northern Idaho, the writer noted what ap- 
peared to be consistent differences in the 
rate of height growth of ponderosa pine 
(Pinus ponderosa Laws.) as it occurs nat- 
As that 
study neared completion, the U. S. Forest 


urally in different habitat types.’ 


Service, through the helpful cooperation of 
Charles A. Wellner and Donald W. 
Lynch, made available height and age de- 
terminations of a series of P. ponderosa 
stands that had been studied in connection 
When these 


stands were revisited and classified by the 


with silvicultural research. 


writer according to habitat types, the exist- 
ence of a correlation with growth was con- 
firmed, but it was so weak that plant com- 
munities seemed to have little value for pre- 
dicting growth rate of the pine. 

In this preliminary study the basic height 
and age data for the different stands were 
extrapolated to a common stand age of 100 
years (the site index), using data from a 
table presented by Meyer (1938) which 
had been offered for this purpose. If such 
data should not accurately reflect the 


1A habitat type (syn.: homoece, equivalent 
environment) is all the area that supports, or 
is capable of supporting, the same kind of cli- 
max vegetation in the absence of disturbance 
grazing, etc. Each 
habitat type will be referred to by the name 


by fire, logging, heavy 


of the climax plant association that character- 


izes it, and the abbreviation h.t. will be used 


freely. 
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growth rate of a tree at different ages, then 
the extrapolated height at 100 years would 
be proportionately inaccurate. With this 
thought in mind the writer subsequently 
tried plotting directly the raw meaurements 
of height and age, keeping the data for dif- 
ferent habitat types on separate scatter- 
grams. When this was done rather sharp 
differences among certain of the habitat 
types were revealed, Thus encouraged, the 
writer began to study additional stands so 
as to show the extent of differences in po- 
tentialities among the habitats more ade- 
quately, and to investigate the reason why 
the site index method is inferior to the use 
of raw data. 

Methods 

Since the Forest Service data were not ob- 
tained with the idea of sampling equally all 
the seven habitat types in which P. pon- 
derosa may be found abundantly in eastern 
Washington and northern Idaho (Dau- 
benmire 1952), it was necessary to expand 
the initial nucleus of data so that all habi- 
tats were adequately represented, Stands 


The author is on the faculty of Washington 
State University, Pullman. 

In addition to the cooperation of the U. S. 
Forest Service and T. S. Russell mentioned in 
the text of this report, Jean B. Daubenmire 
shared in most of the field work. A part of the 
field expenses were met with funds provided 
for biological and medical research by Wash- 
ington State Initiative 171. Valuable sugges- 
tions for improving the manuscript were made 
by D. W. Lynch. Manuscript received Feb. 
27, 1960. 


abnormally open or extremely dense were 
avoided, and the height of at least five 
trees in or projecting slightly above the 
main canopy was measured with an Abney 
level. Age was assumed to equal the num- 
ber of wood layers counted on borings 
made at breast height, plus ten as an esti- 
mate of the time required for the seedling 
to have attained this height. It would have 
been better to have used a larger figure 
for age at breast height for trees growing 
slowly on poor sites, and a smaller figure 
for good sites, but in a study to distinguish 
poor from good sites this would hardly be 
an objective procedure without more evi- 
dence on rates of seedling growth than is 
now. available. 

The data were plotted as scattergrams, 
a separate one for each habitat type. To 
facilitate visual comparisons the outlines of 
one scattergram representing relatively 
rapid growth (Pinus ponderosa/Sym phori- 
carpos albus h.t.) and one representing rel- 
atively slow growth (Pimus ponderosa 
A gropyron spi atum h.t.) were then super- 
imposed on each of the other scattergrams 


(Fig. 1). 


Results and Discussion 


It was found that the habitat types could be 
divided into four unequal groups on the 
basis of rate of height growth of Pimus 
ponde rosa: 


I. BEST POTENTIALITY 


Abies grandis/Pachistima myrsinites 
At 


Il. GOOD POTENTIALITY 


Pseudotsuga menziest/ Physocarpus 
malvaceus h.t. 

Pinus ponde rosa/Physoce ar pus malva- 
ceus h.t. 

Pinus ponderosa Symphoricarpos albus 
h.t. 


Il. FAIR POTENTIALITY 


Pseudotsuga menziest/Calamagrostts 


rubescens h.t. 


IV. LOW POTENTIALITY 


Pinus ponderosa/Purshia tridentata 
h.t. 

Pinus ponderosa A gropyron spac atum 
h.t. 

As pointed out earlier, it seemed neces- 
sary for the sake of objectivity to count 
wood layers at breast height and then to 
add 10 to account for the years required 
for the seedling to have vrown to this 
height. Not only does the standard value 
of 10 mask a tremendous variation, but on 
unproductive habitats the value should be 
higher than 10, and vice versa. It is there- 
fore important to note that if seedling as 
well as older trees grow slowly on poor 
sites, the use of 10 for both good and poor 
sites reduces the apparent difference be- 
tween them. Lynch (1958) estimated the 
following values as appropriate for the 
breast-height age: 


Site Years at breast height 
I 3 

I] 5 

IIT 7 

IV 10) 

V 14 

VI 16 


If it may be assumed that his Site IV is 
roughly equivalent to the Pinus Sym phori- 
carpos h.t. and his Site VI to the Pimus/ 
A gropyron h.t., and the vertical positions 
of the scattergrams compared in Figure 1 
are adjusted accordingly, then about half 
the already small overlap is eliminated. 

The study has shown that the growth 
rate of the pine bears no relation to its suc- 
cessional status, for it grows well as a seral 
tree on burned or logged parts of the 
Pseudotsuga/Physocarpus h.t., and as a cli- 
max tree on Pinus/Symphoricarpos and 
Pinus/Physocarpus h.t. 

Although available moisture is a very 
potent factor governing the distribution 
and physiology (Daubenmire 1943, 1946) 
of trees in the Rocky Mountain region, site 
quality is not related to moisture alone The 
Pseudotsuga/Calmagrostis flora is less xero- 
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Figure 1. Each point on the scattergrams represents the average height and estimated average age 
of a group of at least five trees. The gross outlines of the pattern of points for the Pinus/Sym- 
phoricarpos avd Pinus/Agropyron habitat types have been superimposed on the data for the others 
to facilitate visual comparisons. Data for a few stands fall outside the ranges of the values for 


these C rdinates and abctssas and so have been omitted, 
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phytic than the Pinus/Symphoricarpos and 
Pinus/Physocarpus floras but the growth of 
the pine is inferior in the first. Since Pseu- 
dotsuga/Calamgrostis sites are characteristic 
of coarse glacial till, there iS clear indication 
that lack of fertility more than offsets the 
superior water relations of this habitat type. 

The Pseudotsuga/Calamagrostis habitat 
type as represented in eastern Washington 
and northern Idaho always has a_ fair 
amount of shrubs belonging to the Sym- 
phoricarpos and Physocarpus unions (see 
key to habitat types for taxonomic lists), 
rather than the apparently pure sward of 
Calamagrostis so characteristic of Pseudo- 
tsuga/Calamagrostis stands in the foothills 
of the Cascades, in Montana, or in central 
British Columbia (Tisdale and McLean 
1957). Thus the habitat type in eastern 
Washington and northern Idaho is repre- 
sented chiefly by intergrades between Pseu- 
dotsuga/Calamagrostis and Pseudotsuga 
Physcoarpus h.t. Considering all stands in 
which Pseudotsuga is the sole climax domi- 
nant, it was found empirically that if there 
is a fair amount of Calamagrostis present, 
and the stand is in the glaciated northern 
half of the study area, the growth rate of 
the pine is markedly reduced in comparison 
with sites where this grass is nearly if not 
completely absent. Therefore the hybrid- 
type stands containing at least a fair amount 
of Calamagrostis and located on glacial 
soils are to be kept separate from typical 
Pseudotsuga/Physocarpus stands, whether 
the latter are on glaciated or unglaciated 
soils, and from Pseudotsuga sites that have 
a fair amount of Calamagrostis but are on 
unglaciated soils. 

South of the glacial border where soils 
nearly always include an admixture of loess, 
Pinus/ A gropyron sites are rare, but north- 
ward they are abundant, alternating with 
Pinus/Symphoricarpos growing on what 
may be interpreted as patches of more 
favorable soils. In this northerly sector 
many forests have patches of the two under- 
growth types forming a mosaic of small 
units. It was discovered that even if only 


a moderate share of the mosaic is repre- 
sented by the Symphoricarpos type of un- 
dergrowth, the pine grows as well as if the 
shrub layer were continuous. 

The height of the group of shrubs con- 
sidered as belonging to the Sym phoricarpos 
union varies from about 0.5 to 1.5m, and 
from recorded estimates of shrub height in 
many stands an attempt was made to relate 
this variable to differences in rate of height 
growth of pine within the Pimus/Symphori- 
carpos h.t. No relationship was found. 

Many stands of Pinus ponderosa in the 
area of study are heavily infected with 
campylopodum, 
other areas completely lack this pathogen. 
Before the growth rate study was far along 


Arceuthobium whereas 


it was discovered that forest undergrowth 
is closely correlated with susceptibility of the 
pine to the disease caused by Arceuthobium, 
infection being common on Pinus/ A gro- 
prron and Pinus/Purshia habitat types, but 
not seen elsewhere. It is interesting to note 
that Shepherd (1959) found a_ parallel 
correlation between the susceptibiilty of 
Picea engelmanni to insect infestation and 
the character of the forest undergrowth. 


Validity of the Standard Site 
Index Curves 


As mentioned earlier, raw data show bet- 
ter correlations between growth rates of 
Pinus ponderosa and habitat types than are 
obtained after the site index curves have 
been used to extrapolate to a common age 
(Table 1). Such a situation casts doubt on 
the validity of the curves for P. ponderosa 
in eastern Washington and northern Idaho, 
and calls for testing them in other ways. 

The best opportunity for testing was 
provided at the edges of old burns where 
several decades ago fire had stopped part 
way through a pine stand permitting estab- 
lishment of a new generation of trees ad- 
jacent to the older residual stand on what 
appeared, from considerations of topog- 
raphy, soils and non-arborescent vegeta- 
tion, to be parts of one homogeneous habi- 
tat. If the site index curves are valid, the 
values for the contiguous young and the 
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old stands should be closely similar, but 
among the 11 opportunities to make such 
a comparison there were errors up to 28 
percent in prediction (Table 2). It is to 
be noted that the older stand always has 
the higher site index. The evidence thus 
indicates that the shape of the postulated 
growth curves is materially different from 
the actual growth rates of P. ponderosa in 
the region under study. Within the age 
range studied, growth is more rapid than 
postulated (See also Fig. 2A). Although 
these “‘fire-line”’ data obtained at the edges 
of old burns play havoc with the site index 
curves, they are all included in Figure 1 
where they fit harmoniously with other 
measurements made on the corresponding 
nt. 

Another manner of testing the validity 
of the site index curves is to count wood 
layers of felled trees at different heights 
along the axis, and construct curves to be 


TABLE. 1. Classification of stands of 
Pinus/Agropyron and Pinus/Sympho- 
ricarpos habitat types according to site 
index. Note that modal values almost 
overlap. 


Numbers of stands 
Site index Pinu Pinu 


Class Agropyron  Symphoricarpo 


40- 45 
46- 5( 
§i- 
56- 60 
61- 65 
66- 70 
71-75 
76- 80 
81- 85 
86- 90 
91- 95 1 
96-100 

101-105 

106-110 

111-115 

116-120 

121-125 1 


Totals 2 
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compared with those based on mathematical 
manipulation of height and age data for 
many stands. Only twice was there op- 
portunity to make this type of an analysis, 
but the results (Fig. 2A) of both tests are 
in complete agreement with the results of 
“fire-line” studies, namely that the pine 
grows much faster than the site index 
curves suggest. 

One might easily suspect the validity of 
any site index table that does not rest on 
an ecologic foundation. It has been firmly 
established that the species of the taxono- 
mist is usually not an ecologic unit, but 
rather a mosaic of ecotypes embracing as 
much heritable variation (in growth rate, 
for example) as one finds between species. 
The distinctive combination of climate and 
soil of each habitat type exerts a distinctive 
selective pressure on biotypes and hence 
tends to produce its own ecotype. The clear 
differences in growth data for the different 
habitat types may represent such genetic 
differentiation, or simply ontogenetic modi- 
fication, or a combination of the two. Tests 
under controlled conditions have demon- 
strated that contiguous habitats may sup- 
port genetically differentiated ecotypes 
(Daubenmire 1960, McMillan 1956, 
Bradshaw and Snaydon 1959, Habeck 
1958). The site index concept recognizes 
but one family of harmonic curves repre- 
senting dozens of distinctive habitat types 
ranging from South Dakota to the foothills 
of the Sierras and Cascades, and from the 
Canadian border to Colorado. Such an 
expression of growth ignores an important 
principle of ecology, and it is significant to 
note that others too have found it imprac- 
tical to use the same site index curves for 
contiguous habitats (Carmean 1956, 
Lorenz and Spaeth 1947). 

Spurr (1956) has sounded the extreme- 
ly important warning that nearly all work 
done to establish relations between soil pro- 
file characters and tree growth has involved 
passing the raw data through some stand- 
ard site index manipulation, All such re- 
search stands in jeopardy unless the site 
index curves accurately reflect the growth 


ABLE 2. Comparison of site indexes for 11 pairs of stands of different ages 
growimg on contiguous parts of what appear to be identical habitats. 


Relative \ge 

\duc- H na 

Pinu 104 

4 ere pI ro? 56 

P Ni 134 

1 gropyron 56 

I NM 149 

Pu 37 

Pinu 208 

Pur 61 

Pseudotsug 117 

c 1 g? si) 

Pseudot ug 64 

ld Calamagrost 37 

Fair Pp aos 9 
S@ud iga 

Calamagrost 44 

Pseudot t2a 93 

Calamagrost 41 

Pinu 91 

Sym phoricar pos 56 

' Pinus 97 

Good ? 2¢ 

Physocar pus 39 

Pinus 72 

Physocarpus 46 


rate of the tree on each of the soil categories 
studied! Perhaps this simple but important 
procedural error accounts for some if not 
all alleged failures of the vegetative indicator 
concept in the U.S.A. despite its general ac- 
ceptance in northwestern Europe and south- 
eastern Ganada. As pointed out above 
(Table 1), vegetative indicators do indeed 
appear to have low value in eastern Wash- 
ington and northern Idaho if the raw data 
are converted into site index curves. 


The Possibility of Refinement 


The vegetative indicator method, as ap- 
plied to Pinus ponderosa in eastern Wash- 
ington and northern Idaho, allows ex- 
tremely rapid mapping of the land surface 
according to the expected growth rates of 


Height Difference 
in > ton 
index stand 
x] &() 
+) 6 18 
75 65 
+2 64 l 
9 78 
32 69 9 
102 a 
47 66 7 
104 95 
46 74 >] 
Os 90) 
35 74 16 
80 92 
$3 7 15 
79 83 
34 68 15 
96 102 
51 74 28 
91 93 
36 72 a1 
83 101 
5 77 24 


natural stands of pine, and according to 
the hazard of disease caused by Arceutho- 
bum campylopodum. 

Although prediction of site quality as yet 
plays a small role in forestry in the northern 
Rocky Mountains, economics is forcing 
serious consideration of probable growth 
rates before each site is planted. Some 
think a highly refined and therefore accu- 
rate system of site evaluation is not only 
desirable but possible, not only as a guide 
to planting but for interpretation and extra- 
polation of experimental work. In this view 
the four growth-rate categories suggested 
by the vegetative indicator method are 
rather crude. Although clearly different 
from each other, there is much variation 
with the categories. It is entirely possible 
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that more refined synecology might pro- 
vide a more finely divided scale of indicator 
groupings, but there is good reason to doubt 
that vegetation, soil, or any other criterion 
will permit a very accurate prediction of 
growth rates. 

The data in Table 2 suggest that unpre- 
dictable events in the history of a stand 
have marked influence on the average rate 
of height crowth on good and poor sites 
alike. Such factors probably include (1) in- 
tensity of a forest fire as it alters fertility 
from former levels on the same site, (2) 
timing of the fire, consequently the start 
of the new stand, in relation to major cli- 
matic fluctuations, (3) injury due to insects 
or fungi, grazing damage, wind-breakage, 
frost, winter drouth, etc., which the trees 
may experience then later conceal under 
new layers of wood. A given stand may 
happen not to have had its genetically de- 
termined maximum growth rate impaired 
by any of these factors, but at the opposite 
extreme all the factors may have had cumu- 
lative dwarfing the 
These are the variables that must enter 


influence in trees, 
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into the fire-line differences shown in Ta- 
ble 2. The magnitude of these differences 
in growth rates of paired stands provides a 
critical test for any proposed method of 
predicting growth rates. They suggest that 
the growth of the tree on a site provides a 
highly imperfect base for appraising habitat 
potentialities directly or for calibrating 
either vegetation or soils for this purpose. 
Hence, the development of a highly re- 
fined scale of site evaluation seems impos- 
sible of attainment. The most likely ap- 
proach would be to base the curve repre- 
senting full potentialities of the habitat type 
only on the growth rates of the best stands 
found in each age group. 

Within any one habitat type one could 
scarcely question the theoretical existence 
of infinite small steps along multidirec- 
tional gradients in soil and climate which 
could affect average growth rate, and this 
must be added to the influence of other 
components of environment that operate in 
an erratic fashion (e.g. weather extremes, 
biotic factors). It is therefore remarkable 
that any consistent differences among habi- 
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tevo stands 





is determined by analyses of felled 
! line a Pinus 
e trends according to regression analyses. Num- 
able ee 


tats can be demonstrated. Conformity ot 
the vegetational and tree-growth categories 
lends support to the validity and practicality 
of the vegetation classification. 

In Finland the vegetative indicator con- 
cept has been studied for half a century 
ani has become the standard means of 
forest land evaluation. When this system 
was tested in Sweden it was found that the 
indicator significance of a given forest un- 
dergrowth type differed significantly ac- 
cording to altitude. It is possible that some 
of the variation in height growth shown in 
each of the scattergrams in Figure 1 may 
be accounted for by differences in elevation, 
for within one habitat type there was often 
more than 2,000 feet difference among 
some of the stands studied. Also some of 
the Pinus /Symphoricarpos and Pinus/Phy- 
socarpus stands used were in the narrow 
strips of timber that extend from the for- 
ested mountains outward into the grassland 
of the adjacent basal plain. However, the 
dots on the scattergrams representing these 
stands located below average timberline do 
not fall near the lower margin of the pat- 
tern, showing that there is a useful degree 
of homogeneity within a habitat type re- 
gardless of wide differences in macrocli- 
mate. 

Lynch (1958), working in an area that 
includes northern Idaho and eastern Wash- 
ington, has also been concerned about 
“differences in site index estimates as great 
20) units that are often obtained from 
. adjacent stands which, judging on the 


£ 


as 
basis of soil, slope, aspect, and surface vege- 
tation are growing on identical sites.” He 
studied 24 paired stands, the members of 
each pair differing in density but not age, 
presumably a result of partial thinning by 
fire. After detailed biometric analyses, he 
concluded that the differences in rates of 
height growth are a result of intraspecific 
competition associated with different popu- 
lation densities. His conclusions are un- 
usual in forestry literature in that he finds 
low population density of no significance 
for height growth (see Spurr 1952, Bruce 
and Schumacher 1942), whereas extreme 


crowding alone is claimed to reduce rate of 
height growth. 

Since the writer avoided both open and 
very dense stands, his criticisms of Meyer’s 
data as they apply to this region would 
seem to be uncomplicated by the density 
factor. The records kept show that densest 
stands are not represented among the low- 
est dots on the scattergrams, as would be 
the case if they had been dense enough to 
have had height growth impaired. Thus, 
variable tree density would not seem to be 
a measurable component of the height-age 
variation found within each habitat type. 

There are two significant points of agree- 
ment between Lynch’s conclusions and 
those of the present study. His curves for 
stands of average density, as well as my 
curves based on habitat types, show con- 
siderably more diversity in shape than is 
found in Meyer’s site index curves. Sec- 
ondly, when the stem analysis data in 
Figure 2A of this paper are superimposed 
on the curves representing Lynch’s data 
(his Fig. 7D) there is very close agree- 
ment, both sources of data differing from 
Meyer’s curves in recognizing more rapid 
rates of growth between the ages of 25 to 
70 years, at least. 


Biometric Analyses 


As a check on the interpretation of the 
scattergrams regression analyses were kind- 
ly furnished by T. S. Russell of the Com- 
puting Center, Washington State Univer- 
sity. 

Regression lines were constructed of all 
measurements where the average age of the 
stand was less than 110 years, it being 
judged from the pattern of the scatter- 
grams that the age range of 29 to 109 
years would represent largely the relatively 
straight part of a_ theoretically sigmoid 
growth curve. Regression analyses (Fig. 
2B, Table 3) confirmed the order of rating 
based on visual inspection of the scatter- 
grams, and suggest that definitely less than 
seven useful site qualities are indicated by 
the vegetation of seven habitat types. How- 
ever, precise interpretation of such bio- 
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TABLE 3. 


years of age. 


No. of 

Habitat types Stands 
1. Adtes/Pachistima 19 
2. Pseudotsuga/P/ ysocar pus 20 
3, Pinus Physoc ar pus 21 
4. Pinus/Symphoricar pos 32 
5. Pseudotsuga/Calam agrostts 18 
6. Pinus /Purshia 9 
7. Pinus/Agropyron 17 


metric data is not justified when the actual 
age of the trees is never known, and espe- 
cially when, for the sake of objectivity, a 
value of 10 years for growth to breast- 
height is used for all habitat types. Actual 
growth rates of trees, when stratified ac- 
cording to habitat types, do not form a 
series of harmonic lines such as postulated 
by site index curves. Instead the lines for 
A bies/Pachistima and Pinus/Purshia cut 
freely across the others, and all cross the 
site index curves. 

When the regression lines are extra- 
polated backwards toward zero age, most 
of them approximate zero height at zero 
age. Small deviations from the theoretical 
zero height-zero age point are not impor- 
tant, as the growth rates are not uniform 
starting from germination. The outstand- 
ing feature of these extrapolations is the 
high deviation of +36 years at zero age 
for pine growing on the Abdies/Pachistima 
h.t. This regression line as well as the raw 
data in Fig. 1 indicate that for perhaps 
50 years of its life on this habitat the tree 
grows exceptionally well, but soon there- 
after the rate declines sharply. It appears 
to be this juvenile vigor, which is constantly 
thereafter included in the average age of 
the plant, that alone keeps the height-age 
value above that of other habitats except 
late in the life of the tree. For a 50-year 
timber rotation the 4 bies/Pachistima habi- 
tat is clearly the most productive, but for 
a 100-year rotation three other habitats are 
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Summary of biometric analysis of data for stands less than 110 


95 percent 
Confidence 
Interval (ft./yr.) 


Regr. Equation 


Crt See 


¥ 35.99 + 0.58 X 0.28 — 0.88 
= 3.31 + 1.11 X 0.72 1.50 
x 102+ 1.07 X 0.70 1.44 
Y 6.59 + 0.91 X 0.61 1.21 
7 3.70 + 0:89 ZX 0.72 — 1.06 
Y 1.33 + 0.70 X 0.35 1.05 
¥ 5.10 + 0.85 X 0.55 1.15 


superior (Fig, 2B). 

On the Thuja-Tsuga/Pachistima habitat 
type, which is still cooler and wetter than 
the Abies/Pachistima habitat type, the 
growth rate of planted P. ponderosa (al- 
most none grows there naturally) declines 
so strongly in a few years (Daubenmire 
1952) that plantations are complete failures 
on such sites. There is thus vood evidence 
that, in the region under study, the growth 
of P. ponderosa begins to suffer on sites 
farther up the moisture-temperature gradi- 
ent than the Pseudotsuga/Physocarpus 
habitat type, with decline starting earlier 
and being more marked the farther up the 
sequence, 


Key to Habitat Types Studied 


Most forest stands encountered in the field 
have been burned or grazed, and so sup- 
port an admixture of species indicative of 
disturbance. Yet the valuable indicators re- 
main for those who know what to look for. 
The oft-expressed notion held by many 
American foresters that lower vegetation 
is of no value on disturbed lands is not 
founded on detailed study, but upon im- 
pressions resulting from the great differ- 
ences in physiognomy that are induced by 
disturbance, but which are really super- 
ficial and can be treated accordingly after 
the synecology of a forest mosaic has been 
worked out. A short key can be con- 
structed to the seven habitat types consid- 
ered in this study making use of only a 


() 


few of the most distinctive botanical char- 


icters of each habitat, which usually are 


pparent despite any disturbance. For 
further description of the habitats see Dau- 
nmire 1952. 
The habitat 


ifter the species dominant 


types have been named 


under climax 
conditions, so that the names do not neces- 
reflect each 
site. “hus south of the limits of Pleistocene 


gl mation 


sarily present dominants of 


stand of Pinus ponderosa with 


an understory of Ceanothus sanguineus 


plus a liberal sprinkling of young Pseudo- 


would 


tsuga certainly within the 


occur 
environment range, and exhibit the growth 
characteristics of a Pseudotsuga/Physocar- 
pus habitat type. Therefore names are no 


more absolutely 


diagnostic in synecology 


than in taxonomy, where a Liriondendron 
tree may he called tulipifera for years hbe- 
Al les ari- 


fore it produces a flower, and 


zonica can occur in Colorado, 
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Summary 


The rate of height growth of Pinus pon- 
derosa has been studied in 153 stands in 
eastern Washington and northern Idaho. 
When the data are grouped according to 
habitats, these based on the plant associa- 
tions that are climax for the sites, good 
correlations are obtained, so that vegetation, 
whether climax or not, serves as a rapid 
and useful means of predicting in advance 
the probable growth rate of the pine. Sus- 
ceptibility of the tree to infection by Arceu- 
thobium is also correlated with these vege- 
tative indicators. The data suggest that the 
site index curves that have been used for 
Pinus ponderosa over a wide area are not 
valid in the region under study. A key is 
provided for the identification of the seven 
habitat types supporting natural stands of 
Pinus ponderosa in the region studied. 
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Crown Competition—A Measure of Density 


‘THE MEASUREMENT of stand density has 
been one of the most trying problems in 
forestry. Many measures have been pro- 
posed but a completely satisfactory one has 
not yet been found,’ In an effort to help 
solve the problem, a new measure of stand 
density, developed at the Central States 
Forest Experiment Station, is presented 
here. 

Density, as foresters commonly use the 
term, is the relation between the number 
of trees or some volumetric or areal unit 
In effect, 
then, a density measure is intended to help 
determine the relation of the average tree 


to a specific area, usually | acre. 


in the stand to the maximum growing space 
it could utilize on the one hand, and the 
minimum growing space necessary to live 
on the other. It is within this range that 
the forester wishes to maintain density to 
fully utilize the site for maximum produc- 
tion of desired usable volume. 

The maximum growing space that can’ 


'The authors recognize the vast amount of 


literature on the subject of stocking and 
lensitv, but believe a review of this work 
Ww uld add very little to the presentation of 
new measure given in this report. Good 
ussions of stocking and density measures are 
given by: 


Spurr, Stephen H. Forest inventory. The 
Ronald Press Co. New York. 476 pp. 1952 
Bickford, C. Allen, et al. Stocking, nor- 

lity, and measurements of stand densit 
(Report of S.A.F. committee on stocking). 
J. For. 55:99-104. 1957. 
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be used by a tree of small diameter is, of 
course, less than that for a tree of larger 
diameter. The problem of determining 
the growing-space/tree-size_ relations _ is 
complicated by the effects of present and 
past density on both stem diameter and 
crown size. So for a given species of speci- 
fied d.b.h., the size of the crown may vary 
greatly. ( Although trees compete for mois- 
ture and nutrients as well as for light, it 
would be impractical to measure and re- 
late size of root systems to degree of com- 
petition between trees. In fact, it is prob- 
able that crown area may be an indicator 
of root area and thus a measure of such 
subterranean competition.” ) The only place 
where the crown d.b.h. relations are not 
confounded by competition from other trees 
is in the open-grown tree. So the new 
method described here is based on the rela- 
tion of crown area to d.b.h. of open-grown 
trees.” 


siisgen, M. and Miinch, E. The structure 


€ 
and life of forest trees, page 273. New York, 


\ brief presentation of this subj ct is riven 
in Station Note 108 “Crown development: an 
index of stand density” by Krajicek and 
srinkman. Central States Forest Experiment 
Station, U. S. Department of Agricultur 


For. 


The authors are respectively Research For- 
ester, Silviculturist, and Research Forester, 
Central States Forest Experiment Station, For- 
est Service, U. S. Department of Agriculture. 
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Development of Crown-Width/D.B.H. 
Relations 


To establish the relation between width of 
crown and d.b.h. for open-grown trees, 
34) open-grown trees were measured in 
eastern Lowa: 88 white oaks (Quercus alba 
L.), 60 black and red oaks (Q. velutina 
Lam, and Q. rubra L.),* 35  hickories 
(Carya ovata (Muill.) K. Koch), and 157 
Norway spruce (Picea abies (L.) Karst.). 
As nearly as could be determined, each tree 
had developed without competition through- 
out its entire life. “To maintain uniformity 
in selecting sample trees, only those that 
met the following specifications were in- 
cluded in the sample: 





1. Crown free of competiticn on all sides. 
2. Limbs extending to the ground on small 
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No evidence of pruning, shearing, brows- 
ing, decadence, storm damage, or seri- 
ous insect damage. 

6. Tree apparently not of sprout origin. 

The smaller hardwood sample _ trees 
were found in old fields, on road rights-of- 
way, and on recently cutover, ungrazed 
areas where the sample trees originated 
after cutting. The larger hardwood sampk 
trees were found in’ timbered pastures, 
cemeteries, and, to a limited extent, in 
residential areas. Norway spruce sampk 
trees were found principally in cemeteries 
and parks located in or near small towns, 
in recently established windbreaks, and 
occasionally in residential areas, 

Crown width of each sample tree was 
measured twice, the second measurement 
at right angles to the first, and the results 
averaged. Regression analyses of the data 
showed definitely that the crown width of 
an open-grown tree is closely related to its 


d.b.h. 

Site quality was not determined directly 
for the locations where the sample trees 
were found. For the species examined, 
open-grown trees were somewhat shorter 
than their forest-grown counterparts of 
the same diameter on similar soils and as- 
pects. Most of the hardwood sample trees 
were found on forest soils where site in- 
dices ranged from 45 to 75. The rest were 
found on prairie or transitional soils in 
which the productivity of the similar for- 
est soils was known to range from poor to 
evood, 

Differences in the relation between 
crown and stem diameter among the hard- 
wood species sampled were small (Fig. 1). 
The differences between the red oaks and 
white oak and between white oak and hick- 
ory were not statistically significant at the 
5-percent level. The difference between 
hickory and red oaks was significant, but 
just barely so. For simplicity in applica- 
tion all hardwood data were combined 
(Fig: 2). 

Crown width and d.b.h. for Norway 
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spruce also were highly correlated (Fig. 
3.). The differences between the hard- 
woods and Norway spruce were highly 
significant at the l-percent level. 

Since such a close relationship exists be- 
tween the crown diameter and d.b.h. of 
open-grown trees, and this relationship is 
nearly constant within a species, it may be 
inferred that the crown of a tree of given 
d.b.h. cannot occupy more than a certain 
area even with unlimited growing space. 


Development of an Index of 
Crown Competition 


Knowing the average crown width and 
assuming that crowns of open-grown trees 
tend to be circular, the approximate area 
covered by the crown can be readily com- 
puted. The percentage of an acre occupied 
by a vertical projection of the crown is 
obtained by dividing the area in square feet 
by 435.6. For example, an open-grown 
Norway spruce with d.b.h. of 16 inches has 
a crown width of approximately 26 feet 
(Fig. 3), or crown area of 531 square feet 
(1.22 percent of an acre). This value has 
been termed 


” 


“maximum crown area 
(MCA), because it expresses, im terms of 
percent of an acre, the maximum area that 
could be occupied by the crown of a tree 
of specified d.b.h. 


‘To determine MCA values for each tree 
in a stand would be time consuming. How- 
ever, from the least squares solution of 
crawn width and d.b.h., a formula can be 
derived that simplifies the procedure: 

The MCA value of a tree with a crown 
diameter of CW is 

x (CW)°- 
MCA 
4 
435.6 
0.0018 (CIV )* 
From the regression analysis of the com- 
bined oak-hickory data (Fig. 


as 
CW (1.829D 3.12) 
where D d.b.h. 
and 
(CW )* 3.345 D- 11.413D + 9.734 


Therefore, 
MCA 0.0060 D=- + 0.0205D 


2 0.0175 
Similarly, the formula of MCA values for 
Norway spruce is 


MCA 0.0031 D= + 0.0239D + 0.0460 
With the MCA values determined, it is 


not difficult to compute the sii num- 
ber of open-grown trees per acre that are 


} 
yroduce a com- 


theoretically required to pl 

plete canopy (assuming that crowns ar 
distorted sufficiently to fill all voids in the 
canopy, but without reducing the crown 
area of the individual trees), If a tree has 
a d.b.h. of 16 inches and an MCA of 1.22 


(as in the previous example ), then 82 open- 
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TABLE 1. Determination of the CCF for an even-aged upland oak stand, site 


index 60, age 50 (Schnur 1937). 


D N DN 
(Number of trees 

(d.b.h.) per acre) 
1 10 10 
Z 29 58 
3 50 150 
4 63 252 
5 I55 
6 78 +68 
7 66 462 
8 49 39? 
9 31 279 
10 ws 170 
11 7 77 
12 3 36 
13 1 13 
Total 481 2,752 


CCF = 0.0060(18,346) + 0.0205(2752) + 0.0175(481) = 174.91 


grown trees of that size could cover an 
acre completely (82 X 1.22 = 100 per- 
cent, or sum of MCA values equals 100). 
Using MCA values for other d.b.h. classes, 
any combination of open-grown trees of 
various diameters whose MCA values total 
100 for an acre theoretically could have a 
closed canopy. 

‘To avoid confusion in using one term— 
MC A—to express the area requirement of 
both a single tree and a sum of trees, an- 
other expression, “Crown Competition 
Factor” (CCF), has been coined to rep- 
resent the sums of the MCA values for an 
acre. Crown Competition Factor is defined 
as the sum of MCA values for all trees in 
the stand, divided by the area in acres. It 
is used as an expression of stand density. 
These summations can be expressed by us- 
ing the previous formulas in the following 
form: 


For oak-hickory stands, 
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D=N Check using MCA values 


MCA per tree Total MCA 


Percent Percent 

10 0.044 0.440 
116 082 2.378 
450 .133 6.650 
1,008 .196 12.348 
1,925 .270 20.790 
2,808 356 27.768 
3,234 455 30.030 
3,136 2566 27.734 
2,511 688 21.328 
1,700 822 13.974 
847 .969 6.783 
432 1.128 3.384 
169 1.298 1.298 
18,346 174.905 
(CCF) 


CCF = — 
A 
[0.0060(2D7Ni) + 0.0205(2DiNi) 
+ (.0175(2Ni)] 
and for Norway spruce, 
l 
CCF — 
A 
[0.0031(2DF7Ni) + 0.0239(2DiNi) 
+ 0.0460(2Ni) ] 


where 
Dj; = individual d.b.h. or d.b.h, class 
Ni = number of trees in d.b.h. class 
A = area in acres 


An example of the application of the 
CCF to a specific stand is given in Table 1. 
All the trees in the stand are used to de- 
termine the CCF. In the stand of trees 
shown in Table 1, nearly 60 percent of the 
Crown Competition Factor is contained in 
trees less than 8 inches d.b.h. To consider 


Ol 


lo 


only those trees larger than a specified 
lower diameter limit will result in a low 
CCF and will not represent the actual 
competition in the stand. It is known, of 
course, that the larger trees have an ad- 
verse effect on the growth rate of the 
smaller. It is more difficult to assign an 
effect of the small trees on the growth of 
larger ones. Where moisture is a limiting 
factor at any time during the growing sea- 
son, the smaller trees probably have an 
adverse effect on the growth of the larger 
trees. 

One other precaution is necessary. When 
computing the CCF from stand tables, the 
actual number of trees in each diameter 
class should be used rather than the diame- 
ter of the tree of average basal area and 
total number of trees. If the latter proce- 
dure is followed, excessive errors may re- 
sult when the diameter distribution deviates 
greatly from the normal curve. 


Results and Discussion 


It is obvious that maintaining stand density 
below CCF 100 from time of establish- 
ment to maturity would normally produce 
low-quality trees, especially if they were not 
pruned artificially. Greater density (with a 
higher CCF value) would be accompanied 
by a reduction of the actual or effective 
live crown of the average tree in the stand. 

Density will not increase indefinitely, of 
course, A number of even-aged pure oak 
stands in Iowa consistently showed a ceil- 
ing density of about CCF 200, when all 
trees in the stand were included. Similarly, 
the average CCF for 52 even-aged oak 
plots, each 0.2 acre in size, located through- 
out southeastern Ohio and selected for full 
stocking was exactly 200.° A consistent 
ceiling probably exists for all other species 
also, but it is probably different from that 
found for oaks. 


°These unpublished data on which the com- 
putations are based were collected and fur- 
nished by Willard Carmean and Stephen G. 
Boyce, Soil Scientist and Research Forester, 
respectively, Central States Forest Experiment 
Station. 


At least three things account for differ- 
ences among species in the maximum CCF 
they can obtain. One is their crown area 
development when there is no competition. 
Another is the basic shape of the crown, 
such as conical, parabolic, etc. Also, toler- 
ance may be a factor. For example, meas- 
urements over an 8-year period in heavily 
stocked oak-hickory plots show a slightly 
higher CCF than those found in pure oak 
stands, probably because of the greater tol- 
erance of hickory. Some of the Ohio plots 
cited above also had CCF values in eX=- 
cess of 200, but these too contained some 
of the more tolerant species. 

To further confirm the consistency of 
the values found in field tests, the formula 
for oak-hickory was applied to stand tables 
for even-aged upland oaks. Except for 
high values in younger stands on the poorer 
sites, the CCF values show great uniform- 
ity (Table 2), 

No Norway spruce stands were avail- 
able for field testing of the formula. The 
similarity of crown widths among the oak 
species suggested that perhaps a given 
crown width/d.b.h. relation may be rather 
common within a genus. With this thought, 
the Norway spruce formula was applied to 
stand table data for Sitka spruce and west- 


TABLE 2. Application of CCF to 
stand table data’ for even-aged up- 
land oaks. 


(CCF values by age class and site index) 


Age Site index 
(vears) 40 50 60 70 80 
10 301 249 219 190 166 
20 218 210 188 180 181 
30 184 183 178 176 167 
40 171 170 168 165 160 
50 169 169 175 173 173 
60 171 171 176 175 169 
70 165 168 172 170 165 
80 168 176 173 172 164 
90 175 179 177 174 164 
100 178 182 177 175 166 
1Schnur, G. Luther. Yield, stand, and volume tables 
for even-aged upland oak forests. U. S. Dept. Agric. 
Tech. Bull. No. 560, Table 37, 1937. 


volume 7, number 1, 1961 / 39 








ern hemlock (Table 3). These data are 
primarily for stands in Washington and 
Oregon, yet the CCF values for the age- 
site classes enclosed by the broken line are 
quite consistent, except for age 140, site in- 
dex 160. The latter may be an error in the 
original table. 

A comparison of oak-hickory and Nor- 
way spruce stands of the same density as 
measured by crown competition reveals a 
For 
example, a Norway spruce stand composed 
of 16-inch 
half again the 


striking difference in their structure. 


more than 

number of trees and_ basal 
as an oak-hickory stand with 
trees of like size and the same CCF. 


trees would have 


area per acre 


Comparison With Other Measures 


A thorough comparison of CCF with other 
measures of stand density would require 
much space. A brief comparison serves to 
illustrate the relation between CCF, basal 
area per acre, and Reineke’s stand-density 
index (Table 4). 
that when CCF 
of diameters, basal 
area per acre changes at a rate that is 
typical of most even-aged stands, CCF 
not a consistent relationship 


It is of interest to note 
is held constant over a 
stand 


ran ge 


average 


does show 


TABLE 3. 


with stand-density index. Although a CCF 
of 100 (which might indicate minimum 
density for full site utilization) is approxi- 
mately associated with a stand-density in- 
dex of 110 percent for stands with average 
diameters as shown, the relationship of 
number of trees to average stand diameter 
is curvilinear on log-log paper for any given 
CCF value, whereas stand-density index is 
linear. However, this does not represent a 
deterrent to the application of CCF. Ap- 
proximate full stocking in terms of stand- 
density index (240) results in CCF values 
from 197 to 262 as shown in column 6 of 
Table 4. Again, one could conclude im 
general that a CCF of 200 (which might 
represent maximum for full site 
utilization) is approximately equivalent to 
a stand density 240 for stands 
with an average diameter of from 5 to 15 
inches. 

In comparing “D plus” and “D times” 
spacing® with CCF spacing 


density 


index of 


it is apparent 


®Rules of thumb used as thinning guides: 
Using the “D plus” rule, spacing of trees in 
feet is equal to the diameter of the trees in 
inches plus a constant; using the “D times” 
rule, spacing in feet is diame- 


equal to tree 
1 


ter in inches times a constant. 


Application of CCF to stand table data’ for even-aged Sitka spruce 


and western hemlock. (CCF values by age class and site index.) 





Age Site index 

(vears) 60 80 100 120 140 160 180 200 
40 1184 63 412 345 320 | 292 = 275 269 | 
60 427 351 332 299 288 276 270 267 | 
80 356 327 306 293 284 283 280 276 
100 326 299 296 288 284 285 285 285 | 
120 304 296 291 286 285 85 287 289 
140 29] 289 286 286 286 241 | 290 287 | 
160 286 285 287 287 287 291 =. 296 
180 283 282 281 283 288 288 292 

200 279 280 280 284 289 293 294 

Meyer, W Yield of even-aged stands of Sitka spruce and western hemlocks. U. S. Dept. Agric. Tech. 

Bull. No. 544, Table 36, 1937 
Dash lines indicate limits of basic data. 
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win 


TABLE 4. 


theoretical oak-hickory stands of various average diameters. 








Comparisons of CCF with basal area and stand-density index" for 


When stand-densits Basal area e 
When CCF is 100: index is 240: when: 
Stand- 
l'rees per Basal area density Trees per Equivalent CCF = SDI 
acr per acre index icre CCF 200 240 
\ i??, ber Square tee Per ‘ent Num ér Per ‘ent 5 Juare feet Sguare feet 
2 1,220 26.8 92 3,177 262 S37 69 
5 370 50.3 122 730 197 101 100 
10 122 66.5 122 240 197 132 131 
15 60 73.6 115 125 209 146 153 
20 35.3 77.0 107 78.9 na 154 172 
25 23.4 79.8 102 55.1 236 159 188 
1Reineke, L. H. Perfecting a stand-density index for J. Ag R 46: 627-638 3 
“Calculated from average basal area per tree. (Used in the example to facilitate ¢ \ ind- 
nsity index. As stated in the text, tl actual number of tree ameter cla } ite 


TABLE 5. Spacing of trees for a CCF of 100 and equivalent “D plus” a 
“D times” constants for theoretical stands of oak-hickory and Norway spruce 
various average diameters. 


Oak-hicko Norway spruce 
Average d.b.h.! Cc ‘D plus “D times” CCF “T) plus” “T) times’ 
(inches) spacing constant constant spacing constant onstant 


Feet 


_ 
“~ 
— 


Feet Unit 


na 
of 


2 6.0 4.0 3.00 6.8 4.8 3.40 

5 10.8 5.8 2.16 10.3 5.3 2.06 

10 18.9 8.9 1.89 16.1 6.1 1.61 

15 26.9 11.9 1.79 21.9 6.9 1.46 

20 35.1 15.1 1.76 oral eal 1.38 

25 43.1 18.1 1.72 33.6 8.6 1.34 
1Calculated from average basal area per tre 


that spacing based on CCF will not result 
in a constant factor for stands of different 
average diameters (Table 5). An examina- 
tion of these substantiates the 
findings of Stahelin’ and others that diame- 
ter plus a constant results in too much basal 
area with large trees and too little with 
small trees. 


“constants” 


‘Stahelin, R. 


and slash pine stands to specified densities. 


J. For. 47:538-540, illus. 1949. 


Thinning even-aged loblolly; 


It should be emphasized that CCF is not 
essentially a measure of crown closure. 
Theoretically, complete crown closure can 
occur from CCF 100 to the maximum for 
the species (in oaks, approximately CCF 
200). Instead of estimating crown closure, 
CCF estimates the area available to the 
average tree in the stand in relation to the 
area it could use, if it 


maximum were 


open grown, 
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Conclusions 

There is good evidence to indicate that the 
CCF method of measuring density can be 
applied to even-aged stands of hardwoods 
and conifers. Moreover, although field 
tests of the method have been confined to 
even-aged stands, in theory it can be used 
for uneven-aged stands equally well. The 
method adjusts accurately to variations in 
diameter distribution and to the potential 
crown space requirements for the indi- 
vidual trees. 

Traditionally, researchers have _ insisted 
that the ideal measure of density should be 
correlated with growth and yield and inde- 
pendent of stand age. On this premise CCF 
has shown up well in preliminary tests. In 
stand density studies in even-aged oak net 
basal area growth (accretion less mortal- 
ity) is at a maximum when CCF is 135. 
This holds true for stands 25, 50, and 90 
years old, with corresponding average 
basal area per acre of 71, 86, and 96 square 
feet respectively. The growth data cover 
an 8-year period for the two younger age 
classes and a 4-year period for the other. 

Several important points should be em- 
phasized: 

1. Relatively few measurements of 
truly open-grown trees are required to pro- 
vide a basic CCF equation for a species. 

2. The ceiling of about CCF 200 for 
oak density should not be assumed for other 
species, even when crown widths are found 
to be similar. Relative tolerance and shape 
of crown (conical, parabolic, etc.) deter- 
mine the maximum densities attainable by 
a species, 

3. In stands having mixtures of species 
with different crown characteristics and 
tolerance, the determination of CCF may 
need to be handled by proportional repre- 
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sentation. 

The CCF is as simple to apply as basal- 
area tables, but, unlike basal-area control, 
interpretation of values does not depend 
upon diameter distribution of a_ specific 
type or range. “The Crown Competition 
Factor, based on a verified relationship be- 
tween crown width and d.b.h., should 
Measure competition more acurately than 
the methods presently used. 


Summary 


Measurements of several hundred oak, 
hickory, and Norway spruce disclosed a 
high degree of correlation between crown 
width and d.b.h. for open-grown trees of 
the same species. Apparently this relation is 
independent of age and site quality. From 
the open-grown-tree data, formulae were 
derived which express the maximum per- 
centage of an acre the crown of an open- 
grown tree of a given diameter and species 
could occupy. 

Theoretically, competition for crown 
space between trees in a stand begins when 
all crown space is occupied and each tree 
crown is equal in area to that of an open- 
grown tree of the same d.b.h. This condi- 
tion is taken as a base of 100 for expressing 
competition among crowns for growing 
space. From this base, stand density may 
be expressed as a percentage. For con- 
venience, this percentage is termed the 
Crown Competition Factor (CCF), be- 
cause it is an estimate of the competition 
among crowns for growing space. Pre- 
liminary trials indicate that the method is 
applicable to even-aged stands, and theo- 
retically can be used to advantage in 
uneven-aged stands. ‘These trials show 
strong evidence that site and stand age do 


not influence CCF. 


The Influence of Gibberellic Acid on Growth 
and Fiber Length of Quaking Aspen 


GIBBERELLIC ACID, the metabolic product 
of the fungus Gibderella fujikuroi, is a 
growth-promoting chemical that has been 
tried widely on agricultural crops and re- 
cently has been the subject of a number 
of studies having forestry implications. 

Some of the earliest tests of this chemical 
on forest tree species were made by Marth, 
et al. (1956). Their work produced stem 
growth increases of from 0 to 400 percent, 
with the conifers showing the least response 
to the treatments. Since this first work, 
additional information has been published 
regarding the influence of gibberellic acid 
on tree growth by Nelson (1957), Hac- 
skaylo and Murphey (1958), Knight 
(1958a, b), Schaedle (1958), Swan 
(1958), and Westing (1959). Most 
trials have been conducted on young trees 
with the practical aims of improving early 
growth and development. Results ob- 
tained from these studies range from very 
promising to disappointing, and are quite 
variable because of the number of tree 
species inyolved and testing methods em- 
ployed. 


Preliminary Trial 


The present study was established to inves- 
tigate the influence of gibberellic acid on 
growth and fiber properties of quaking as- 
pen (Populus tremuloides Michx.). The 
study was divided into two phases. The 
first was a greenhouse trial of limited scope 
designed to give preliminary information 
regarding methods of application, treat- 
ment levels, response time, and the fre- 
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quencies of reapplication necessary to main- 
tain rapid growth whe: gibberellic acid 
was applied as a foliage spray. The second 
phase was a field trial using the levels of 
chemical established in the preliminary 
work, 

Quaking aspen seed from an _ experi- 
mental cross made in the greenhouse was 
planted in 11 8-inch standard clay pots. 
After five weeks each pot was thinned to 
EZ good seedlings and at eight weeks of 
age the treatments listed in Table 1 were 
applied. ‘This phase of the study was not 
replicated. The chemical’ used consisted of 
a combination of gibberellic acids A and 
X. 

The 500 ppm, 100 ppm, and the 1 ppm 
treatments were applied as an aqueous spray 
using Tween 2() to improve coverage and 
uptake. The spray treatments were reap- 
plied at 2, 4, and 6 week intervals. Spray 
treatments consisted of a 10-ml solution 

1The gibberellic acid used was supplied by 
Abbott Laboratories, North Chicago, Illinois. 
Additional information regarding chemical 
composition and methods of handling can be 
obtained by writing to that company. 


The authors are Silviculturist and Geneti- 
cist, respectively, at The Institute of Paper 
Chemistry, Appleton, Wisconsin. They are 
greatly indebted to Abbott Laboratories. 
North Chicago, Illinois, for financial support 
of the study and to Marlin T. Leffler and 
Norman Leffler of that company for technical 
assistance in handling and applying gibberellic 
acid. Manuscript received Mar. 14, 1960. 
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TABLE 1. Summary of height 
growth response of aspen seedlings 
treated with gibberellic acid in the 
preliminary trial. 


Frequency Height growth 


of reappli- response (treat- 
G.A. treatment cation ment ws. control) 
Weeks Percent 
500 ppm 2 216 
4 451 
6 62 
100 ppm 2 214 
4 146 
6 98 
10 ppm 2 —14 
4 —14 
6 40 
1% G.A. in 
Lanolin - 110 


applied to each pot of 12 seedlings as a 
fine mist, using a small electric paint 
sprayer. The 1 percent gibberellic acid in 
lanolin paste was applied by encircling the 
seedling stem with a thin band about %4 
inch above the ground line. A toothpick 
was used for making the application. 

The seedlings in this trial were fertilized 
weekly with a soluble fertilizer and grown 
in the greenhouse on a 14-hour day. Meas- 
urements and observations on trcated seed- 
lings were made twice each week and 
treatments were reapplied as scheduled. 


Results of the Preliminary Trial 


Table 1 summarizes the height growth re- 
sponse obtained when the treated seedlings 
were compared with the control seedlings 
eight weeks after treatment. Figure 1 
illustrates the type of growth curves ob- 
tained from height measurements and Fig- 
ure 2 the growth response and the vari- 
ability of the treated individuals in several 
of the applied treatments. 

Little or no response was obtained with 
the 1 ppm level at any of the treatment 
reapplication frequencies. The 1 percent 
gibberellic acid in lanolin paste gave a re- 
sponse similar to the 100 ppm and the 500 
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ppm treatment levels. The best treatment 
in this preliminary trial appeared to be the 
500 ppm applied every four weeks. 
Reapplication every six weeks resulted 
in the least response, both at the 100 ppm 
and the 500 ppm levels. Reapplication 
every two weeks at 100 ppm and 500 ppm 
levels resulted in a number of individuals 
in each pot with long, narrow leaves, light- 
green leaf color, and spindly, crooked, ab- 
normally Response after 
treatment was rapid and changes in the 


formed stems. 


slope of the growth curves usually were 
evident within four to six days after appli- 
cations. 


Field Trials of Gibberellic Acid 


The second phase of this study involved 
the application of the results obtained in 
the preliminary trial to a study aimed par- 
ticularly at measuring the influence of gib- 
berellic acid on growth and on fiber length 
of quaking aspen under field conditions. 
Fiber length was of special interest because 


MM 


HEIGHT, 


TOTAL 


Figure 1. Height growth curves of green- 
house-grown quaking aspen treated with a 
100 ppm aqueous gibberellic acid spray. 
The several curves represent seedlings treat- 
ed every two weeks, every four weeks, 
every six weeks, and the control. The ar- 
rows indicate the time the treatments were 


reapplied, 








—— 


—— 


nt 
the 


gt 


of its importance in influencing certain 
strength properties of the paper produced 
from aspen. Also of interest in this second 
phase was the further refinement of treat- 
ment levels and frequencies of reapplica- 
tion, 

Quaking aspen seed from a controlled 
experimental greenhouse cross was again 
used as a source of experimental material. 
The seed was sown early in March and 
the seedlings transplanted into plant bands 
at four weeks of age. At 11 weeks, the 
treatments listed in Table 2 were applied 
to a total of 1,760 trees, and the trees 
then planted in the field. The gibberellic 
acid treatment levels employed were 250 
ppm, 100 ppm, and 10 ppm applied as an 
aqueous spray and 1 percent gibberellic 
acid in a lanolin paste. The reapplication 


ce 


frequencies used were every 3, 4, and 5 
weeks, The treatments were applied as 
previously described, the aqueous spray with 
an electric paint sprayer and the gibberellic 
acid-lanolin paste as a small band approxi- 
mately 1 inch above the ground line. Re- 
application in the field was made on sched- 
uled days with enough spray used per in- 
dividual to wet the leaves thoroughly. The 
gibberellic acid-lanolin paste treatment was 
not reapplied. 

Each treatment was replicated twice, 
and each replication contained 80 seedlings 
lined out in a single row spaced within 
the row at 8 to 10 inches. Weekly height 
measurements were made on the first 20 
individuals in each replication and addi- 
tional observations (survival, leaf size and 
color, and dieback) were made periodic- 
ally during the growing season on all 80 
seedlings. Some information regarding the 
effect of gibberellic acid on root develop- 
ment was obtained by oven drying and 
weighing the root systems from the first 
ten individuals from each of the replica- 
tions. 

The time-consuming aspects of fiber 
measurement work made it desirable to 
check only certain representative treat- 
ments. Two trees of average size from 
each replication or a total of 4 trees from 





Figure 2. Growth response of aspen seedlings 


treated with gibberellic acid in the prelim- 
inary greenhouse trial, 8 weeks after treat- 
ment. These aqueous spray treatments were 
applied every four weeks. 


each selected treatment were sampled by 
taking a 6- to 8-inch section from the lower 
portion of the stem. The samples were sub- 
sequently chlorited according to the meth- 
od of Spearin and Isenberg (1947), and 
stained with congo red. The stained fibers 
were washed, suspended in the proper con- 
centration and mounted on microscope 
slides. From each sample 100 fibers 
(tracheids) were measured on a projector 
at a magnification of 75X. 
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Results of Field Trial with 
Gibberellic Acid 


Weekly height-growth measurements were 
made during the growing season and a 
summary of the measurements is presented 
in Tables 2 and 3. Several treatments are 
further illustrated in Figures 3, 4, and 5. 
Table 2 also provides information regard- 
ing survival, percentage of trees showing 
tip dieback,” diameter growth and root de- 
velopment of the treated trees. 

Analyses of variance were made on the 
original data summarized in Table 2. Re- 
sults of the analyses presented in Table 4 
are useful in interpreting treatment effects. 
The “F” 
using the error mean square term except 
in those instances where the level x fre- 
quency (L x F) interaction was significant. 
When the L x F 


cant the “F” 


values listed were computed 


interaction was signifi- 
values for levels and fre- 


“The dieback of the growing tip appears 
similar to a disease caused by 


t 


N a piclad 1u0n 
tremulae in which young shoots bend over, 
darken in color, and dry out. 


TABLE 2. 
berellic acid in replicated field trial, 


Frequency 


of reappli- Trees 
cation planted Survival 
Treatment (weeks) (No.) (percent) 
Gibberellic 
acid 
250 ppm 3 158 44.9 
+ 155 S22 
5 157 69.4 
100 ppm 3 159 79.2 
+ 58 89.2 
5 159 88.7 
10 ppm 3 159 91.8 
+ 148 98.1 
5 156 96.8 
] percent 
G.A,. in 
lanolin - 158 94.9 
Control - 159 92.4 


1Percent of trees showing growing tip dieback 
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quencies were computed using the L x F 
interaction mean square term. 

Height growth of the treated individ- 
uals in the field trials during the first 8 
weeks followed a pattern quite similar, al- 
though much less spectacular, than that 
of the preliminary greenhouse trial, At 
8 weeks, as Table 3 illustrates, height 
growth response of the treated trees when 
compared with the control trees ranged 
from 3.4 to 48.9 percent, with the 100 
ppm and the 250 ppm treatments showing 
the greatest response. During the growth 
period from the 8th to the 14th week, the 
growth rate of the control trees was quite 
rapid and at the end of the 16th week, the 
average height of the control trees was ap- 
proximately equal or greater than any 
treatment, except the 100 ppm five-week 
treatment. 

The “F” values for the average height 
data indicate no significant differences in 
seedling growth between control and spray 
treatments, lanolin and spray treatments, 
and height growth differences resulting 
from levels or frequencies of reapplication 


Growth, survival, and dieback of aspen seedlings treated with gib- 


Ave. Stem 
Dieback! Ave. ht. Dia. Ave. o.d. wt 
(percent) (cm) (mm) roots (g) 
14.6 54.5 33 0.9 
15.5 56.8 .28 0.5 
11.5 64.1 43 a9 
13.2 60.9 34 1.0 
11.4 52.0 36 ye 
10.1 78.7 .67 7.4 
9.4 64.2 oF 7.6 
2:5 59.7 56 6.3 
1.9 60.6 37 Tan 


.0 49.1 47 4.7 
1.9 62.7 OF 8.9 


on 7/14/58, after a period of rapid growth. 


ae 


rrr 6 





TABLE 3. 
end of o and 16 weeks. 


Frequency of 


reapplication Ave. ht. 
‘Treatment (weeks) (cm) 
Gibberellic acid 
250 ppm 3 36.8 
+ 39.6 
5 33.0 
100 ppm 3 36.5 
+ 33.4 
5 35.6 
10 ppm 3 32.4 
4 30.0 
5 29.0 
1% G.A, in 
lanolin ~ 27.5 
Control = 26.6 


Ip 


of gibberellic acid. A separate analysis using 
data from the control and 100 ppm 5- 
week spray treatments only gave a highly 
significant “F”’ test. Since in this case the 
difference includes the effect of the treat- 
ment and the influence of soil differences, 


HEIGHT GROWTH OF GIBBERELLIC 
ACID TREATED ASPEN 





CONTROL VS. 10 PPM 


AVERAGE HEIGHT, CM 





0 2 4 © = 3 i2 14 16 


TIME SINCE TREATMENT, WEEK 


Figure 3. Height growth of quaking aspen 
treated with 10 ppm gibberellic acid, based 
on measurements of seedlings in the repli- 
cated field trial. 


At 8 weeks 


ercent height growth increase or decrease, treated seedlings ¢ 


Average height and response of aspen seedlings to treatments at the 
é & / e 


At 16 weeks 


Response! Ave. ht. Response! 
(percent) (cm) (percent) 
38.3 54.4 —13.2 
48.9 56.8 - 9.4 
24.1 64.1 + 2.2 
37.2 60.9 2.9 
25.6 52.0 17.1 
33 s 78.7 +25 5 
21.8 64.2 + 2.4 
12.8 59.7 4.8 
9.0 60.6 3.4 
3.4 49.1 —21.7 
62.7 


ompared with control seedlings. 


one should be cautious in interpretation. 
The differences, however, are of sufficient 
magnitude to warrant further research 
using similar concentrations and frequencies 
of application. 


Diameter growth was not checked as 
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Ficure 4. Height growth of quaking aspen 


treated with 100 ppm gibberellic acid, 
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FiGurE 5, 


Height growth of quaking aspen 


4 


ees aa oe oe a ee sy 
tréatéd wit! gibberellic actd, 


closely as height growth and only stem 
diameter at the ground line was measured 
at the end of the growing season. General 


field observations indicated that 


spindly individuals were produced by the 
3- and 4-week reapplication frequencies of 
These 


were confirmed by the diameter measure- 


100 and 250 ppm. observations 


ments. Diameter growth was also reduced 
considerably in the 250 ppm_ five-week 
treatment and the lanolin paste treatment, 
and “F” 
diameter further substantiate these obser- 
vations. 


values obtained for average stem 


The statistical analyses of survival and 
root weight data are interesting because 
Significant 
“FEF”? values were obtained in control vs. 


they follow a similar pattern. 


spray, levels and level x frequencies com- 
parisons and indicate important reductions 
in root development and survival were 
produced by the higher treatments levels. 
Dieback of seedlings during periods of rapid 
growth also appears to be influenced by 
treatment level. 

The influence of frequencies of reappli- 
cation fails to be significant when it is tested 


many against the interaction term, In a number 
TABLE 4. Analysis of variance of survival, dieback, height, diameter, and 
root development data. 
Source of variation 
Potal Control vs. Lanolin vs. Levels Frequencies LxF inter Error 
(21)! spray (1) spray (1) (2) (2) ction (4) (11) 
Survival 
SS 10461.3 386.0 599, 1 7310.9 615.8 991.6 537.81 
MS 386.0 399.1 3655.4 307.9 247.9 30.7 
F 7.6 11.8 14.7 1.2 4.9 
Dieback 
SS 797.6 111.8 16.7 273.0 63.4 33.4 299 4 
MS 111.8 16.7 136.5 31.7 8.4 are 
F 4.1 0.6 5.0* ee 3 
Av. Height 
Ss 2549.4 14.6 250.4 108.7 286.2 532.4 1357.2 
MS 14.6 250.4 54.3 143.1 133.1 123.4 
F 0.1 2.0 0.4 1.2 1.1 
Av. Stem 
Diameter 
Ss 433 059 00] .166 O88 069 050 
MS 059 00] 083 043 017 0045 
F or 0.2 4.8 25 3.8% 
Av. o.d, Wt. 
of Roots 
Ss 202.7 43.7 1.1 101.7 25.7 26.5 3.9 
MS 43.7 1.1 50.9 12.8 6.6 35 
F 124.8** 3.1 PP 1.9 18.9** 


1The numerals in parentheses indicate degrees of freedom. 
*Significant at 5 percent level; **Significant at 1 percent level. 
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of cases this is due to the interaction term 
being highly significant. Examination of 
the data often shows obvious reasons. Stem 
diameter and ovendry weight of the roots, 
for example, are little influenced by the 
frequency of reapplication at the 250 and 
10 ppm levels. At the high treatment level 
stem diameter and root growth are much 
reduced regardless of the frequency of ap- 
plication whereas the low level fails to in- 
fluence growth significantly even when ap- 
plied every three weeks. In contrast to this, 
at the 100 ppm level of gibberellic acid, 
stem and root growth appear to be very 
sensitive to frequencies of reapplication. 
This situation is reflected in the analysis of 
variance by a highly significant interaction 
term. 

The 250 ppm gibberellic acid treatments 
at all frequencies of reapplication were 
characterized by greatly reduced survival, 
a high percentage of trees showing evi- 
dence of the dying back of the crowing 
point, reduced diameter growth and little 
or no height growth response. Root devel- 
opment was greatly reduced. Leaves were 
commonly long and narrow, much reduced 
in total area, light yellow in color and 
many exhibited marginal discoloration and 
small scattered necrotic areas. Table 2 
presents data which indicates the maegni- 
tude of the treatment effects and Figure 
6 illustrates the typical reduction in root 
development that was observed in treated 
individuals at the higher treatment levels 
and frequencies of reapplication. 

Trees in the 100 ppm three-week and 
the 100 ppm four-week treatments, al- 
though having much improved survival 
over the 250 ppm treatments, exhibited 
similar symptoms of overstimulation. The 
100 ppm five-week treatment, although 
causing slightly reduced survival, greater 
dieback, and a somewhat reduced root sys- 
tem as compared with control, was the only 
treatment in which the treated trees exceed- 
ed the height growth and equaled the diam- 
eter growth of the control trees. Leaf 
size, shape and color, although slightly 
modified, was very similar to that of the 








control trees, 


The trees in the 10 ppm treatment, al- 
though showing a slight reduction in root 
development failed to differ from the con- 
trol trees in leaf size, color, and growth 
characteristics. The trees treated with 1 
percent gibberellic acid in lanolin paste 
survived well but failed to outgrow the 
control trees in height or diameter growth 
and exhibited a sizable decrease in root de- 
velopment. 


Fiber measurements fail to substantiate 


TABLE 5. Fiber lengths of seed- 
lings treated with gibberellic acid. 


Fibe: length in millimeters! 


Treatment Replication 1 Replication 2 Average 
Control 42 «41 36 (39 394 
Lanolin 44 «(O41 41 38 410 

10 ppm 

3 weeks .42 «41 40 36 388 
100 ppm 

3 weeks a sae oe .346 
250 ppm 

3 weeks .36 .35 ae sa aaa 
100 ppm 

5 weeks 36 =6.40 ae. «ot .378 
Control 

(re- 

check) AP 37 ae at .390 

IStandard deviation of the mean of one treatment: 


.007 mm; least significant difference: .020 mm. 
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original hopes that height growth increases 
would be accompanied by increases in fiber 
length (Table 5), Close scrutiny of the 
data reveals that only the trees from the 
three-week treatments (both the 250 and 
100 ppm) were significantly different from 
the control, and in these cases the fiber 
lengths were shorter. The trees in the 100 
ppm five-week treatment, although taller 
than the control trees, did not differ sig- 
nificantly in fiber length. 


Discussion 


The late summer slowdown in growth of 
the trees in the 250 ppm treatments, and 
in the three- and four-week treatments with 
100 ppm apparently was related to the de- 
creased root development and the decreased 
total leaf area of the treated individuals. 
Dieback of the growing tip, greatest during 
periods of rapid growth and under condi- 
tions of maximum moisture stress, was 
most serious at the higher treatment levels. 
This apparently was related to decreased 
root development and the presence of an 
increased amount of succulent, rapidly 
growing stem tissue. 

The lanolin paste treatment did not 
produce the response that might be ex- 
pected based on the preliminary trial. The 
reduced growth of the trees in this treat- 
ment in part seems to be related to reduced 
root development. The height growth data 
for the trees treated with lanolin paste is 
difficult to interpret in relation to the other 
treatments because it is the result of re- 
duced growth response, high survival, and 
a moderate amount of dieback; while in 
the case of the 250 ppm treatments, for 
example, the average height information 
listed is the result of excessive response, 
poor survival, and a fairly large amount of 
dieback. Survival was important in_ this 
latter instance because it was necessary to 
omit a rapidly growing tree from the 
measurements when it died back com- 
pletely, and to base the average figure pre- 
sented in the table on the remaining trees. 

Data fail to show any increase in fiber 
length due to the treatments employed. 
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This seems to be in conflict with the re- 
sults of Swan (1958) who obtained fiber- 
length increases of approximately 19 per- 
cent when working with hybrid poplars 
treated with gibberellic acid. His fiber- 
length increases, however, were associated 
with height growth response of approxi- 
mately 65 percent. 

Aspen is apparently quite sensitive to 
dosages and frequencies of application of 
gibberellic acid, Evidently a fairly fine bal- 
ance exists between stimulation and over- 
stimulation. With additional work on 
treatment levels and frequencies of appli- 
cation it seems entirely possible that greater 
height growth response could be obtained 
especially if fertility levels and moisture 
levels are maintained at nearly optimum 
conditions. Examination of fiber lengths 
under conditions of greater growth re- 
sponse would doubtless give a better in- 
sight into this aspect of the problem. 


Summary 


A preliminary greenhouse trial and subse- 
quent replicated field trials were carried 
out to study the effects of gibberellic acid 
on the growth and on fiber length of quak- 
ing aspen. 

1. The preliminary trial, conducted 
for an 8-week period, indicated that gib- 
berellic acid, applied as a spray or as a 
lanolin paste, affected growth under green- 
house conditions. Response was rapid and 
overstimulation resulted when spray treat- 
ments at concentrations of 500 ppm and 
100 ppm were applied every two weeks. 

2. Replicated field trials revealed that 
early height growth response was not as 
pronounced in the field as in the greenhouse. 
Height measurements after 8 weeks in the 
field indicated that seedlings treated with 
the higher levels of gibberellic acid (250 
ppm and 100 ppm) had the greatest aver- 
age seedling height. Later measurements 
revealed that late summer growth of the 
control trees was quite rapid and that there 
was little difference in the total height be- 
tween control trees and the treated trees 
by the end of the growing season. The 


)() 
r- 
its 


treated seedlings in the 100 ppm 5-week 
treatment were the only seedlings which 
exceeded the control trees in average 
height and average diameter. The diam- 
eter growth differences for this treatment 
were not statistically significant. 

3. Field survival after 6 weeks and at 
the end of the growing season was poor- 
est in the 250 ppm treatments and great- 
est in the control, the | percent gibberellic 
acid in lanolin, and in the 10 ppm level. 

4. Most dieback of the growing tip 
was observed in the higher treatment 
levels and the least in the control and the 
10 ppm treatments. 

5. Root development was best on the 
control trees and those in the 100 ppm 5- 
week treatment. Root development was 
greatly impaired in all treatments with 250 
ppm, the 100 ppm 3-week, the 100 ppm 
4-week treatments, and the gibberellic 
acid-lanolin paste treatment. 

6. Leaf discoloration at the margins, 
vellowing, and the presence of necrotic 
spots were noted at the 250 ppm and the 
100 ppm concentrations. Leaves of trees 
treated with gibberellic acid at higher levels 
were longer and narrower and had a total 
surface area much lower than the control 
and the 10 ppm concentrations. The leaf 
area differences along with the root de- 
velopment of some of the heavier treat- 
ments apparently accounted for the growth 
advantage of the control plants late in the 
summer. 

7. Fiber length determinations were 
made on a limited number of seedlings, 
four sets of measurements being taken for 
each treatment. Seedlings from the lano- 
lin, and the 10 ppm treatments had the 
longest fibers but did not differ significant- 
ly from the controls. The 100 ppm and 
the 250 ppm treatments when applied 
every three weeks, produced seedlings with 


significantly shorter fibers than the con- 
trol. Trees from the 100 ppm 5-week 
treatment had fibers that were slightly 
shorter but not significantly different from 
the controls. 
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Properties of Some Random and Systematic 
Point Sampling Estimators 


THe BrrrercticH method of sampling 
forest stands has won widespread ac- 
ceptance since it was first introduced 
(1948). The original conception of esti- 
mating stand basal area by counting trees 
was extended by Grosenbaugh (1952) to 
provide for estimates of volume, number of 
trees, or other stand variables, based on 
measurements of the sampled trees. He 
showed (1958) that such sampling was 
done with probability proportional to size 
and was capable of giving unbiased esti- 
mates of means and variances. Much work 
remains to be done on the variance and 
distribution of these estimates so that mean- 
ingful confidence statements can be made, 
the effect of changes in the sampling angle 
studied, and the relative effectiveness of 
this method and traditional plot sampling 
compared for particular forests. “This paper 
is presented as a necessary preliminary to 
the further consideration of these impor- 
tant questions. It will deal with the un- 
biasedness and the variance in the class of 
Bitterlich-type samples called point samples. 
Both random and systematic sampling plans 





FicGurE 1. 
rule. 


Geometry of the point sampling 
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A modified version of 
point sampling involving a subsample de- 
scribed by Bell and Alexander (1957) will 


also be cc ynside red ° 


will be examined. 


Sampling Procedure 


In the Bitterlich method of drawing a 
sample the primary sampling unit is a clus- 
ter of trees associated with a locus of origin. 
The locus of origin is a point in the case of 
point sampling, a line in the case of line 
sampling, an alternative but related pro- 
cedure proposed by Strand (1957). A tree 
is included within the cluster in point sam- 
pling if the distance, $, from the tree to the 
sampling point, divided by the radius of the 
tree, 1/2 D, is less than some constant 
value, A. The condition for including a 
tree in a particular sampling unit then is 


S/5D <K 


which may also be written as 

cosec 1/2 A < Kk 
where 4 is the angle formed by lines of 
sight from the sampling point to the oppo- 
site edges of the tree. Figure 1 pictures the 


situation in the plane of the breast-height 


The authors are, respectively, Lecturer in 
the School of Forestry, University of Califor- 
nia, Berkeley, California; and, formerly Statis- 
tician, School of Forestry, University of Cali- 
fornia, Berkeley. Appreciation is expressed to 
Elizabeth L. Scott, R. Zasepa, and W. G., 
O’Regan for reading versions of the manuscript 
and making helpful suggestions. Manuscript 
received Jan. 3, 1960. 
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of the tree, 4.5 feet above the ground. 

This angular relationship is the key to 
the simplicity of Bitterlich’s proposal for 
distinguishing the trees associated with each 
sampling point. If a fixed angle @ is so 
chosen that @ is twice the angle whose 
cosecant is kK, then a tree 1s selected for 
inclusion in the sampling unit if the angle 
formed by lines of sight to the edges of the 
tree is larger than the reference angle a. 
In this way, neither the distance from the 
sampling point to the tree nor the diameter 
of the tree need be measured to determine 
its inclusion. It will be seen that this type 
of sample will include more large trees and 
fewer small trees in the sample than their 
proportionate representation in the forest. 
In Figure 2, trees 1, 3, and 5 are included 
at point P, while trees 2 and 4 are not 
included. In practice, reference angles are 
small, of the order of one to four degrees. 
Typically a small number of trees, up to 
perhaps 30 in number, is associated with 
a sampling point. 

A number of convenient instruments 
have been devised by Bitterlich and others 
for projecting the reference angle outwards 
from the sampling point. Some of these 
compensate for the effect of sloping ground. 
Others such as the simple hand-held glass 
prism, may be adapted for such special 
circumstances. A good summary of the 
instrumentation of the method is given by 


“7 4 e 3 
eo Cy Sample point 
4 


Ficure 2. Association of trees with a sample 
point. 





Finch (1957). Adjustments to take ac- 
count of sloping ground and the elliptical 
shapes of tree cross sections are dealt with 
by Grosenbaugh (1958). 


Assumptions and Notation 


The plan in this paper is to find the ex- 
pectations and variances of some of the 
estimators suggested for use with point 
sampling. ‘The analysis which follows as- 
sumes a level forest of known extent, whose 
land area is denoted by L. There are N 
trees distributed over the forest in a fixed 
manner and so located that no tree center 
is nearer to the edge of the forest than the 
product of its radius, r, and a constant, K. 
The basal area (cross-sectional area 4.5 
feet above the ground) of tree 7 is sym- 
bolized by 4:. ‘The trees are all assumed 
to stand vertically and to be circular in 
cross section. There are no obstructions to 
visibility owing to distance or intervening 
vegetation. All measurements are assumed 
to be made without error. The size and 
arrangement of trees is perfectly definite 
and fixed. The only random elements are 
the locations of the sampling points making 
up a particular sample. 

The notation which is to be used is sum- 
marized here for the convenience of the 
reader: 

A =angle with sampling point as 
vertex, tangents to the tree as 


sides. 


«= twice the angle whose cosecant 
is K. (cosec a/2 = K) 


K =a numerical constant 


L = total land area in the forest 


f™ 
., 
yy 
to 

| 


= basal area factor 

ri — radius of the 7th tree 

6; = basal area of the ith tree 
zi = volume of the ith tree 


gi = xi/bi = ratio of volume to basal 
area of the th tree 


N = total number of trees in the for- 
eS sO 
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G= = 4; = total basal area in for- 
é=si 
est 
N 

V= 2 x = total volume in the 
i= 1 
forest 


n=number of points in a sample 
(j=1,2, ... 5) 

m ==number of points in a subsample 
(j=1,2, .. « 98) 

p=n/m, reciprocal of subsampling 

fraction; mm, and p all integers 

Yij=1 if ith tree included at jth 


sample point, zero otherwise 
N 


Nis = yi = number of trees in- 
é==1 : 
cluded at the jth sample point 


fa » ; 
£ — L/K*n ' 2 N j = Bitter- 
j=1 


lich’s estimator for total basal 
area 
, N 
v=L/K°n ( = 2 gyi = Gro- 
j=1 i= 
senbaugh’s estimator for total 
volume 
Oni=area of overlap between the K- 
circles of the Ath and 7th trees 
O'ni =overlapping co-incident areas 
between the A- circles of Ath 
and ith trees 
g¢ estimator for total basal area in 
a systematic point sample 
nm number of points in one syste- 
matic subsample 
n 
C =— = the number of independ- 
n 
ent systematic subsamples form- 
ing one sample 
¢g estimator for total basal area in 
a systematic point sample with 
C random starts 


The Basal Area Estimate 
Bitterlich’s estimator for total basal area 
(1948) is formed by multiplying the aver- 
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age number of trees included at a sample 
point by a constant basal area factor asso- 
ciated with the angle gauge used. ‘To write 
an expression for this estimator we define 
a set of indicator random variables Yi; 
which take on the value 1 when tree 2 is 
included at sample point 7 and the value 
0 when it is not. At each sample point we 
consider every tree in the forest and ask 
the question: “Is this tree included at this 
sample point?” If it is included, we assign 
a value of 1 to the random variable asso- 
ciated with that tree and that sample point, 
that is, we count the tree. In a random 
sample of 2 points with the basal area factor 
denoted by L/K* the estimator for total 
basal area! may be written 
n N 
os ee BD 2s (1) 
" Kn j=1 t=1 
where yij is the value (either 0 or 1) taken 
on by Yi. 


Sampling Properties. We wish to investi- 
gate whether g is in fact an unbiased esti- 
mator of G and whether a proposed esti- 
mator for the variance of g is an unbiased 
estimator of the true variance of g. It will, 
therefore, be necessary to find the expecta- 
tion of g, the true variance of g, and the 
expectation of a suggested estimator for the 
variance of g, under the assumption that 
the points making up a sample are ran- 
domly and independently located. 

As a step in finding the required expecta- 
tions we look first at the expected value of 
the indicator variable Yij, the only random 
element in the sampling situation. Since this 
variable can take on only the values 0 or 1, 


its expectation is simply the probability that 


1The quantity to be estimated here is con- 


ceived to be total basal area in the forest, 
» 


G = = Jj. This is closely related to average 
=} 

basal area per unit area in the forest, G/L. 

In the latter framework the basal area may 

be expressed either in square feet per acre or 

square meters per hectare, the basal area factors 


being, respectively, 43,560/K* and 10,000/K?. 


f 
n 


IS 


it takes on the value 1, that is the probability 
that the ith tree will be included when the 
jth sample point is located at random in the 
forest. ‘The area within the forest which 
meets the inclusion condition 8; < Kri is 
the circle traced by a radius Kri units in 
length and centering on the center of the 
tree. This enlarged circle or ring which 
may be conceived to surround each tree in 
the forest has been termed the K-circle by 
Grosenbaugh (1955) and is pictured in 

Figure 3. In symbols we may write: 

E(Yy) = Pr (Yy = 1)(1) 

+ Pr (Yi = 0)(0) 
__ ring area of ith tree 


land area 


m(Kri)* K*b; 
_ L ~ 
using the relationship that mri*7 = di, the 


basal area of the ith tree. 
This expression may be used in finding 
the expectation of g: 


E a ee, 

iss ane E(> E Yu 

, K-n j=1 #1 / 
I 


” N 
Sta ge a BC ¥e). 
K°nj=1 i=1 
Substituting the expected values of the Yy 
and summing over all trees and over the 
” points: 


‘ 
2 > 
nk" 2 bj 
EB i=) \ 
E(g) =——:————-= 2 &=G. 
K°n L or 
Bitterlich’s g is thus an unbiased estimate 


of G, the total basal area in the forest, 
under conditions of random sampling. 


Variance of g: It is con venient to define 
Ny 


N; = = Yi. Then the variance of gy 


4 
using the theorem on the variance of a 
linear combination is 
2 n 
2 __ L \ jos 2 
i" = —— << é 
Kn \ j=1 N; 
n—l n 
he eS \ 
—~ _ —_ . 
j=1k=itt NjNxf 


where o NN; is the covariance between 
. J+ . 


(2) 






K-circle of 
jt) tree 













Overlapping 
area of hth 
and ith trees 






Land area 
of forest 


Figure 3. K-circles and overlapping area, 


the tree counts at the jth and the 4th points. 
Since the sampling is at random, all the 
variances of the Nj are equal and the co- 


variances are zero. Then 


= if -L ); ‘ «+ 
o* = — |} no~ = - N, 
g K°n Nj nKk* 
Recall that, by definition 
a 
‘“ z _~.4 
o~ . =on — oC .. 
Ni Syy = Yo 
is} 
v—1 N Gow 
+22 Zz Yai Yai ( 


“y, = E(Y’y) — [E(%u) 2? 


= Ey 
~~ 


under random sampling, because Yij being 


an indicator variable, E(Y*i;) = EC Yu). 


And 
ee = 
YrjiYij — E(YnjYi) a, 
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‘Thus 
\ r2} \ 742 
on S A 2 A 
\ kK 1 L* 
vV—1 N 
2 2 = EY anjYi 
| rT i : . r4y 
a S Kt bnbi 
,: h 1 +1 L* 


which, combining the second and fourth 
terms, reduces to 


N; ce” Gee 


ohn 2 > = E( Y) jYay). 
} h+1 


To find E(Ynj¥ij) notice that YrayYis = 
1 if and only if the jth point falls in the 
overlap of the A-circles corresponding to 
the Ath and ith trees. The cross-hatched 
area in the two-circle case illustrated in 
Figure 3 is this overlapping area, which 


will be denoted by Oni. It follows that 


Pee ale Oni 
E(YnjYij) = 


Therefore 


es Esta 8a] 
4 1 


t z > Oni. (4) 
Lax 
And 
je iL oi 2 
o* KA a ~ 4-04 + 
é 1 7 L I 
a Mant N _ 
2 ae ; aL k- 
z 2 Gu>- (i) — —-@) 
Do a= t=h41 j K-n : 
21 V l N 
r— 2 2 Ox (5 
K*% nai i=n4t ) 


The variance of g is seen to be propor- 
tional to the land area and to the total basal 
area in the forest and to be inversely pro- 
portional to #, the sample size, as might be 
expected. An increase in the constant K 
will tend to decrease the variance of g, but 
this decrease may be offset by the resulting 
increase in the sum of the overlapping areas. 
The spatial distribution of the trees in par- 
ticular forests will affect the nature of the 
change in variance accomplished by an in- 
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creased value of K, 

With appropriate modifications, it is pos- 
sible to rewrite the expressions for estimat- 
ing total basal area and the variance of 
total basal area so that the more general 
case of forests with trees near the edge 
may be included. The iormulas reflecting 
the adjusted probabilities occasioned — by 
trees whose K-circles overlap the edge of 
the forest are not cited here because of 
their added length and complexity. 

It is shown in the appendix that the 
standard form for the variance of an esti- 
mator 

(6) 
” ” = 
. pin Bne- ( 3 v'] 
j=1 j=1 


z K*n 


s 


n(n 1) 
is an unbiased estimate of o*). An exam- 
ple of the computation of sy” is worked out 
below in connection with the discussion of 
the volume estimate. 


The Volume Estimate 


When the volume and basal area are deter- 
mined for each tree included at every point 
in the sample, total volume may be esti- 
mated by the following expression (Gro- 
senbaugh 1952): 

«Kn ai sane’ ay 4) 


-— , the ratio of volume to basal 


i 


where gi 
area of the ith tree and yi is the particular 
observed value of the random variable Yj;. 
The same procedure used to determine the 
expectation of g may be used to show that 
v is an unbiased estimator of V, the total 
volume in the forest. Similarly, 


2 fin N > 
ot =| ; = > 0:-0,,~ (8) 
K*n \, —1 i=1 ‘i ‘a 
N—-1 N 
-+ z a a qnqisy y 
h=1 ixhA+l hj ij 
2 r2 N r4 yy N 2 
_ & ja wa kK eo. 
Kn UL j=1 i Vent : 
> N-!1 N 
i. S . 


ut 
of 


The unbiased estimate of this variance is 
<a 


K*n 


§ —— 


x 


(9) 


\ (92 ] ) 


Computations. A numerical example will 
show the computations needed to estimate 
volume and basal area and the variance of 
volume and basal area estimates under 
random sampling. The data given in 
Table 1 are taken from some point sam- 
pling trials conducted by Cox in mixed 
conifer stands on the Latour State Forest, 
California.” For each of 44 sampling points 
the number of trees included, Nj, and the 


sum of the volume to basal area ratios for 
N 


these trees, 2 Givi, is listed. Other terms 
i 1 


needed are computed from the data in 
Table 1 as follows: 


44 

SNP PPP + P+. + 1P 
10,259 

44 N 2 

== avi 2228? + 18967 + 
I 

+ 30277 = 358,471,162 


The value of Z/K, the basal area factor on 
a per acre basis, is given as 9.372. Using 
formulas presented above, the following 
estimates” are obtained. 


Basal area per acre: 


9.372 
fe = —— - 617 131.42 sq. ft. 
44 


Volume per acre: 


“Data kindly provided by Paul Cox, Cali- 
fornia State Division of Forestry. 


8The symbols gz and eg stand for estimates 
of basal area and volume on a per acre basis; 
? } 
d 


Sg ~ and s, ~ are the estimate 
a a 


variances of 


basal area per acre and volume per acre. 


TABLE 


l. 


Tree count and sum of 


ratios at 44 sampling points, Latour 
State Forest. 


Sample Pt 


) 


19 
17 
16 


19 


volume 7, number 1, 


Sum of ratios 
\ 
~ 


1 


2228 
1896 

991 
2737 
2899 
2560 
3469 
+746 
2048 
762 
1649 
2798 
4466 
3343 


3230 


ait 


3067 
3658 
1685 
4268 
3344 
609 
93 
958 
1789 
2144 
1021 
3710 
3134 
1299 
284 
1329 
OSS 
3000 
2646 
1835 
2661 
4258 
3225 
1818 
1192 
2486 
4860 
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Variance of the estimate of basal area: 


2 __ (9.372)* 44(10,259) — (617)* 
— 44-43 
= 74.60 


Variance of volume: 
(9.372)? 


ce ae 
44(358,471,162) — (112,832)? 
hh 43 ar 
= 3,209,246 


Volume estimate of Bell and Alexander. 
Bell and Alexander (1957) have proposed 
a modified procedure for estimating vol- 
umes in which the number of trees Nj, in- 
cluded at a sampling point is determined 
at all the sample points, but the sum of 
volume to basal area ratios is determined 
at only a fraction of the sample points. 
Their estimator for total volume may be 
written 


n 
v N; 
L ia m  N 
a — "$= a5 >> = 4ivu (10) 
Rs & x, ” 
DN; 


where x = the total number of sample 
points and m is the number of points at 
which sums of ratios are found. A numer- 
ical example showing the computation of v’ 
is given below in connection with the com- 


putation of the mean square error of v 


Us 


Sampling properties of v: The charac- 
teristics of v are investigated under the 
assumption of rancom sampling. Multiply- 
ing both numerator and denominator of 


(10) by L/K*m gives 











ot ads 
© ; I ¢ 
Kkr-m L - 
i ems aa N; 
Se K*n j=1 
— = N; 
Kn 
oa . - 
5S 
ao 


in which v and g are the estimates formed 


from the subsample and g’ is the estimate 
of basal area based on the total sample. Bell 
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and Alexander’s v is seen in this form to 
be a ratio estimate. Following Hansen, 
Hurwitz and Madow (1953, vol. II, 
Chap. 4) the expectation of v and its 
mean square error are: 


E(v) =VX 


cr 9 + 
| 1 ' Ov gt Orig i Og Og,q' | 
+o et a (1) 
\ J G VG Go G- J 
and 
o*, a Oo", 


M.S.E.(o’) = V2) 72 42% 428 
(v) ete te 


4 
L071 go" 200s 


+ 


2Gaa 
oe *. (12) 


VG VG CG 


The variances and covariances linking the 


large and small samples are as follows: 


. zz 
Og: = — Og 
p 
_ 
Cg.g¢ — O~ gq 
l 
Ov gt —_ — Orig 


where the integer p = n/m. These rela- 
tionships are derived in the Appendix. The 
expression for the expectation of vy (11) 
indicates that it is a biased estimator of 
total volume. Making use of the variance 
and covariance relationships linking the 
main and the subsamples, this bias may be 
written as 

E(v'—V = v(@*) (Se — N13) 

p G* J Gj 

or in terms of the correlation coefficient be- 


tween wv and g, Ppr,g: 


E v—V) — V (A 
p 


——* 


9 
(= pr oF Fy \ 
= | 


I\G? VG | 


=v (A) tee — erate) (14) 
\ p/IG\G Vv 


This form of the expression for the bias 
indicates that if 
7 o,/G 
Pvig o./V 
the bias of v is sure to be small. 

Because of the bias of the ratio estimate 
v, we work with the mean square error 
of v rather than with its variance. The 


to 
en, 
II, 
its 


+) 


las 


expression for the mean square error of v 
(12) may be written: 





Ps 
MS.E.(v)=V?4o2 po? 0% 
iv’ Gg eG 
2 a4 
p? Trg tl 
re Ve eI 
and combining terms further 
: , oe | Fv 2(p—1) Ov,9 
M.S.E.(v ) = V* LV = p VG 
gi! oy" \ 
» GJ 
In relative terms, we may write: 
(15) 
MS.E.(v ) oy p—1 [org o6° 
— vw Pe 4 ee “3 


It follows from this expression that there 
will be no advantage in using vy rather 
than v whenever, in terms of the correla- 
tion coefficient, 
2p, a calG 

‘a 

In order to make an informed choice be- 
tween a single sampling plan to estimate 
volume using v and a double sampling 
plan using v, a knowledge of the relative 
costs of counting trees and measuring their 
volumes and basal areas is needed, as well 
as some information on the variances and 
covariances in the particular forest under 
consideration. 

In order to estimate the relative error 
of v, we follow the suggestion of Hansen, 
Hurwitz and Madow (1953, vol. I, p. 
176). The simplest estimator of this rela- 
tive error follows directly upon Equation 
CES}: 

S(v)_ sv_ >? “1 2Qsr.9 att . (i 

u u Pp \ vg g 
The necessary estimates of variances, co- 
variances, volume and basal area are based 
on the subsample. The large sample is used 
only for the estimate of volume, v. Table 
2 gives the numbers of trees included at 
each of 44 additional sample points on the 


TABLE 2. 


tional sampling points, Latour 


Forest. 


Tree count at 44 addi- 


No. trees 





Sample Pt. included 
7 Nj; 
45 19 
46 19 
47 12 
48 22 
49 14 
50 12 
51 11 
52 10 
53 12 
54 25 
55 10 
56 9 
57 9 
58 6 
59 9 
60 16 
61 7 
62 9 
63 9 
64 7 
65 13 
66 6 
67 14 
68 10 
69 19 

0 22 
71 16 
-5 9 
73 9 
74 Ss 
- 7 
76 24 
77 4 
78 10 
=9 2 
80 24 
81 13 
82 14 
83 10 
84 24 
85 25 
86 22 
87 2 
88 20 

Total 617 


State 
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Latour State Forest. The a points given 
in Table 1 are considered to be the sub- 
sample, the 88 points listed in ‘Tables 1 and 
2 to be the large sample in a numerical 
example illustrating the computation of y 
and of s*(v)/v*. In terms of sample size 
we have SS. 44. and p= ‘a 
The estimate of volume by Equation (10) 
is 
9.372 1234 
Va . °112,832 = 24,033 bd. ft. 
88 617 

We note that in this particular case v’ U; 
because it happened that in the example 
chosen the number of trees included at 
88 points was exactly twice the number 
included at 44. 

To estimate the mean square error of 
v we make use of the quantities computed 
for the numerical illustration given above. 
In addition, we need the estimated covari- 
ance of v and g based on the data given 
in Table 1 and the total number of trees 
included at sample points 45-88 of Table 2. 
An unbiased estimate of the covariance of 


v and g is: 


L? n XZ X qiyy (Nj) 

Srig . 
K+» m(m—1) 
Go qi) a 
yi 1) 
(17) 

Referring to Table 1, note that 
m if N 
> Z qil (N;) 


= (2228)(9) + (1896)(9)+ 


+ (3027 )(17) 
= £575,216. 


We have then 


(9.372)- 
44 
(44)(1,875,118) — (112,832)(617) 


44-43 
13,596.38 


60 / Forest Science 


and by (16) 


P(v') 3,209,227.9 [2-1 
u™ (24:033.22)* {| 2 
| (2)( 13,596.38) 74.585 
(24,033.22) (131.42) (131.42)? . 
= (),0055562 5 1 2(0.0043048 ) 


(0.0043185) | 00034106. 


In percentage terms, in the illustration given 
(7) 


© 100% = V0.0034106 « L00% 


> 5.8% 
based on 88 points and a subsample of 44, 


compared with 
s(v) 


Jt 


62 . 1! O% 


© 100% V0.0055 


5 ¢ 
ua Pe 


based on just the 44 points without the use 
of double sampling. 

The proposal which is offered by Bell 
and Alexander (1957, p. 17) for esti- 
mating the relative error s(v )/v can only 
be thought of as an approximation, It 
consists of two terms, one of which is 
closely related to s(g¢)/z, the other of 
which expresses the relative error of tree 
ratios with respect to the average of tree 
ratios. There is a confusion here, since in 
point sampling, the measurement we are 
concerned with attaches to points (i.e. the 
sum of the tree ratios at a point) rather 
than to trees. The formula for the relative 
error prese nted above (16) is offered as 
being conceptually sounder than the ex- 
pression of Bell and Alexander. The sub- 
sampling fraction and the effect of the 
covariance between v and 
1 


are explicitly 


recogeniZ 


Systematic Sampling 


To this point the discussion has assumed 
that random sampling is followed. We 
now consider the effect on the foregoing 
analysis of drawing the points making up 
a sample in a systematic way. This exten- 
sion of the argument is appropriate because 
of the widespread application of systematic 


+S ) 


06, 


ven 


44 


, 


use 


sel] 
sti- 
nly 


ferest sampling in practice. 

Specificially, we wish to investigate the 
estimation of total basal area under condi- 
tions of systematic point sampling with a 
single random start. The origin of a syste- 
matic sample of 7 points may be any point 
within the forest. The other »—1 points 
of the sample are related to the initial 
point in the following manner: the second 
point is p units distant from the first, and 
its orientation with respect to the first point 
is in a direction which is drawn at random. 
This direction defines the orientation of a 
rectilinear grid coextensive with the forest 
and with points p units apart in one direc- 
tion and qg units apart in the direction nor- 
mal to that first drawn. Points are located 
at the intersections of this grid outwards to 
the forest edges until the required number 
n has been located. It is assumed that fp 
and g are chosen large enough that a given 
tree will not be selected more than once 
in any given systematic sample. A syste- 
matic sample of » points drawn in this 
fashion may be thought of as a random 
sample of size one out of an infinite num- 
ber of possible systematic samples. 


As in the random sampling case, we 
construct the random variable Yyi which 
is capable of taking on only two values: 
vy 1 if the sth tree is included in the fth 
systematic sample and yy; 0 otherwise, 

The probability that tree 7 will be in- 
cluded in any systematic sample will be » 
times the ratio of the area its K-circle to 
the total land area of the forest. That is: 


nk*! 


PCY; 1) L 


For each systematic sample, an estimator 
for total basal area is: 


L tl «= 
K* 
which corresponds with g used previously 
in the random sampling case (1). The 
estimator for total basal area under sys- 
tematic sampling is unbiased: 
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The variance of gr may be found using the 
identity 

ogi = E(¢*) [E( gr) ]*. 
The first term may be expanded as fol- 
lows: 


i : 
E(2*;) : E( = yn)" 
K*n 
2 = Eyy?+2 
K*n 
\ l V 
yy 


- — E yriyyn 
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The indicator variable ysiyyn takes on the 
value 1 when trees 7 and /: are both in- 
cluded in systematic sample f, otherwise 0. 
It is not necessary that both be included at 
the same sampling point, only that each 
one be taken in at one of the » sampling 
points making up the systematic sample. 
The probability that two particular trees 
will be included in a systematic sample 
may be visualized in the following man- 
ner, Consider a rectangular grid whose 
sides equal the spacing of sample points in 
the system and let the grid be superim- 
posed on the forest in such a way that the 
K-circle of tree A lies completely within 
one rectangle, as pictured in rectangle 3 of 
Figure 4. Then consider the A-circle of 
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FiGureE 4. Systematic sampling grid in relation 


to two represent itive trees. 
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tree i which lies partly in rectangle 1 and 
partly in rectangle 2. Conceptually pick 
up rectangles 1 and 2 and make them co- 
incide with rectangle 3 in the manner rep- 
resented in Figure 5. The area common 
to the A-circle of tree A and the imposed 
portions of the A-circle of tree 7 is termed 
the “overlapping co-incident area” and de- 
noted by O'n. The sum of this overlap- 
ping co-incident area in relation to the total 
land area of the forest is the probability 
that two trees will be selected at a point 
of the systematic sample. 

Oni 

L 


Using this expectation in the equation 


E(ysiym) = Prlyriym = 1] = n- 


above we have 
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‘Thus, we may write an expression for 
the true variance of the estimate formed 
under systematic sampling. We know of 
no unbiased estimator for this variance, but 
it may be useful to compare it term by 
term with the corresponding expression for 
random point sampling. We may rewrite 
the variance of a random point sample esti- 
mate (5): 


2 | 1 
o-(¢) random = —— hy 7 
AH 1 n 
( 2 bi) 
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Note that the first term is common to both 
expressions and that the negative term in 


the expression for systematic sampling 
variance is ” times greater than the corre- 
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The concept of overlapping co- 
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sponding term in the random sampling 
variance. The overlapping area term, how- 
ever, will always be larger for systematic 
sampling because co-incident overlaps at 
different points are considered as well as 
direct overlaps at single points. We can 
make no general statement about the rela- 
tive magnitude of the two variances be- 
cause we do not know the effect of the 
spatial arrangement of trees or the size of 
the overlapping area terms. A disturbance 
may be introduced into an attempt to esti- 
mate the variance of the systematic sam- 
pling estimator if there should be an un- 
foreseen periodic distribution of variability 
in phase with the systematic sampling pat- 
tern. 


CO=- 


W- 
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‘To guard against such a possibility and 
to help assure an estimator with lower 
variance under systematic sampling, con- 
sider an alternative plan of systematic sam- 
pling with C random starts, C being a sim- 
ple integer of the order of 3 or 4, Each of 
the systematic samples now contains » 
points, the total number of points in the 
sample being the same as with a single start, 
n = Cn. The spacing of points is Cp and 
of lines Cg, compared with the p and qg of 
the single random start case. The orienta- 
tion of the first systematic sample governs 
that of the remaining samples. The esti- 
mator for total basal area, g, is now the 
mean of the estimates provided by the C 
samples: 


. «2 , Jb = « 
ae See ee ee 
i G 3 1 C K*°n 1 f=1 


It differs from the g of random sampling 
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only in the manner of selection of the 


points. We may show that it is unbiased: 
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We note that this expression is valid also 
for the limiting cases of random sampling, 
setting C = m and w 1, and for sys- 
tematic sampling with a single start, that 
is, C = land » =n. The unbiased esti- 
mator for this variance is 


- | oe ise 
Mea= C2 C1 
| Sy C <¢ - 
K*C7n? (C-1 2 ts 
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Since C is a small number, this estimator 
tends to have high variance. We may there- 


fore alternatively compute s*(g), random, 


as an approximation to the desired variance 
and argue that the taking of additional sys- 
tematic samples has rendered it unlikely 
that we will underestimate the variance of 
systematic sampling when we treat it as 
though it were random. 

‘The same approach may be applied with 
necessary changes in the expressions to the 
estimation of total volume, using systematic 
sampling. 


Summary 


We have shown that estimators g, g/, and 
g which have been proposed for basal area 
in random and systematic point sampling 
are unbiased; and similarly that the corre- 
sponding estimates for volume are unbiased. 
We have seen that the estimate for volume, 
', formed by sampling for basal area at 
points and sub-sampling for volume at a 
smaller number of points out of » is 
biased, but with the bias small when speci- 
fied relationships hold for basal area and 
volume parameters. 

Expressions for the true variances of g 
and £ have been derived and it has been 
shown that unbiased estimates of these 
variances sg” and 5s . exist, and similarly 
with the corresponding variances for vol- 
ume. The mean square error of vy has 
been derived and an estimator for this mean 
square error has been suggested. 

The conditions under which sub-sam- 
pling for volume is likely to be advanta- 
geous have been formulated. A procedure 
for doing systematic point sampling based 
on several random starts has been sug- 
gested whereby some of the advantages of 
both systematic and random sampling may 
be retained. 

Supplementary Discussion 

Appendix 1. Unbiased estimator for the 
true variance of total basal area under 
random point sampling: 

In this appendix we demonstrate that 
Grosenbaugh’s estimator for variance of to- 
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tal basal area is an unbiased estimator of the 
true variance. We evaluate the expectation 
of sy”, formula (6), and compare it with 
oj, formula (5). The expectation we are 
seeking 1s: 


ae ay 
Es*, : 
kK ay 
: i ‘i 
: 2 NF Pe s 
E ‘ 1 \ 1 ( 18) 
(9 ] ) 
First, we find the expectation of Nj and of 
N 
N;*, using the definition Nj = = Yi: 
I 
= 5 F i leg 
EN > EY a eS G 
i l Z: i 1 i, 


as shown above in the derivation of for- 
mula (1); and 


EN? =E (3 Yu) 


N N 1 N 
= ( =>y%4+22 2 SY; ru] 
h=! h+1 

as seen in the development of formula (4+) 
above. To find the expectation of (27_, 
Nj)" we square the sum and take the ex- 
pectation of the cross-products, remember- 
ing that Nj and Nx are independent: 


E (2) 


rr aa ‘ 12 es 


hen, summing and substituting into equa- 
tion (18) we have: 


siang hy 
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which is equivalent to o,°, the formula for 

the true variance (5). 

A ppendix Zz; 

expression for mean square error of 7 
Following Hansen, Hurwitz and Madow 

(1953), vol. 1, p. 514) a good approxima. 


Evaluation of terms in the 
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tion for the mean square error of v is 
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Expressions for 0“, and o*, are given above 
(5), (8). Similarly, the covariance of v 
and g may be written 
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The variances and covariances involving Z> 
the estimate of basal area based on the large 
sample of size 7, may be expressed in terms 
of those for the small sample of size m, 


since 7 mp, p being an integer: 


; L I = 2 


= TOs 


é . ro ee 
, K- mm) jo 1=1° 
m N 
.£ 1.22 
; >> _ = Vi 
k- mp i=1 = 1 
2m Vv 
+  Zyyt... 
j m+1 i 1 
mp N 
_ Bs 
7 - — Vij. 
j (p—1)(m+1) i 1 


Each term of this sum is of the form g/p. 
The variance of this sum of p independent 
terms is 


o* ( g ) = — ; : 

In finding the covariances, use is made of 
the theorem that cov a($+c) = cov (ab) 
— £OV (ac): 
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+ other covariances of this form. 
All the covariances but the first are zero, 
since the others are based on a different 
random set of points and the associated 
measurements are independent of one an- 
other. It follows that 


and, by the same argument, it can be shown 
that 


Orgs = - Org. 
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Organic Layers in Minnesota Aspen Stands 


And their Role in Soil Improvement 


QUAKING ASPEN (Populus tremuloides 
Michx. ) a post-fire pioneer species, is the 
most widely distributed broadleaf tree spe- 
cies found in North America. In the three 
Lake States, it occupies about 19.8 million 
acres of land, or approximately 39 percent 
of the total forest area (Chase 1947). It 
is especially prevalent on burned-over up- 
land forest soils formerly occupied by pine 
or northern hardwoods. Although not as 
valuable as pine or the northern hardwoods, 
it is important economically and is now the 
leading pulpwood species in this region, 
supplying about one-third of the total an- 
nual cut (Horn 1950). Formerly con- 
sidered a weed species, aspen now occupies 
a significant position in the production of 
fiber, lumber, and specialty products. 

Less appreciated is the important role of 
the species as a soil improver, and as a 
pioneer crop preceding the establishment 
of later more valuable stages of more tol- 
erant sub-climax and climax types, notably 
balsam fir (Abies balsamea (L.) Mill.), 
and northern hardwoods, such as sugar 
maple (Acer saccharum Marsh.). This 
paper will document the role of aspen as a 
soil improver after forest fires, particularly 
its role in redistributing nitrogen to the 
surface layers, whose nutrient status is 
vitally important in ecological succession. 

Daubenmire (1953) found aspen foliage 
to have the highest nutrient analysis of any 
of 12 western tree species studied, leading 
all others in average content of total nitro- 
gen and potash, and being second in phos- 
phorous and calcium content. 

Alway and Rost (1928) came to the 
conclusion that forest fires were detrimen- 
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tal in upsetting the overall nutrient bal- 
ance, since severe burns can destroy all or 
much of the organic layer, ranging from 
7 to 26 tons per acre of oven-dry material, 
and representing a loss of nitrogen of 450 
to 1500 pounds per acre. The mineral soil 
itself was not appreciably affected by burn- 
ing, being only .02 percent lower in nitro- 
gen content than unburned areas. Rowe 
(1953) states that intense fires may con- 
sume the entire humus layer and reduce 
vigor of aspen suckers. 

In connection with a study of soil prop- 
erties and site index of aspen in northern 
Minnesota and Wisconsin (Stoeckeler 
1960), there was opportunity to make ob- 
servations on the weight of organic layers 
of quaking aspen stands and to assess the 
role of aspen as an improver of the surface 
layers, especially with reference to nitro- 
gen, Nitrogen appears to be the nutritional 
key in advancing plant succession toward 


more valuable climax types in the Lake 
States. 
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Study Methods 


Separate fractions of the organic layers in 
29 quaking aspen stands in nothern Minne- 
sota were sampled by laying out from 3 to 

sampling areas 6 by 6 inches in size. 
‘These were sharply defined by a steel fram: 
+ inches high. Samples of F and H layers 
were kept separate. On 14 plots the L 
layer was also sampled. The L layer was 
defined as recent undecomposed leaf fall, 
I as a partially decomposed litter where th 
ives were as yet not fully decomposed, 
and H as the layer of well decomposed lit- 
ter not recognizable as to origin (Hoover 
nd Lunt 1952), 


Sticks,. stones and large roots were re- 
moved and the samples placed in numbered 
ags, later spread out at room temperature 
to air drv thoroughly. Once air dried, the 
samples were weighed, a small portion 
oven-dried, and the air-dry weight recal- 
culated on an oven-dry weight basis. These 
data were evaluated by correlation analysis 
(Snedecor 1937). 

The mineral portion of 21 of the profiles 
was sampled by horizons, five of which 
were sampled in detail for total nitrogen 


analysis. “Total nitrogen was determined 
by the Kjeldahl method (A.O.A.C. 
1950); cation-exchange capacity and re- 
placeable bases by the method of Chapman 
and Kelly (1930); available phosphorous 
was determined by the Truog method 
(1930); and, pH was determined electro- 
metrically with the Beckman pH meter. 


Site index was determined by sampling 


} 


from 5 to 7 representative trees of th 
dominant and codominant class for total 
height and age, and predicting all sites to 
age 50, using curves developed by Kitt- 
ridge and Gevorkiantz (1929). Volum 
data were calculated from tables prepared 


t 


by Gevorkiantz and Olsen (1955). 
Weight of Organic Layers 


The combined oven-dry weight per acre of 
all organic layers varied from 19 to 56 
thousand pounds (Fig. 1). Based on th 
average line, the values for 10-year periods 
from the 10th to the 40th year are ap- 

2 2 thousand 
Projecting the line 


back a few years to zero age 


' ee . 
proximately 35, 43, and 5 


pounds, respectively. 
it would in- 
tersect at about 19 thousand pounds. ‘This 


indicates that on the average these plots 


had probably almost 10 tons of oven-dry 
material present at time of establishment 


of the aspen stand. This weight of organi 
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Figure 2. Oven-dry weight per acre f com- 
bined F + H layers of 29 Minnesota plots 
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Ficure 3. Oven-dry weight per acre of litter 
layer of 14 Minnesota plots as 
of stand. This 
nual litter fall. 


rélated to age 


represents essentially the an- 


matter was hence left from the previous 
stand of pine, hardwood, or balsam fir- 
spruce after it had burned, 

On 29 plots where only the weight of 
F and H layers was determined, the values 
ranged from about 17 to 62 thousand 
pounds, The residual amount of these 
layers at zero age of the aspen was about 
23 thousand pounds (Fig. 2), or in close 
agreement with the more limited sample 
from 14 plots (Fig. 1). 

The weight of the litter only, for 14 
plots, shows an increase at least to age 35 
(Fig. 3). Based on the curve values, the 
plots approximately 33 years old had about 
3200 pounds per acre of litter, while in 10- 
year-old stands it was slightly over half 
that amount. There is some evidence that 
it tends to level off gradually, 

The weight of F layers observed in 29 
plots, ranging from 9 to 70 years, shows a 
substantial increase with age (Fig. 4). 

The H layer has a very similar trend 
(Fig. 5), but the amount is substantially 
heavier, constituting as much as 70 to 90 
percent of the entire combined weight of 
the Ao. Here again, there is evidence that 
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the residue from the previous stand was 
largely H layer and that it must have been 
about 17 thousand pounds per acre. 

There is apparently a rather high per- 
cent of the organic layers which becomes 
completely decomposed or is incorporated 
into the soil by the action of fungi and 
other soil micro-organisms. “The combined 
organic layer production curve (Fig. 1) 
shows an apparent gain of some 16.5 thou- 
sand pounds in weight between ages of 10 
and 30, or about 825 pounds per year. In 
this 20-year period the average annual lit- 
ter production was probably two to three 
times that amount, ie. 1614 pounds 
(Owen 1954) to about 2550 pounds 
(from Figure 3 rather than litter traps). 
The total calculated increment for the 20- 
year period is thus probably in the range of 
32.3 to 51.0 thousand pounds, but the 
actual gain was only 16.5 thousand pounds, 

The 51 to 68 percent which disappeared 
is thus a rough estimate of overturn rate. 
This rate for aspen is somewhat higher 
than that estimated for hardwood leaves by 
Maki (1951) who used metal containers 
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Figure 4. Oven-dry weight per acre of F 
layer of 29 Minnesota plots 
of stand, 
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en weights of leaves which were reweighed a 
as f «a50 
periodically. wz 
T= a2 ° 
1es Nitrogen Distribution in the Profile “a 40 
Tu 
ed Nitrogen analyses of each of 5 profiles —— 
P 
nd rr? 11 . t " on 30 
(Table 1) indicates that the F layers had 3 
ved substantially higher values than H_ layers = 3 
1) wee ally gne alues ae ayers. Ow 20 
Ihe Ai layers averaged (0.327 percent, ws 
10 while Az averaged only 0.038 percent. = 10 
( or = : 2 ~ 
I The calculated amounts of total nitro- 
n patel : ¢ saan ie < 
gen in terms of pounds per acre in these 5 
lit- Sails pana acca die tadeeiiie Ae 0 10 20 30 40 50 60 70 
ae profiles considering oven dry weight, depth, AGE OF STAND (YEARS) 
and bulk density were as follows: F = 
nds 188 lbs.; H = 363 lbs.; Ar = 473 lbs.; Ficure 5. Oven-dry sceight per acre of H 
a lavee nt ; »nta plote as slated tn ave 
nds Az = 294 lbs. The average oven-dry on aad fo? Minnesota plots as related to age 
. 6 . / la lad. 
m4 weight of the F and H layers was 11,042 — 
()- > : 
- ¢ and 34,032 pounds per acre respectively 
: . for these 5 plots. 
“4 Since forest fires, from the evidence at in northern Wisconsin and Michigan. 
uC 7 hand, often destroy all the F laver and These species have rather high nutrient re- 
— about half the H layer, the loss of nitro- quirements. 
aoe gen in stands of the type discussed here Quaking aspen gradually restores not 
we would be about 370 pounds per acre. If all only the nitrogen fraction in surface layers 
ou, of it were destroyed, the amount would be by its accumulation of humus, but also bases 
wars 551 ls, in the range indicated l such ; ash, calci and magnesi 
51 pounds, in the range indicated by Al- such as potash, calcium, and magnesium. 
way and Rost (1928). These minerals are brought by the roots 
The loss of a substantial amount of or- from greater depths, particularly the B and 
ganic matter and nitrogen in forest fires C horizons and through leaf fall they en- 
and reduction of water-retaining capacity rich the H layer which contains an in- 
1S considered a prime factor in setting back tensive network of fine feeder roots, These 
plant succession of the economically valu- re-absorb some of the nutrients to continue 
able climax species, such as balsam fir in the cycle. 
Minnesota, and sugar maple and balsam fir There are large differences in the quan- 
Pr 
) 
TABLE 1. Nitrogen distribution in five Minnesota soil profiles. 
1 Plot Site Per cent of total nitrogen in each horizon 
) No. index Soil t pe! F H Ay \o By Bs Bs Cy Cy 
\ I 58 Nebish L.f.s. 1.520 580 .210 .024 .018 .023 — 021.018 
4 62 Nebish 1.f.s. 1.650 1.200 025 .015 .025 .027 .028 .024 
8 66 Dora Lake f.s.]. 1.790 .840 .280 .047 .039 .031 .009 .043 .022 
11 66 Beltrami si.]. 1.710 1.220 .390 .055 .042 .056 .055 .029 .021 
14 66 Nebish loam 1.840 1.470 .430 .041 .059 .046 026 .024 
70 ) Average 1.702 1.062 .327 .038 .035 .036 .030 = .029 = .022 
Ratio to Cs 78 48 15 2 2 2 1 l 1 
) ; a i a a i 
of F tl 
a 1Abbreviations are as follows: l.f.s. = loamy fine sand; f.s.l. = fine sandy loam; si.l. = silt loam. Recent 
s, ake 


changes in soil series nomenclature are such that the two Nebish loamy fine sand profiles would probably be mapped 


as Kinghurst. 
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TABLE 2. Comparison of pH and of nutrients in several horizons of 21 Min- 


) “7 


nesota soil profiles, some with lime-rich parent materials and others with acid 


a acl naar ae 
pare WL *atleriats. 


Paren Total N 
na ] H« mn pH ( en 
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tities of nutrients in the various horizons. 
The ratios of the values for nutrients in 
H layer to the Ae layer is especially strik- 
ing (Table 2). In comporing lime-rich 
with lime-poor profiles (based on efferves- 
cence of the Ce horizon with dilute hydro- 
choloric acid) there was a considerable con- 
trast in calcium and magnesium content 
for the H and Ag layers and a very marked 
difference in the Be layers.’ Calcium was 
4.4 times more abundant, and magnesium 
5.4 times more abundant in the Bz horizon 
of the profiles where the parent material 
showed effervescence with the acid, The 


latter are in the gray-wooded classification 
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Avail. P= Repl. K—-Repl.Ca_—_—Rep!. Mg 
(ppm) (ppm) (ppm) (ppm) 
170.0 86.0 779.0 81.0 

“ia rat 39.6 Sek 
34.0 17.4 190.8 39.0 
147.0 94.0 521.0 64.0 
29.0 9.3 33.3 1.4 
32.0 9.8 43.9 72 
161.0 89.0 680.0 75. 
28.0 g.0 37.2 +.9 
31.0 14.5 134.8 26.9 

5.8 Pi. 18.3 $5.3 

1.0 1.0 1.0 1.0 

1.1 1.8 3.6 5.5 
reas was 72 a 62 respect V ’ ’ 

65 and 58 cubic feet per acre for lime-rich and lime- 

Mg refer to phosphorous, potash, calcium, and ma 
n i 

Summary 


In stands of quaking aspen in northern 
Minnesota the litter per acre varied from 
| 


about 320( 


) pounds of oven-dry material in 
33-year-old stands to slightly over half that 
amount in 10-year-old stands. 

The combined weight of the F and H 
layers varied from 17 thousand to 62 thou- 
sand pounds per acre and was related to 
age of stand. 

The H layer constituted as much as 70) 
to 90 percent of the weight of the entire 
organic layer. It is deduced that the forest 
fires which established the aspen stands did 
not destroy all of the organic layer of the 
previous stand but left a residue of around 
17 thousand to 20 thousand pounds of or- 
ganic material per acre, or about half the 
weight of the combined F and H layers 
found in the sample aspen stands at age 40. 

The net annual increment in weight of 
organic layers between ages 10 and 30 for 
29 plots was about 825 pounds per acre. It 
was calculated that from one-half to two- 
thirds of the organic matter produced was 


lost by incorporation into the soil or disap- 
peared due to action of organisms, and 
leaching. 

The calculated av erage amount of nitro- 
gen in the organic layers of 5 well estab- 
lished aspen stands was 188 pounds per 
acre in the F layer, and 363 pounds per 
acre in the H layer. It was estimated that 
repeat burns in stands of the type discussed, 
often consume about 370 pounds of nitro- 
gen per acre, thus causing a setback in the 
ecological succession of the economically 
valuable climax species having rather high 
nutrient requirements. 

Quaking aspen gradually restores the ni- 
trogen in the surface layers of burned-over 
sites by deposition of leaf litter as well as 
restoring bases such as potash, calcium and 
magnesium brought up by root action from 
lower levels. 

The H layers in aspen stands growing 
on soils with an abundance of lime car- 
bonate had a considerably higher content 
of available calcium and magnesium than 
those on the more acid soils. The former 
had a site index of 10 feet more or 16 per- 
cent greater, and a mean annual increment 
of 7 cubic feet or 12 percent more, per 
icre, than the more acid soils. 
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Equations for the Description of 
Past Growth in Even-Aged Stands 


Of Ponderosa Pine 


OvER THE YEARS many mathematical 
models have been proposed to describe bio- 
logical growth processes. Some of these 
have been developed by foresters; others, 
originated by workers in other fields, have 
been applied to forest trees. There is still 
no agreement on the one best expression 
for the description of growth phenomena. 

A group of Japanese investigators, study- 
ing growth and intraspecific competition 
in herbaceous plant communities, have 
proposed several equations to describe the 
relationship between mean plant weight 
and stand density at various ages. These 
equations give a satisfactory fit to data de- 
rived from several even-aged stands of 
ponderosa pine (Pinus ponderosa Laws.) 
in northern Arizona, and appear to offer 
promise as a means of analyzing some as- 
pects of growth in forest stands. 


Collection of Data 


Field work for this study was carried out 
in the fall of 1957 at two locations on the 
Fort Apache Indian Reservation in east 
central Arizona (1) near Bog Creek, about 
10 miles east of the town of McNary and 
(2) on Odart Mountain, some 7 miles 
west of Maverick. These two points are 
about 40 miles apart, but are closely simi- 
lar in elevation, topography, and soils. Both 
are in the heart of some of the best and 
most productive pine forest in Arizona, 
Sampling was limited to pure stands of 
even-aged ponderosa pine of pole size. In 
the dry climate of the Southwest, only in 
occasional years are favorable weather and 
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a good seed crop so synchronized as to per- 
mit establishment of abundant reproduc- 
tion. In a single stand of sapling or pole 
size, every tree almost without exception 
will have begun its growth in the same 
year, so that age differences among trees 
can be ignored, 

Selected stands at each location were 
sampled, using 1/10-acre, 1/40-acre, and 
1/100-acre circular plots placed to include 
as many as possible of the typical conditions 
of stand age and stocking density in the 
area. Plot size was chosen to include ap- 
proximately 30 to 50 trees on each. On 
every plot, diameter at breast height of 
each tree was measured; and on most, all 
trees less than about 5 inches d.b.h. were 
felled and a circular cross section taken 
from each. Increment cores, at about 9 
inches above the ground, were taken from 
the few trees that remained uncut, On 4 
plots only a random sample, comprising 
about one-fourth of all the trees, was cut 
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or bored. Of the 14 stands included in this 
analysis, 3 originated in 1903, 7 in 1909, 
and 4 in 1914. Size of individual trees on 
the plots ranged from less than 0.5 inch to 
about 11 inches d.b.h. Stand density varied 
from the equivalent of 360 to 8500 stems 
per acre, 

The cores and stump sections were sur- 
faced, and ring counts made to determine 
the diameter inside bark of each tree at 
several earlier periods. Because of the need 
to determine tree ages precisely for this 
and other concurrent studies, sections were 
taken as close to the ground as possible. 
Unless otherwise stated, all diameters re- 
fer to measurements inside bark at or near 
ground level. A regression equation based 
on 94 trees over 4 inches d.b.h, indicated 
that diameter inside bark near ground level 
is equal to 0.93 d.b.h. The standard devia- 
tion of this regression coefficient is 0.014. 

Cross-sectional stem areas in square feet 
per acre were computed from the diameter 
data for each stand. The term “‘stem area” 
will be used to denote the cross-sectional 
area derived from diameter measurements 
at ground level. Stem areas were computed 
at the time of harvest and at each earlier 
ten-year interval in the life of the stand. 
Additional interpolated measurements of 
the 1909 stands made possible the com- 
parison of stem areas of all stands at a 
common age. 

Most studies of growth in forest stands 
have dealt only with trees of merchantable 
size, or with potential crop trees. Analysis 
of total production of plant material, how- 
ever, requires inclusion of every tree, in- 
cluding’ those that will not contribute to 
final harvest. ‘Total stem area per acre for 
each stand was divided by number of trees 
to find the average cross-sectional area, 
The latter was used as the primary de- 
pendent variable in fitting the data to equa- 
tions previously used for herbaceous plants. 

The past stem areas calculated here ob- 
viously include only those trees that were 
present in the stand in 1957; the areas for 
previous years include no allowance for 
trees that have since died. For the present 





analysis, however, this is of no importance, 
since this study is concerned with the past 
growth of those trees that make up the 
present stand. Natural mortality in these 
ponderosa pine stands is exceedingly low in 
the absence of fire, and is concentrated al- 
most wholly in the smallest size classes 


( Cooper 1960 ). 


Relations of Stem Area to Age and 
Stand Density 


Kira et al. (1953) planted soybeans at five 
different grades of stand density, and meas- 
ured the total weight per plant at six sep- 
arate times during the 119-day growing 
season. “They found that the relationship 
of mean plant weight to density at any 
given time could be represented by 


we = K (1) 


where w is mean weight per plant, p is 
stand density, and a and K are constants. 
They showed that this equation is appli- 
cable to a wide range of data, their own 
and that of other researchers. 

They found that the exponent a in- 
creased from zero at time of germination 
to approximately 1 at time of equilibrium, 
when dry weight per plant was inversely 
proportional to stand density. The value 
of a, from germination to equilibrium, in- 
creased as a linear function of the logarithm 
of time. Kira et al. considered that the ex- 
ponent a represented the degree of space 
utilization by an individual plant at a given 
stage of growth. They named this expo- 
nent the C-D (Competition-Density) Jn- 
dex, and termed the relationship symbolized 
in (1) the C-D Effect. The C-D rela- 
tionship assumes that final yield per acre at 
maturity is a constant regardless of plant 
density, an assumption generally supported 
by studies of pasture plants and by forest 
yield studies. 

The soybean experiments dealt with 
plants distributed uniformly on the plots. 
Further investigations (Kira et al. 1956) 
showed experimentally that the C-D effect 
holds equally well for plants non-uniformly 
distributed. The C-D relationship likewise 
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proved applicable to natural communities 
of weeds and to young stands of Japanese 
red pine (Pimus densiflora). In the same 
paper it was demonstrated that the identi- 
cal form of equation holds regardless of 
whether plant some indirect 
measure such as basal area is used in the 
calculation, and 


weight or 


whether 
mean plant weight is based on initial plant 
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density or on plant density at time of sam- 
pling. Of course, if natural mortality has 
reduced stand density in the time since 
planting, the numerical values of a and K 
in (1) will depend upon whether initial or 
present density is used to compute mean 
plant weight. Evidently, then, there is no 
barrier to the application of equation (1) 
to the 14 pine stands, where the individual 
trees in each stand are distributed irregu- 
larly, and the only density data available 
are those at time of final sampling. 

In order to fit observed data to the C-D 
equation, the latter can be expressed in the 
form 


log os log kK a log p (2) 


Mean stem area per tree plotted against 
stand density should give a straight line on 
log-log paper if the above equation holds; 
likewise, the fit of equation (2) can be 
tested by regression. Separate regressions 
of stem area on stand density, in the form 
of (2) were calculated for the 14 pine 
stands at each of 4+ ages. Interpolated 
measurements were used where necessary 
so that all the mean stem areas used in any 
equation were based on the same stand age. 

In all 4 equations, the regression coefh- 
cients were significant at the 1 percent 
level. Furthermore, as predicted by the 
C-D theory, the four calculated values of a 
fell on a straight line when plotted over 
the logarithm of time (Fig. 1), indicating 
an increasingly efficient utilization of space 
with the progress. of time. 

Upon further analysis, the same group 
of Japanese investigators who first proposed 
the C-D equation suggested another form 
of the relation between mean plant weight 
and stand density (Shinozaki and Kira 
1956). They pointed out that data on 
mean weight (or stem area) and stand 
density at any given time can be fitted to 
the equation 


1/w = Ap +B (3) 


where 4 and B are parameters whose val- 
ues change with time, and w and p are the 
same as before. 
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Four separate regression equations were 
calculated relating mean stem area per tree 
at four ages in the past to present stand 
density of the 14 pine stands. This time, 
(3) was used as the form of the regression 
(Fig. 2). Once again, the regression co- 
efficients were significant at the | percent 
level. When plotted against the reciprocal 
of time, the logarithms of both 4 and B 
fell very nearly on straight lines (Fig. 3). 


The Logistic Growth Equation 


In their study of several different growth 
curves, Shinozaki and Kira (1956) showed 
that field data conforming to the recipro- 
cal equation (3) can be fitted to the well- 
known logistic equation, Before attempting 
to do so, it is necessary to consider some of 
the properties of the logistic equation, and 
the assumptions underlying its derivation. 

The rate of increase, in absolute num- 
bers, of any population growing in unlim- 
ited space may be expressed as 


dN 





= rmN (4) 
dt 


where N is the size of the population at 


} 


time ¢, and r» is a parameter that may best 


be designated the imtrinsic rate of natural 
The rate of increase of the same 
population growing in limited space is 


iN 


picrease. 


(nr) cN)N. (5) 
it 

Here N is the density of the population, or 
number per unit area, while ¢ is the con- 
stant amount by which the intrinsic rate of 
natural increase is reduced by the addition 
of one more individual to the population 
The instan- 
taneous relative rate of increase, or rate of 


already present in unit space. 


increase per unit of population, is found by 
dividing both sides of the equation by N 
so that 


) 


1 dN 
, cN . (6) 
N at 


Equation (6) is the differential form of the 
logistic equation (Andrewartha and Birch 
1954). 

In the last 30 years the logistic equation 
has been used to describe a variety of 
growth phenomena, and much has been 
written about its applicability. The logistic 
theory of growth presupposes that every 
species has a definite intrinsic rate of nat- 
ural increase which is a constant in a given 
environment and which is as characteristic 
of the species as any other taxonomic or 
morphological feature (Andrewartha and 
Birch 1954). For any given species, there- 
fore, rm is considered to have a definite 
value which is dependent upon the particu- 
lar set of environmental conditions cur- 
rently operative. It is, in short, a theoreti- 
cal exponential rate at which a population 
might be expected to increase in a given 
environment if unimpeded by the influence 
of other organisms of the same kind. An- 
drewartha and Birch (1954) have pre- 
sented data from studies of insects and 
other animal populations to support the 
conclusion that there exists an attribute of 
animals corresponding to the intrinsic rate 
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of natural increase. 
The logistic equation is more familiar in 
its integrated form: 


N=————__ (7) 


1 + e*-"m! 


Here & is the constant of integration, the 
constant A is the asymptotic value toward 
which N tends as time ¢ approaches infin- 
ity, and ¢ is the base of natural logarithms. 
A biological interpretation of K suggests 
that it is the density at which the particu- 
lar space under consideration becomes “‘sat- 
urated” with the organisms involved, 

The simple logistic equation has been 
successfully fitted to the growth of animal 
populations. It apparently does not ade- 
quately describe the growth of plant com- 
munities, although MacKinney et «al. 
(1937) successfully applied it to loblolly 
pine stands. A_ basic assumption of the 
logistic theory is that rm remains constant 
with time. Shinozaki and Kira (1956) 
have suggested that the intrinsic rate of 
natural increase may actually change as 
plants grow older, so that rm = f(t). The 
simple logistic equation (7) can then be 
modified by substituting 


t 
[ f(T) dT (8) 
for rmt. This integral is essentially equiva- 
lent to the sum of all the infinitesimally 
changing values of rm from time zero to 
the time at which the value of the integral 
is desired. Shinozaki and Kira refer to the 
modified equation as a generalized logistic 
equation. 


Fitting Field Data to the 
Logistic Equation 


Shinozaki and Kira (1956) used a com- 
plicated graphical method to relate their 
data to the logistic equation. A simpler 
method is based on the well-known fact 
that the instantaneous relative rate of in- 
crease r in any population can be estimated 
by 
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N dt t 


1 dN IlnNe— 1lnN;, 
| an fl 


where Ni and Ne are the sizes of the popu- 
lation at any two known times and ¢ is the 
time interval between measurements. It is 
evident that the shorter the time interval ?#, 
the better is the estimate of r, provided that 
there are no great year to year variations 
in growth due to varying environments. 

The instantaneous relative rate of in- 
crease in stem area was calculated for the 
14 pine stands from equation (9), using 
mean stem area per tree as the estimate of 
population size. An estimated value of r 
was calculated for 3 ten-year periods in 
the life of each stand: 1928-1937, 1938- 
1947, and 1948-1957. The rate of increase 
so calculated is actually the integral of all 
the changing rates that prevailed over a 
ten-year span, The year-to-year variations 
in growth that became more noticeable as 
the period between measurements is short- 
ened precluded the use of a shorter interval 
than 10 years in the calculation of r. 

Tt was apparent that the calculated in- 
stantaneous relative rate of growth declined 
from decade to decade. According to the 
simple logistic theory, the only factor that 
should have contributed to the decrease in 
observed growth rate was the continual in- 
crease in tree biomass present per unit area. 
Plotting of the data indicated, however, 
that this was not a sufficient explanation. 
There was also a nearly linear relation be- 
tween instantaneous growth rate r and the 
reciprocal of stand age. This suggests that 
the intrinsic rate of natural increase, rm, in 
these stands is not a constant but may be 
made up of a component independent of 
time and a component that varies as the 
reciprocal of time, so that 


1 
rm = ro a ri. (10) 
t 


According to the logistic theory, the in- 
trinsic rate of natural increase is lowered 
an amount c by the addition of each new 
increment to the population. The proper 


nO =e ty 


measure of this increment is obvious in the 
case of animal populations — it is clearly 
just the number of animals added to unit 
space. In the plant community, though, 
the number of organisms is less important 
in its influence on growth rate than is some 
other measure of population size such as 
weight or volume of vegetation. 

The measure of space occupancy most 
commonly used by foresters is basal area, 
which is closely related to the stem area at 
ground level used in this study. Plotting 
the instantaneous relative rate of growth, 
r, against stem area per acre yields a curvi- 
linear relationship whose points show con- 
siderable scatter. This line can be rectified 
by assuming that the parameter ¢ varies 
inversely with time, 

A simpler hypothesis gives a better fit to 
the data, however. Rate of growth bears a 
linear relation to population size if the sum 
of all tree diameters per acre is used in- 
stead of stem area or basal area. This is 
logical when it is considered that stem di- 
ameter is almost directly proportional to the 
surface area of the cambium, whereas stem 
cross-sectional area is related to cambial 
area by a square root function. Lexen 
(1943) has suggested that total bole sur- 
face per acre is a better measure of stock- 
ing in forest stands than is basal area, since 
the bole area is actually the base upon 
which new growth is laid down. 

Significant quantities of carbohydrates 
that might otherwise be stored as new 
wood are utilized by the cambium in ,res- 
piration, The extent of this respiration is 
dependent upon the total surface area of 
cambium present. In stands of comparable 
height, one composed of many small trees 
will have a greater area of respiring cam- 
bial surface per acre than will one with the 
same basal area concentrated in a few large 
trees. The stand with the higher total rate 
of respiration might be expected to exhibit 
a slower growth rate. The direct relation- 
ship between stem diameter and surface 
area of cambium makes it appear that the 
sum of tree diameters per acre is a better 
measure of the influence of stand density in 


impeding growth rate than is basal area. 
This conclusion was confirmed in the pre- 
sent study. 

Substituting (10) in the differential 
form (6) of the logistic equation, 


l 
rT — fo Trim cS ( I] ) 
t 


where § is the sum of all tree diameters per 
acre and ¢ is stand age. This equation can 
be solved by least squares, and its applica- 
bility to the data at hand tested. Since the 
calculated values of r used in this solution 
are estimates of the mean _ instantaneous 
relative rate of stem growth over a ten- 
year period, the value of § should be an 
average of the diameter sums at the begin- 
ning and the end of each decade, Mathe- 
matical considerations dictate the choice of 
the geometric mean of S$; and Ss, or \/ S82, 
as the proper estimate of § to insert in equa- 
tion (11). The correct value of ¢ is the 
age midway between the time of measure- 
ment of consecutive values of S§. 
Based on data from the 14 pine stands, 
] 
r= 2.91 — 
t 


- 0.0000359 § . (12) 


Both regression coefficients were significant 
at the | percent level; the constant term 
ro was not significantly different from zero. 
The calculated regression accounted for 
about 83 percent of the variance in r 
around its mean value. 

Attempts were made to fit the data to 
several other equations, using total stem 
area per acre instead of the sum of tree di- 
ameters, and using the logarithm of time 
instead of its reciprocal. ‘The expression 
proposed in equation (11) gave the best fit 
and the smallest error variance. A still 
better fit could probably have been obtained 
if, instead of the sum of tree diameters, 
actual bole area per acre had been used as 
one of the independent variables. 

The intrinsic rate of natural increase, 
assumed in the simple logistic theory to be a 
constant, is shown to vary inversely with 
time, as suggested by Shinozaki and Kira 
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(1956). Equation (12) states that the 
instantaneous relative rate of stem growth, 
or rate of growth per square foot of stem 
area present, is inversely proportional to 
stand age, and that it decreases as the sum 
of tree diameters per acre increases. 

Substituting general variables for the 
specific coefficients found by regression, 
equation (12) can be rewritten in differen- 
tial form: 


{t t 


where, as before, zw 1s mean stem area per 
tree, S is the sum of all tree diameters per 
acre, ¢ is stand age, and ri and ¢ are con- 
stants. ‘his equation can be integrated, but 
the integration is laborious and the result- 
ing expression somewhat unwieldy, The 
graph of the integrated form of equation 
(13) is a more or less bell-shaped curve, 
with a rising sigmoid portion which reaches 
a maximum value of zw and then declines 
to zero. The rising limb of this curve ade- 
quately describes the growth of the trees in 
the ponderosa pine stands studied here. The 
data used in fitting the equation fell near 
the point of inflection of the sigmoid curve, 
making it relatively easy to fit tree areas to 
the assumed curve. The data could as well 
have been fitted to any of several other 
equations, 


Discussion 


Any acceptable growth equation must of 
course fit the data which it purports to de- 
scribe. Mere closeness of fit to field obser- 
vations is no guarantee of its validity, how- 
ever. Unless the parameters of a proposed 
growth equation can be shown to have some 
biological meaning, its application is little 
more than an exercise in mathematical 
curve fitting. The coefficients of a purely 
empirical equation have no definite physio- 
logical basis, and there is no way of know- 
ing how they might change with variations 
in stand conditions. 

The existence of real attributes of ani- 
mal populations corresponding to the pa- 


rameters of the ‘vgistic equation has often 
been suggested. In insect populations, and 
perhaps in other animal populations as well, 
there is evidently an intrinsic rate of nat- 
ural increase that is characteristic of the 
specie Ss. ‘This theoretical rate is never real- 
ized in the field, but it can be calculated 
from the birth rate, death rate, and age 
structure of the population (Andrewartha 
and Birch 1954). The observed popula- 
tion growth rate is the intrinsic rate of 
natural increase reduced by an amount 
proportional to the number of individuals 
present in unit space (Equation 6). 

The logistic equation adequately describes 
the stem area growth of the ponderosa pine 
stands studied here if it is assumed that the 
intrinsic rate of natural increase is inversely 
proportional to stand age rather than a con- 
stant. It is clearly possible to assign a real 
biological meaning to the coefficient ¢ in 
equation (12). It is the amount by which 
the instantaneous growth rate is reduced 
by the addition of each new increment of 
respiring tree surface to the stand. 

A rational interpretation of the intrinsic 
rate of natural increase itself is less easy, 
for unless some logical reason can be shown 
for the variation of this rate with time, it 
becomes merely a means of fitting an em- 
pirical equation to observed data. A possible 
explanation lies in the fact that as a tree 
grows older and increases in size, it re- 
quires a greater quantity of soil nutrients. 
The root system can effectively expand in 
only two dimensions while the above- 
ground portion is growing in three, and 
even lateral expansion is restricted by the 
roots of other trees. Thus the root system 
becomes progressively less able to supply 
all the tree’s needs as the tree grows older, 
and the growth rate declines due to the 
lowered root effectiveness per unit of tree 
biomass, 


Other Growth Equations 


Beverton and Holt’s (1957) comprehen- 
sive treatise on the growth and dynamics 
of fish populations deserves careful study by 
students of tree growth, for fishes and trees 
have much in common. Both continue 
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crowing more or less throughout their 
lives, and the growth of both is remarkably 
responsive to environmental influences. 
Beverton and Holt (1957:97) reviewed 
many of the published growth equations, 
and concluded that most were unsuitable 
for description of the growth of fishes. 
Either their constants had no_ biological 
meaning, the equations were based on 
faulty physiological assumptions, or they 
did not adequately fit the observed pattern 
of growth, The latter criticism was direct- 
ed against the simple logistic equation, 
which gives a symmetrical s-shaped curve 
that does not correspond to the actual 
crowth of either trees or fishes. ‘The point 
of inflection in both classes of organisms 
comes relatively early in life. Pearl and 
Reed (1923) introduced the necessary 
asymmetry into the simple logistic equation 
by adding a polynomial in ¢, but this modi- 
fication is wholly empirical, with less bio- 
logical meaning than equation (13). 


The equation that Beverton and Holt 
found most approprate to the growth of 
fishes was that of von Bertalanffy (1938). 
This is of the same general mathematical 
form as the logistic, but Beverton and Holt 
considered that its coefficients conformed 
more closely to the physiological charac- 
teristics of fish growth. Their exhaustive 
treatment of population dynamics is largely 
built around this equation. It may be 
equally applicable to the growth of trees. 

The von Bertalanffy equation requires 
more data for its satisfactory fitting than 
were available here. A graphical method 
for the handling of this equation was given 
by Beverton and Holt (1957:282). Pre- 
liminary plotting of the data from the 14 
pine stands indicates that the Bertalanffy 
equation may give a good fit. More re- 
cently, Parker and Larkin (1959) have 
criticized certain features of the Berta- 
lanffy equation as it applies to the growth 
of . shes. Their criticisms appear to be 
equally pertinent to the problems of tree 
growth. They have suggested a parabolic 
equation which, together with some of the 
methods of Beverton and Holt, may prove 


highly suitable for analysis of the growth 
of trees, 


Summary 


Two equations relating mean basal area 
per tree to stand density have been found 
applicable to even-aged stands of ponderosa 
pine in northern Arizona, ‘These equations 
proved equally applicable to past areas ob- 
tained from ring counts. Data on mean 
basal area per tree at several past ages have 
been fitted to the well-known logistic 
curve. The logistic equation gives a satis- 
factory fit to the grow th of these stands if 
it is assumed that the parameter commonly 
known as the intrinsic rate of natural in- 
crease is not a constant but varies inversely 
with time. 

The measure of population size most 
closely related to the instantaneous relative 
rate of growth in these stands was the sum 
of all tree diameters per acre rather than 
the conventional basal area per acre. This 
is logical when it is considered that diame- 
ter bears a more direct relationship to the 
surface area of the cambium than does 
basal area. 

Additional equations that have proved 
satisfatory for desc ription of the growth of 
fishes are suggested for analysis of the 
crowth of forest trees. 
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Effects of Photopertod and Source on 
Seedling Growth of Eastern Hemlock 


MANY PROVENANCE TEStTs of trees have 
been made in the field, but few indoors 
under controlled conditions of photoperiod 
and temperature. Seedlings of eastern hem- 
lock, T’suga canadensis (L.) Carr., have 
low light requirements and are therefore 
well suited to indoor experiments on the 
effects of both genetic and environmental 
factors. Such experiments have been car- 
ried on in Connecticut since 1952.1 The 
present paper shows how genetic differ- 
ences among seedlings from 30 sources of 
hemlock, combined with photoperiod and 
thermoperiod, account for relatively large 
differences in the duration and amount of 
seedling stem elongation. 


Review of Literature 


Racial variation. The importance of seed 
source for forest tree growth was first sug- 
gested by Duhamel du Monceau in France 
(Syrach Larsen 1949). But it is mainly in 
the last 50 years that the implications of 
genetic differences have been tested and 
applied by foresters. The earliest studies 
with native American species in North 
America were established by Zon with pon- 
derosa pine, Pinus ponderosa Laws., in 
1911 and with Douglas-fir, Pseudotsuga 
menziesti (Mirb.) Franco, in 1912. Both 
tests clearly showed the importance of seed 
origin (Munger and Morris 1936, Weid- 
man 1939, Munger 1947). 





1Background and results of this and other 
research, and acknowledgments of aid from 
many individuals, have been brought together 
in Bulletin 620 of the Connecticut Agricul- 
tural Experiment Station, New Haven, 1959. 


BY 
HANS NIENSTAEDT 
J. S. OLSON 


Since these early experiments, prove- 
nance tests have demonstrated the impor- 
tance of seed origin for a number of native 
species. Among these are: red pine, Pinus 
resinosa Ait., by Rudolf 1947, Hough 
1952a, b; jack pine, P. banksiana Lamb., 
by Schantz-Hansen and Jensen 1952; lob- 
lolly pine, P. taeda L., by Wakeley 1944, 
1954; lodgepole pine, P. contorta Doucl., 
Critchfield 1957; white spruce, Picea 
glauca (Moench) Voss, by Rudolf 1953; 
and sugar maple, Acer saccharum Marsh., 
by Kriebel 1957. All show very striking 
effects of seed origin, but, except for the 
work by Critchfield and Kriebel, they give 
little information regarding the relative 
plasticity of the different species and prac- 
tically none regarding the delineation of 
races within the species. 

Recently large-scale studies have been 
established. They include in the southern 
states the four important pines: slash, Pinus 
elliotti Engelm.; loblolly ; shortleaf, P. 
echinata Mill.; and longleaf, P. palustris 
Mill.; and in the Lake States region, jack 
pine and white spruce. Other species under 
investigation are red pine, white pine, Pinus 
strobus L., northern red oak, Quercus 

The senior author is Geneticist, in charge 
Northern Institute of Forest Genetics, a field 
unit of the Lake States Forest Experiment 
Station, Forest Service, U.S. Department of 
Agriculture. The junior author, while work- 
ing on this study, was Associate Forest Ecolo- 
gist at the Connecticut Agricultural Experi- 
ment Station; he is now Geobotanist at the 
Oak Ridge National Laboratory, Oak Ridge, 
Tenn. Manuscript received Dec. 21, 1959. 
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rubra L., various species in the genus Pop- 
ulus, and the more important western 
conifers. 

Field testing will always be a part of 
seed source studies, and the ultimate selec- 
tion of desirable sources must be based on 
field tests. Certain types of racial variation 
may, however, be studied at the seedling or 
even cotyledon stage by growing the plants 
under artificially controlled conditions 
(Schmidt 1951). Kramer (1936) noted 
that knowledge of the optimal photo- 
period should aid in predicting whether or 
not a given species is suitable for growth at 
any given latitude. His statement may be 
extended to include different races within 
the same species. 


Photoperiod and growth responses m forest 
trees. Photops riod—the durations of the 
daily exposure to light and dark—has a 
marked effect on the cessation of growth, 
leaf fall, and to some extent on the break- 
ing of dormancy (Wareing 1948, 1956). 
Different species may react very differently 
to a given dayleneth. Moshkov and Bog- 
danov, working independently of each oth- 
er, studied the photoperiod responses of sev- 
eral species from 17 different tree genera 
(Gevorkiantz and Roe 1935). A shorten- 
ing of dayleneth resulted in a decrease in 
the length of growing period, that is, pre- 
mature cessation of growth and more rapid 
hardening of the young shoot. Some of 
the southern species, if grown on natural 
dayleneth (Leningrad, maximum 20 
hours), continued growing into the fall 
and suffered severe frost injury. If grown 
on short daylength, they would ripen early 
and survive the winter without injury. 
Kramer (1936) studied 10 native North 
American species, including one conifer. 
On short-day treatments (8™% hours) all 
species except one grew less and became 
dormant sooner than the plants on normal 
daylength. When exposed to long days the 
responses differed among species—ranging 
from continued growth throughout the 
winter to little or no effect of increased 
daylength. Long-day treatments also re- 
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sulted in earlier breaking of dormancy 
than did short-day treatments. A light in- 
tensity of 2 foot-candles, or even less, was 
sufficient to bring about the type of re- 
sponses observed (Kramer 1937). 

The variation between species in response 
to different photoperiod depends on the 
time at which growth stops in nature 
(Downs and Borthwick 1956). Horse- 
chestnut, Aesculus hippocastanum  L., 
which stops growth in June and July, was 
the least sensitive to long days; elm, Ulmus 
americana L., which grows until frost, was 
the least sensitive to short days. However, 
long days will result in some extension, 
short days in a shortening of the period of 
elongation. 


Photoperiod and racial variation in forest 
trees. Photoperiod responses of trees of dif- 
ferent provenances of the same species often 
follow a definite clinal variation pattern 
(Langlet 1942-1943, Sylven 1940). At 
a given latitude local types of aspen and 
Scotch pine, Pinus sylvestris L., have their 
annual growth rhythm adjusted to the 
prevailing daylength regime. If trees are 
moved northward, the growth period 
lengthens, with corresponding late shoot 
ripening and possible injury from _ fall 
freezes. Provenances moved southward 
will cease growing early in the season even 
though other climatic factors are favorable 
for growth. 

Pauley and Perry (1954) have reported 
similar results from a study of black cot- 
tonwood, Populus trichocarpa Torr. & 
Gray, and eastern cottonwood, P. deltoides 
Bartr., and showed that the variation in 
the date for cessation of height growth is 
definitely clinal in nature, although a local 
diversity of photoperiodic responses occurs 
among clones independently of latitude of 
origin. ‘They advanced the hypothesis that 
adaptation to any particular length of 
growing season results from the natural 
selection of the genotypes which have a 
suitable response to the prevailing day- 
length regime of the latitude. 


Preliminary results with Douglas-fir by 


Irgens-Moller (1957) suggest that high- 
altitude genotypes may be more sensitive 
to photoperiod than low-altitude types be- 
cause of the greater survival value of the 
individuals that are well adapted to the 
more rigorous climate. 


Therm pe riodism nM trees. Kramer’s 
(1957) review indicates that surprising- 
ly little attention has been paid to tem- 
perature responses, particularly the  so- 
called “thermoperiodic” responses, i.e., re- 
sponses resulting from differences in day 
and night temperature. His results with 
loblolly pine and the following results with 
hemlock suggest that diurnal temperature 
changes and photoperiod may both have 
important effects on seasonal growth and 
geographic distribution, and that adapta- 
ions to both these factors differ within 


wide spre ad forest tre¢ specie s. 


I] . al — 
Photopertodic and thermopertodic responses 


in hemlock. Preliminary experiments with 
2-year-old seedlings from a Connecticut 
source (Olson and Nienstaedt 1953, 1957) 
showed striking differences in growth under 
8, 12, 16, 20, and 24 hours of illumination. 
As in Scotch pine (Wareing 1950), a com- 
bination of 20-hour days and 4-hour nights 
showed the best growth. Interruptions of 
the dark period by dim light also confirmed 
Wareing’s demonstration that tree growth, 
like many other photoperiodic responses, 
seems more aftected by duration of the 
dark period than of the light period; hence 
photoperiodic effects will be described in 
this paper by nightlength. The favorabk 
effects of high day-low night temperature 
combinations were also indicated in 10 
combinations of day and night temperature. 

The present study, with a more closely 
spaced series of nightlength treatments be- 
tween 8 and 12 hours and a wider range 
of temperature combinations, has confirmed 
these relations. It incorporated 30 seed 
sources and consistently showed differences 


between seed 


sources that appear to repre- 
sent genotypic variation resulting from 


adaptation to local climatic conditions. 


Materials and Methods 


Sources of material, Eastern hemlock cov- 
ers roughly a triangular area in eastern 
North America from Nova Scotia and 
Quebec (lat. 48° N.) to Minnesota and 
south to northern Alabama (34°30 N.), 
In the north it grows from sea level to an 


elevation of pe rhaps 1,600 to 1,800 feet 
and in the south from about 1,000 to about 
4,800 feet, with an occasional outlier at 
700 to 800 feet and at least one as low as 
Sti) feet. 


In 1953 seed was collected from 30 
locations. They spanned most of the lati- 
tudinal, longitudinal, and altitudinal rang: 
of hemlock. Several sources were from 
different altitudes in the same general lati- 
tude. Figure 1 shows the location and 
‘Table 1 gives pertinent data on each seed 
source. In the table the sources have been 
arranged in four groups according to the 
average length of the frost-free season, 
based on data from a_ nearby official 
weather station adjusted to allow for differ- 
neces in elevation, etc. This adjustment is 
only a rough one, but errors due to micro- 
climate or lack of nearby data are probably 
small relative to the whole span of growing 
season over hemlock’s range—from about 
100 to 200 days. 

Collections were represented by equal 
amounts of seed from 4 to 10 trees, except 
for two single tree coll ctions (Nos. 15 
and 18) near the limit of the range of the 
species. 


After extraction th seed was cleaned 


and empty seeds removed in an airstream. 
Samples of 100 seeds w weighed (Ta- 


ble 1) and tested for viability by cvermina- 
tion. 


Treatments. After about 90 days of strati- 
fication at 5° C, the seed was planted in 
pots. Germination was generally high and 
pots were thinned at random to secure a 
gxood spacing of 10 to 20 seedlings in each 
pot. 

When all the seedlings had reached the 
cotyledon stage and just begun to show the 
first new leaves (44 days after planting) 
they were moved to rooms or boxes where 
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Figure 1. The geographic range of eastern hemlock, showing source numbers of seed collections 
§ 5 / ’ S J 


used in indoor seedling experiments, and geographic trends in elongation averag 


° ? ? 
nightlength treatments. 


light and temperature could be controlled 
according to conditions given below. Plants 
were fairly unife*m at this time although 
some from the more northern sources, 
which had partially germinated during 
stratification, were slightly more advanced. 

All the seven photoperiodic treatments 
were carried out in one room, partitioned 
into compartments. Each compartment 
was furnished with banks of fluorescent 
daylight tubes. Light intensity was main- 
tained at approximately 450 foot-candles 
(4850 lux) at plant level. 

The plants were all exposed to 12 hours 
of the high-intensity light. To provide 12, 
11, 10, 9, 8, 6, and 4 hours of total dark- 
ness in each 24 hours, additional light was 
furnished from 15-watt incandescent bulbs, 
giving from 4 to 8 foot-candles (about 40 
to 80 lux) at plant level. The additional 
light was given following the full intensity 
illumination except in the 6-hour treatment 
where 6 hours of dim light preceded the 
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aged over several 


12 hours of full illumination. All these 
photoperiodic treatments were carried out 
at a temperature of 17° C (+2° C). 

In addition another series on 16 hours 
of light and 8 hours of darkness was main- 
tained at constant temperatures: 7, 17, 22, 
27, 32, and 37 degrees C. By moving some 
plants between the different temperature 
conditions, a third series of alternating (day 
and night) temperature combinations was 
obtained (Table 5). These plants were 
exposed to high temperature approximately 
8 hours and to low temperatures for 16 
hours. The “day-temperature” treatment 
was during the period of high-intensity 
illumination. All seed sources were in- 
cluded in the daylength and constant- 
temperature series, but only 15 in the alter- 
nating-temperature series. All sources on 
all treatments were replicated as two ran- 
domized blocks. 


Collection of data. The numbers of ac- 
tively growing and the numbers of plants 


wy 


e ws © 


showing buds were determined after the 
plants had been under treatment for 18, 
33, and 88 days days or a total of 62, 77, 
and 132 days after sowing the seed. At 
the same time the number of plants show- 
ing a second flush of shoot elongation was 
noted; a fully matured bud stage had usu- 
ally not been reached by the time this 
second elongation was resumed. In order 
to eliminate errors in measurements due 
to changing the soil line and to minimize 
seed weight effects, only elongation above 


TABLE 1. 


hemlock under controlled conditions. 


x 


Length Weight 


Monthly of of 
Seed North mean frost- 100 
source latitude Altitude temp. °F free seed, 
no, (degrees) (feet) Jan. July season grams 
153 35°45" 1,000 40 77 195 .342 
7 41°15’ 20 30) 73 195 357 
13 35°40" 1,300 40 76 192 .357 
17 35°40’ 1,200 41 75 188 396 
33 42°00’ 600 26 72 184 302 
12 35°40’ 2,300 38 74 184 374 
10 40°18" 100 33 75 181 .386 
20 38°15’ 1,900 32 69 164 341 
26 42°25’ 950 25 71 158 346 
183 35°45’ 4,800 34 68 156 311 
32 40°45’ 1,000 30 74 153 286 
99 44°40’ 100 2 65 150 .243 
8 41°45’ 500 26 71) 150 377 
35 46°47" 300 12 68 147 277 
4 44°00’ 400 20 70 139 294 
9 41°58’ 1,050 25 71 131 .329 
24 44°50)’ 120 18 68 130 362 
6 42°20’ 2,300 23 67 130 308 
29 46°15’ 1,000 9 70 128 382 
37 45°50’ 300 14 67 126 320 
38 45°50’ 600 8 67 122 272 
1 42°55’ 500 21 69 117 320 
a9 46°35’ 800 16 66 114 295 
34 46°20’ 1,500 7 65 111 306 
36 44°07’ 1,550 16 65 111 .264 
2 44°15’ 1,950 15 66 109 275 
30 46°00" 1,600 9 66 109 .247 
28 46°35’ 800 14 65 107 295 
3 43°20 1,950 17 64+ 106 -247 
5 44°00’ 2,100 15 63 98 281 


1After 62 days, all plants growing on nightlengths of 10 hours or less were still active. 


y 


the cotyledons were measured. Records 
were taken at intervals during the growth 
period of the plants, but only the total elon- 
gation at the end of the experiments is pre- 
sented here. The heights were recorded to 
the nearest millimeter. 


Results 


Duration of growth. Within 62 days after 
planting, some seedlings were already 
forming terminal buds under 11l- or 12- 
hour nightlengths (‘Table 1). Many more 
formed buds in the next 15 days, under all 


Source of seed and photoperiod as related to period of growth of 


Days after planting 
621 = 
Nightlength in hours 
12 11 12 11 10 9 8 6= + 


Percentage of plants growing 
100 100 0 20 5 39 45 0 100 
100 100 0 0 20 4 9] 22 67 
100 100 4 50 52 55 96 9? 97 
100 100 0 + 7 30 76 59 «100 
100 100 0 0 0 4 56 56 9? 
100 100 0 54 55 48 100 73 90 
100 100 0 33 4? 71 100 52 58 
100 100 0 0 0 0 71 56 28 
100 100 0 4 7 17 71 33 81 
96 06 0 0 0 18 83 39 9? 
100 100 0 0 0 0 29 0 47 
100 100 0 0 0 0 26 31 24 
94 97 0 12 0 0 35 5 39 
90 9? 0 0 0 0 24 0 7 
86 100 0 0 4 0 5 5 18 
100 100 0 0 0 8 30 39 74 
96 88 0 0 0 9 48 0 19 
94 100 0 0 0 0 60 31 63 
100 90 0 0 0 0 38 0 7 
86 95 0 0 0 0 7 0 40 
&8 75 0 0 0 0 0 0 
100 100 0 0 0 0 32 4+ 12 
&8& 75 0 0 0 0 12 4 44 
85 72 0 0 0 0 17 21 40 
86 90 0 0 0 0 26 0 12 
97 93 0 0 0 3 13 21 11 
84 66 0 0 7 4 52 0 19 
96 77 0 0 0 0 48 9 40 
96 100 0 0 0 0 32 11 0 
9? 82 0 0 0 0 28 0 27 





2In this treatment, the weak-light period preceded the period of high intensity light; evidently this made it less 
effective than in cases where it followed the high-light period. 


3Single tree collection. 
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nightlengths, 

Both nightlength and seed source (en- 
ironment and heredity) affected this timing 
of bud formation, which essentially stopped 
terminal elongation. For any given source, 
eedlings tend to form buds earlier under 
long nightleneths than under short night- 
lengths. But for nightleneth, 
seedlings of northern parentage tend to stop 
than of 


any given 


eTroOWINE earlic r those southern 


parentage. 


For example, after 77 days more than 


half the plants of low-elevation North 
Carolina sources (12 and 13, estimated 
frost-free seasons near 190) days) still had 


active stem tips under nightlengths as long 
A low 
showed a regular increase in the percentage 
of 
nightleneth to 100 under 4 hours. 


as 11 hours. ‘Tennessee source (17) 


active plants, from (© under 12-hour 
On the 
other hand, for most sources from regions 
with growing seasons of less than 164 days, 
all plants had stopped growth under night- 
lengths of 9 hours or greater; many plants 
had gone dormant, even under 8- and 4- 
hour nightleneths. 


Stem eli nati 
tion of growth were mainly responsible for 
the total height above the 
cotvledons, although differences in the rate 


of elongation 


n, The differences in dura- 


differences in 


also affected height growth. 
Figure 2A shows the gradation of decreas- 
ing total height with night- 
length. Comparisons with other parts of 
Figure 2 show the striking differences re- 
Under either 4- or 12 
nightlengths, plants of southern 
sources (2A, 2D) and sources from mid- 
dle latitudes near the Atlantic coast or 
Lake Michigan (2C), all of which have 


relatively long frost-free seasons, grew con- 


increasing 


lated to source. 
hour 





siderably more than plants from generally 
northern sources having short 
Ip 

(20). 


Similar ditferences 


growing 
seasons 
shown for seed- 
lings of sources in the same latitude, but 
different altitude. Figure 2E shows seed- 
lings under 4-hour nightlength, for the 
Great Smoky Mountains in the south and 
the Adirondack Mountains in the north, In 
each case the low-elevation (long-season ) 


are 


sources grew more than the high-elevation 
(short-season ) sources. 

Figure 1 shows the influences of both 
latitude and altitude of 
the quantitative trend in seedling growth, 
all Average 
elongation exceeded 14 mm only for low- 
elevation North Carolina sources. It ex- 
10 mm for coastal Connecticut, 
New Jersey, southern New York, southern 
Michigan, Ohio, and Tennessee sources, 
the exception of the high elevation 
source from the Great Smokies. Seedlings 
most of the 
throughout the main belt of hemlock for- 
ests from Nova Scotia to Minnesota, had 
mean elongation between 6 and 10 


source region on 


averaging over treatments. 


ceeded 


with 


of this source, and sources 


mm, 
while a few Canadian sources and that 
from Franconia Notch, N.H., averaged 
less than 6 mm. 

Since these seedlings were grown to- 


gether in essentially identical environments, 
it would seem natural to conclude that any 
differences that were statistically significant 
were due to genetic diversity within the 
species. This conclusion would seem all the 
in view of the natural 
geographic pattern outlined in Figure 1. 
The differences among 
sources might, however, be due merely to 


more convincing 
in elongation 


the increase in growth with seed weight, 
which is conspicuous in Figure 3A. The 
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apparent relation to frost-free season might 
conceivably be an illusion brought about by 
the tendency for seeds to grow larger in 
regions of longer seasons, as shown in Fig- 
ure 3B and Table 1. On the other hand, 
the relation of growth to frost-free season 
seems (Fig. 3C) to be even closer than its 
relation to seed weight, and this could not 
be explained completely by the correlation 
of seed weight and season. 

‘To clarify this point, the following sta- 
tistical approach was used. ‘Table 2 shows 
that the separate regressions of the log- 
arithm of mean growth to 10(0-seed weight, 
by; of weight to frost-free season, byr; 
and growth to season, bz,, are all statis- 
tically significant; in each case, variances 
accounted for by regressions, B’, are sig- 
nificantly greater than variances due to the 
scatter around these respective lines, as 
might be expected from observing Figures 
3A, 3B, and 3C respectively. 

To estimate effects of season independ- 
ent of seed weight, and of seed weight in- 
dependent of season, we need the partial 
regressions, bec ona ben F Dave Ss). Th 
the case of the 8-hour night treatment illus- 
trated in Figure 3C, we see that the first of 
these, 

by= bee y= [zx | +e12[zy ]=.004,27 
is indeed very much larger than the second, 
b= bey.a=c22[ zy | +e12[ 2x] =.000,773 


where the c’s are the elements of the in- 
verse matrix given in Table 3. Compared 
with the corresponding simple regressions, 
regression on seed weight now appears con- 
siderably less important than it did before. 
The sum of squares accounted for by both 
partial regressions, 
2 B°=B*,+ B22=h1 [zx] +b2[zy] 

700,358 

is only slightly greater than that accounted 
for by the simple regression on frost-free 


season alone, 
B* 2= br [zx ]=.678,949 

The difference between these, .021,409, is 
not statistically significant compared with 
the chance variance °=.013,446 due to 
scatter around the line, even in the case of 
the 8-hour night treatment where this term 
is relatively large. 

Thus partial regression analysis did, in- 
deed, reveal a distortion in the apparent 
importance of the two factors of season 
and seed weight, due to the fact that they 
are correlated with one another, But in- 
stead of unfairly exaggerating the effect of 
season because of a real effect of seed 
weight, it seems now that their correlation 
resulted in an exaggeration of the appar- 
ent influence of seed weight because of the 
real effect of season. This result thus 
eliminates the most obvious objection to the 
interpretation of the differences in seedling 


TABLE 2. Simple regressions relating elongation, seed weight, and frost-free 
season of seed source region.’ (Sixteen-hour day—S8-hour night photoperiod 


treatment used as example.) 


Frost-free Wt. of 100 seed 
season mg 
Mean x = 142.833 days 


Variance [x-] = 27708 [v-] = 57185 
Cross-products [vx] = 24346 
Slopes buy = .878;663 
Variance due to regression Beye = 21392 


Sums of squared deviations 


316.4 mg. 


Log mean 
epicotyl elongation 


a= 1.158,253 
[22] = 1.063403 


[zy] = .148,198 [=x] = 137.158 
bey = .001,551 bz = .004,950 
B2.y = .229,915 B22 = .678,949 


around regression [v2] — Bye = 35 [7] — B2zy = .833,488 [22] — B22. = .384,454 
Variance around regression s2, = 1278 2, = .029,767 s2, = .013,730 
Variance ratio (all signifi- 

cant at 1% level, F .01 

(1, 28) = 7.64) 16.738 vie 49.45 

1Brackets in table denote sums of deviations around means, e.g., 

[x2] = 2(x—x)2 = Tx2 — (2x)? N, fasi = Z(s—s) (x—x) = Zsx — (Zz) (Zx)/N. 
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54 
30 


TABLE 3. 


Regressions and variances for separate photoperiodic treatments. 


Nightlength* in hours 


Statistical computations! + 8 9 10 11 12 
Mean elongation 2 1.265 1.158 889 831 836 .766 
Simple regress 

On season bes 004,523 .004,950 004,301 .003,294 .002,777 .002,447 
On 100 seed weight bey .001,297 001,551 001,225 .001,047 .000,934 .000,712 
Partial regressions% 
On season by = beey .004,184 .004,271 .003,998 002,802 002,247 002,230 
On 100 seed weight bs = bey.z 000,385 .000,773 .000,346 .000,560 .000,603 .000,240 
Total scatter around mean [=~] 789,940 1.063,403 .7 86.688 538,917 504,741 371,130 
Due to season and weight 
=B2 = by [zx] + ba [ev] 572,087 .700,358 516,921 311,942 .226,788 .167,138 
Added by seed wt. B? — dz2[ =x] 005,315 .021,409 004,276 011,220 013,068 001,656 
Residual scatter = .008,069 .013,446 .009,99] 008,406 010,295 005,555 


1Brackets indicate sums of products or squared deviations around the means. 


“Nightlengths in these treatments varied only in duration 
12 hours of bright light common to all treatments. 

3Computed from the following elements of the inverse 
[xvy]/D = .000,024,548,419, coo [x-|/D 000,027,938, 


growth as being due primarily to some real 
genetic difference that is systematically re- 
lated to the climate of the region where 
the parents had grown. 

As would be expected from the differ- 
ences in timing of bud formation (Table 
1), there were differences in growth for 
seedlings of all sources, related to night- 
length. This is clearly shown by the shift 
in position of the regression lines in Figure 
3D. The main differences were in the 
general level of growth, but there seem to 
be differences in slope as well. 

Partial regression computations were 
made also for the other photoperiodic treat- 
ments and are shown in the columns of 
Table 3. In these, variances accounted for 
by seed weight, over and beyond that ac- 
counted for by frost-free season alone, were 
even smaller than in the preceding example 
for an 8-hour night. Asa result, the effect 
of seed weight is clearly not significant 
when tested in an analysis of variance (Ta- 
ble 4) for all six regression lines. 

This overall analysis also provides a test 
of whether the regressions of growth on 
season for different photoperiodic treat- 
ments, graphed in Figure 3D, represent 
essentially parallel lines or differ signifi- 
cantly in their slopes. The very large vari- 
ance accounted for by a single common 


matrix: ¢n [v-]/D 
7”) 


372, where D = —[xvy]* 


of a supplementary dim light follocing a period of 


4 .000,057,660,452, ec» = 


[x2] [y2] = 991754264. 


slope, i.e., the “main effect” of season, cer- 
tainly suggests a strong approach to paral- 
lelism. Nevertheless, the interaction of slope 
with treatment, based on the five degrees 
of freedom between treatments, is signifi- 
cant at the level of .025 and indicates a 
real heterogeneity among slopes. This is 
relevant to the later discussion of whether 
all sources can be considered as having the 
same “critical nightlength” for bud forma- 
tion, 

The preceding tests are based on the 
pooled estimate of scatter or chance varia- 
tion based on the 27 degrees of freedom 
for each of 6 treatments, for a total of 
162 degrees of freedom. Bartlett’s tests for 
heterogeneity among variances indicated 
that they were not significantly different 
for these 6 treatments. This conclusion 
suggests incidentally that the transforma- 
tion of the arithmetic mean elongations to 
logarithms was quite effective in eliminat- 
ing the large difference in variance that 
would have arisen if the analysis had been 
made on the untransformed means on an 
arithmetic scale. (Later analyses have 
reached essentially similar conclusions us- 
ing the means of the logarithms of the in- 
dividual data, instead of the logarithms of 
the arithmetic means which were used 
here. ) 
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The tendency for seedlings of northern At the very low temperature of 7 , the 








origin to make less growth than do those of whole process of growth was so slow that 
southern origin, which was clear for all the plants remained active a long time while 
nightlength treatments at 17°, was also making relatively little height growth, 
evident at various temperature treatments even in the case of long-season sources. In 
that were all maintained with an 8-hour this one case, there was a significant par- 
nightlength. However, the quantitative tial regression of growth on seed weight, 
differences in growth were generally not as beyond what could be acounted for by sea- 
great as at 17 for a variety of reasons. son alone. Probably the reasons for this 
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cotyledon stage, the unfavorable tempera- 
ture for photosynthesis of new food, and 
the « 


creater reliance on stored 
food, which would be supplied mainly by 


nseque ntly 


\lso, at relatively high temperatures, es- 
high night 
o ° ° ’ -O 
formation of buds was later than at 17 


pecially at temperatures, the 


constant temperature. While seedlings from 








short-season source regions tended to go 
dormant earlier than the others, the dif- 
ference was less striking than at 17°; evi- 
dently cool temperature as well as long 
nightlength enhances bud formation, and 
plants from short-season sources are more 
sensitive to this influence. On the other 
hand, very high temperatures like 32° or 
37° tend to suppress the amount of elonga- 
tion, even though plants continue to grow 
a longer time than at 17 While seed- 
lings of long-season origin grew taller than 
the others, this stunting tended to diminish 
their superiority in height growth. ‘Thus 
differences among sources in both the dura- 
tion and amount of elongation are smaller 
under these conditions, 

The spread in elongation between long- 
and short-season sources for varying com- 
binations of day and night temperature is 
summarized numerically by the regression 
coefficients (slopes) in Table 5A. These 
are greatest at 17° constant and 27° day- 
17° night temperatures, and decrease for 
both low and high temperatures, as just 
stated. 

The mean elongation, averaging over 
all sources, is given in Table 5B. This 
shows that daily temperature alternations 
between 17° and 27° are somewhat more 
favorable for elongation than either 17° 
or 27° constant temperatures, as is typical 
for plants showing a thermoperiodic re- 
sponse. Higher and lower temperature 
both tended to decrease total elongation. 


Discussion 

Eastern hemlock seedlings illustrate espe- 
cially well the influence of the different 
factors which regulate the growth and dor- 
mancy of trees. First, stem elongation even- 
tually ceases, even under the most favorable 
environmental conditions, probably because 
of some inherent physiological mechanism 
that becomes effective as the season pro- 
gresses. Secondly, nightlengths exceeding 
8 or 9 hours tend to accelerate bud forma- 
tion and hardening off of stems both in 
our experiments and in the field. Nursery 
seedlings of most sources commonly de- 
velop terminal buds around mid-August 
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(in southern Connecticut), at the time 
when natural nightlength is just increasing 
beyond 9 hours, including an hour of twi- 
light. Thirdly, cool temperature also ap- 
parently helps to hasten bud development. 
In addition, it is well known that many un- 
favorable environmental conditions such as 
drought or nutrient deficiency can further 
hasten the termination of growth. 

The differences in the duration of 
growth of the plants from different seed 
origins undoubtedly are the results of adap- 
tation to a combination of these environ- 
mental factors rather than one alone. Fur- 
thermore, an additional controlling factor 
is indicated by the fact that the regression 
of elongation on frost-free season is similar 
under all the perimental environments 
used. This factor presumably acts through 
the inherent physiological mechanism men- 
tioned, and is independent of the control by 
photoperiod or other external factors. It 
must be similar to the “endogenous mech- 
anism”? emphasized by Wareing (1956). 

On the other hand, this does not exclude 
the possible importance of differences 
among provenances in the so-called “criti- 
cal nightlength.” The term implies a night- 
length at which a small increment in length 
has a relatively greater effect on bud for- 
mation than changes at either longer or 
shorter nightlengths, without necessarily 
implying a very abrupt threshold, In Ta- 
ble 1, we noted that some seedlings of 
long-season origin tended to remain active 
on nightlengths exceeding 9 hours even 
after almost all of the short-season seed- 
lings had gone dormant under this light 
regime. The differences in slope of regres- 
sion lines (Fig. 3D, Tables 3, +) also sug- 
gest that there is relatively more spread in 
growth between 9-. 10-, 11-, and 12-hour 
nightlengths among long-season seedlings 
and less between 4- and 8-hour night- 
lengths than there is for short-season seed- 
lings. Further statistical analysis elsewhere 
(Olson et al. 1959) bears out the impli- 
cation of a slight but real difference among 
sources in “critical nightlength.” Night- 
lengths only slightly longer than 8 hours 
seem io have greatest effect in hastening 
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bud formation on seedlings of northern 
origin, while nigthtlengths approaching 9 
hours may be most effective for seedlings 
of southern origin. 

Such a small difference might seem in- 
consequential if it were not true that night- 
length changes rather rapidly in late sum- 
mer. Under a given changing photoperi- 
odic regime, a difference of a quarter or 
half an hour in the “critical” nightlength 
might make the short-season genotypes go 
dormant 5 or 10 days earlier than the long- 
season genotypes. Thus differences in re- 
sponse to photoperiod may be reinforcing 
whatever natural tendency exists for these 
genotypes to differ in the timing of bud 
formation and hardening off irrespective 
of photoperiodic control. But we cannot 
be sure of the relative importance of these 
two modes of control. 

In natural conditions, the additional ef- 
fect of cool temperature would also hasten 
bud formation in northern plantations as 
compared with southern ones.  Further- 
more, the greater “spread” among sources 
in time and amount of growth at 17° C as 
compared with higher temperatures, which 
is represented by the higher regression co- 
efficients in Table 5A, suggests that north- 
ern genotypes are more likely to go dor- 
mant at these moderately cool tempera- 
tures than southern ones even when photo- 
periods are favorable. This additional fac- 
tor further reinforces the genetic differ- 
ences related to inherent physiological con- 
trol and to photoperiodic regime. 

In view of the ever-present competition 
under field conditions, we can readily un- 
derstand why seedlings which tend to make 
the most growth would survive natural se- 
lection and pass this trait on to future gen- 
erations. Presumably hemlock, as well as 
other species (e.g. Bates 1930, Langlet 
1938, Sylven 1940), has by this selection 
become adapted to the most complete utili- 
zation of the long frost-free season of 
southern and coastal regions. 

Of course, any benefit from excess 
growth would be cancelled if plants re- 
mained active until the tender stem tissue 
was exposed to early autumn frost. Such 


damage has indeed been evident for several 
years in southern seedlings grown near 
New Haven, Conn., and dieback is begin- 
ning to counteract their initial advantage 
in height growth in many cases. Seedlings 
grown in Wisconsin from this nursery 
showed even more striking differences in 
frost injury. ‘The amount of damage had 
a highly significant positive correlation with 
the growing season at the point of origin 
and with the amount of elongation under 
the 9-hour night treatment of the present 
experiment (Nienstaedt 1958). In the 
north, early bud formation and hardening 
off has presumably been more advanta- 
geous than the extra height growth of 
southern seedlings. Thus any local popula- 
tion is subject to natural selection by op- 
posing tendencies, and it is the balance of 
these which conditions the normal duration 
of the growth period. 

The gradation in the timing and amount 
of growth shown in Figures 1 and 3 sug- 
gests that the pattern of genetic change is 
a cline and follows the gradual shift in cli- 
mate over the species range. In some places, 
as for example in the Great Smoky Moun- 
tains, large climatic differences occur over 
short distances. The hemlock population in 
some of these areas is almost continuous so 
that a free exchange of genes is possible 
along the intervening mountain slopes. 
Nevertheless, genetic differences may have 
been maintained by continuous natural se- 
lection which counteracts the dilution of 
genes in the population at one end of the 
climatic gradient by those from the other, 
even where short distances are involved. 

The sampling of the diversity within the 
species in this study hardly is intensive 
enough to prove definitely the existence of 
a single cline or to exclude the possibility 
of abrupt genetic discontinuities along some 
portions of the geographic gradient within 
the distribution. Hemlock in particular 
often occurs in groves which are isolated 
by physiography or disturbance history. 
Genetic discontinuities are apt to arise in 
such isolated populations. Squillace and 
Bingham (1958) have recently shown very 
localized genetic isolation of western white 
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pine in areas of abrupt contrasts in topo- 
graphic relief and in site quality. Quite 
possibly, intensive studies would find similar 
differences in hemlock, especially in those 
parts of its range in New England and 
certain higher portions of the Appalachian 
Mountains where it may occupy poor ridge 
sites as well as favorable coves. 

As long as a strong possibility remains 
that both gradual and relatively abrupt 
genetic changes occur within a species pop- 
ulation, it seems premature to draw a rigid 
contrast between clinal and “ecotypic”’ 
variation, ecotype implying discontinuity. 
Until much more is known about genetic 
patterns of a large number of character- 
istics of hemlock and of most other species 
that have been studied intensively, it may 
be advisable to use the term “ecotype” with- 
out a strict connotation of isolation. 

The correlation between seed weight 
and early growth of seedlings has received 
considerable attention in the literature and 
has been summarized by Hough (1952b). 
Often seedlings from large seed tend to be 
larger and heavier than those developing 
from light seed. However, the correlation 
between seed weight and the growing sea- 
son at the point of seed origin is also a 
(Langlet 1938, Bates 
1930, Baldwin 1942). Perhaps some of 


common one 


the differences that have been attributed 
to seed weight belong in small or large part 
to other genetic factors, which evolved 
through an adaptation to the 
frost-free pe riod, 


leneth of the 


Such confounding of variables can in 
1 
i 


part be overcome by the partial regression 


analysis used here. For example, this tech- 
nique applied to Hough’s red pine data 
showed a significant effect due to growing 
season (F 22.00). In his data, how- 
ever, unlike the hemlock data, there was 
also a significant regression due to seed 
weight (F = 5.1). 

Aside from the statistical, physiological, 
and genetic interpretations of these differ- 
ences, one can readily see the practical im- 
plications for those who are interested in 
planting hemlock for ornamental or for- 
estry purposes. Probably most commercial 


hemlock seed comes from relatively north- 
ern sources, in the main belt of hemlock- 
northern hardwood forests. In states with 
generally milder climate, as in Connecticut 
where more than half a million hemlocks 
are currently being grown in nurseries, 
stock from these sources is likely to grow 
more slowly in nurseries, and perhaps later 
in life, than seedlings of local origin. South- 
ern Appalachian seed, occasionally mar- 
keted in the north, is likely to produce 
plants that are more subject to frost injury 
and foliage discoloration than local stock. 

In mountainous regions and near large 
bodies of water where considerable  cli- 
matic differentials occur within relatively 
small areas, transfer of stock over small 
distances may result in considerable loss of 
production and damage from frost. This 
of course has been demonstrated for other 
tree species; the magnitude of such losses 
remains to be determined by field studies in 


each case, 


Summary 


Eastern hemlock seedlings from 30 sources 
were grown under controlled daylength 
and temperature conditions. They showed 
that: 

1. On any given regime of day- and 
nightlength, seedlings from a source region 
with a long frost-free season tend to form 
buds and stop elongation later than those 
from a region of short frost-free season. 
The 30 sources showed a distinct clinal 
variation in photoperiodic responses, but this 
does not exclude the possibility of relatively 
abrupt genetic changes either on a broad 
or local geographic scale. 

2. Seed weight was also dependent on 
the leneth of growing season at the point 
of origin. This correlation tended to exag- 
gerate the apparent regression of height on 
seed weight. A partial regression analysis, 
however, shows that the effect of seed 
weicht was generally negligible. 

3. Hemlock has a continuous gradation 
in response to varying nightlengths, shorter 
nightlengths permitting longer growth. 
The “critical” nightlengths showing the 
createst influence on growth seems to lie 





nt 


in the range of 8 to 9 hours, probably 
nearer 8 hours for seedlings of regions of 
short frost-free seasons and nearer 9 hours 
for seedlings of long-season sources, 

4. Height growth also is influenced by 
temperature in several ways. Constant 
high or low temperatures retard growth. 
The unfavorable effects of high day tem- 
peratures are counteracted by lower night 
temperatures, with 27° C day and 17° C 
night nearly ideal. 

5. Although field progeny tests of for- 
est tree ecotypes are necessary to assess the 
ultimate value of the various genotypes of 
trees, short-term studies with seedlings un- 
der controlled conditions may give valuable 
early indications of responses to environ- 
ment and the amount and type of genetic 
variation within a species; they can separate 
effects of heredity and environment and 
show their interactions in a way that is not 
possible in the field. 
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Suggestions to Contributors 


An article submitted for publication should 
be the best the writer is capable of producing, 
with al] statements, tables, quotations, names 
and formulas verified before submission. The 
careful author will have his typescript reviewed 
by colleagues. 

Articles should be typewritten, with 1%4- 
inch margins, on one side only of white bond 
paper, size 8% by 11 inches, or 8 by 10% 
inches. Carbon copies are not acceptable. All 
copy should be double-spaced. Pages should be 
numbered consecutively, preferably in the up- 
per right-hand corner. Only words intended 
to be set in italics should be underlined. 

Material such as equations that exceed 42 
spaces on the typewriter cannot be carried in a 
single line. Such material should either be set 
so that it can be broken into two or more lines. 
or should be presented in tabular form. 


Tables 


Each table should be typewritten on a sepa- 
rate sheet, given a title at the top, and should 
be numbered consecutively. Even if only one 
table is submitted, it should be designated 
Table 1. 

Footnotes used in tables should be designated 
by numerals, never by asterisks or other typo- 
graphical syrabols. These footnotes should be 
typewritten as part of the table. 


_ IMustrations and Figures 


An illustr-tion — whether a photograph, a 
line drawing, a map, or a graph — is desig- 
nated as a figure. A caption should be sub- 
mitted for each illustration, but the caption 
thould be typewritten on a separate sheet. 

All figures should be numbered consecutive- 
ly. If only one is submitted it should be 
designated Figure 1. On the margin or back of 
each illustration should be written lightly in 
soft pencil the number of the figure, the name 
of the author, and an abbreviated title. 

Photographs should be sharp and clear, print- 
ed on glossy paper. They should not be smaller 
than 4 by 5 inches. Sizes 5 by 7 inches and 


» 8 by 10 inches are preferable. Photographs 














should never be rolled or bent. Care should be 
taken in fastening photographs with paper clips 
which often make indentations that show up 
in reproduction. 

Drawings, maps, charts, and graphs are 
produced as line engravings. They should Sees at 
least twice as large as they are to appear when 





reproduced, and should be drawn in india ink 
on heavy white paper, or tracing cloth. 


Footnotes 


To indicate a footnote, place a superior fig- 
ure after the word that refers to the note. 
Consecutive numerals should be used, never 
asterisks. Separate the footnote from the text 
by running a line about one inch inward from 
the left margin of the type. 

Use footnotes to give credit to unpublished 
material and communications. If only a few 
references to literature are made (less than one 
per 1000 words), complete literature citations 
may be given in footnotes rather than in a 
separate listing. 


Literature Cited 


In Forest Science practice, references to 
literature citations are designated by the au- 
thor’s name and year of publication inserted in 
parentheses at the appropriate place in the text. 
If there are more than two authors, list only 
the senior author’s name in the text with the 
abbreviation ¢# al, Example: (Smith ef ai. 
1954). Only published references should be 
given in Literature Cited. Periodical abbrevi- 
ations should follow Guide to the use of For- 
estry Abstracts, Commonwealth Forestry Bu- 
reaa, Oxford, England, 1950. 


Nomenclature and Terminology 


When a species is first mentioned in a paper 
its common name may be immediately followed 
by its italicized scientific name in parentheses, 
but this latter name need not be repeated. 

As the authority for the exact spelling of 
tree names, Forest Science follows Check 
List of Native and Naturalized Trees of the 
United States (including Alaska), Agriculture 
Handbook 41, Forest Service, U. S. Dept. 
Agric., Washington, D. C., 1953. 

For the spelling of other plant names, con- 
tributors should follow Standardized Plant 
Nemes, 2d Edition, J. Horace McFarland 
Company, Harrisburg, Pa. 1942. 

Technical usage in forestry and allied fields 
follows Forestry Terminology, 3rd Edition, So- 
ciety of American Foresters, Washington 6, 
D. C. 1958. 

There are many style manuals available for 
the guidance of writers. One of the best is the 
Style Manual of the U. S. Government Print- 
ing Office, Washington 25, D. C. 1959. 

Webster's New International Dictionary is 
the accepted authority for general spelling. 
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